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Abstract 

4 

Near and Far-Field Acoustic Scattering through and from Two Dimensional 

Fluid-Fluid Rough Interfaces 


Recent experimental results [F. E. Boyle and N. P. Chotiros, J. Acoust Soc. Am. 91, 
2615-2619,1992; N. P. Chotiros, J. Acoust Soc. Am. 97, 199-214, 1995] reveal acoustic 
penetration from water into sandy sediments at grazing angles below the compressional 
critical angle in relation to the mean surface. These authors interpret their results to 
indicate the excitation of a Biot slow wave in the sediment. Another explanation is 
considered here. Modeling the ocean as a homogenous fluid and the sediment as a lossy 
homogenous fluid, computer simulations of these experiments based on analytical 
derivations in this work show that roughness of the water-sediment interface causes 
propagation of acoustical energy from water into the sediment at grazing angles below 
the compressional critical grazing angle; these simulations indicate that the experimental 
results can be explained in terms of diffraction of an ordinary longitudinal wave. These 
simulations use an analytical expression for the time-dependent mean square incoherent 
field scattered through (and from) a rough 2-D fluid-fluid interface that is derived in terms 
of the bistatic scattering cross section per unit area per unit solid area (differential cross 
section) of the rough interface. First-order perturbation theory is used to derive an 
expression for the differential cross section. The coherent field is calculated using the 
flat-surface result (zero-order perturbation theory) and compared to the coherent 
component of the second-order perturbation theory result. Effects of sound-speed 
gradients on the field scattered from the rough water-sediment interface are also shown 
using the first-order perturbation derivations. 
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CHAPTER 1 


Introduction 


1.1 Purpose 

A current area of considerable interest in ocean acoustics is the penetration 
of acoustical energy from a source in water into the seafloor (Williams, 
Satkowiak, and Bugler, 1989; Chotiros, 1995). Recent experimental results 
(Boyle and Chotiros, 1992; Chotiros, 1995) reveal acoustic penetration from 
water into sandy sediments at grazing angles below the compressional critical 
angle in relation to the mean surface. In addition, energy appears to propagate 
in the sediment at speeds slower than the compressional wave speed, in a 
direction concurring with Snell's law. These authors interpret their results to 
indicate the excitation of a Biot slow wave in the sediment. 

A primary goal of this work is to demonstrate that their experimental results 
can be explained in terms of diffraction of an ordinary longitudinal wave, rather 
than refraction of a slow wave. The effect of roughness on penetration of the 
seafloor is demonstrated by developing a rough surface scattering model that 
provides a solution to scattering into the rough seafloor. In order to provide a 
tractable model, simplifying assumptions about the ocean and seafloor are made, 
and an approach for solving the corresponding time-dependent rough surface 
scattering problem is taken that leads to simple analytical results. The ocean 
acoustic model used and the methods of solution are discussed in the following 
sections. Although the main problem is acoustic scattering through the rough 
fluid-fluid interface for an incident field in the water, analytical results derived in 
this work also include scattering from the rough fluid-fluid interface. Furthermore, 
these analytical results are general, and can be applied to other rough surface 
scattering problems. 
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1.2 Ocean acoustic model 

The propagation of sound in the ocean is very complicated (see for example, 
Urick, 1975; Brekhovskikh and Lysanov, 1991; Burdic, 1991). Sound speed is 
not constant, but varies as a function of temperature, salinity, and depth. In 
addition, acoustical energy propagating in the ocean is attenuated by absorption 
and by scattering due to inhomogeneities in the water. However, for the 
examples considered here, the propagation distances in the water are short 
enough so that the sound speed does not vary appreciably. Furthermore, at the 
frequencies of interest (300 Hz - 60 kHz) attenuation of acoustical energy over 
these short distances is also negligible. For this reason, the water is modeled as 
a lossless homogenous fluid. 

The ocean seafloor is in general an inhomogeneous porous, viscoelastic 
medium, possibly supporting the propagation of Biot slow waves (Biot, 1956a, 
1956b; Biot, 1962a, 1962b) and transverse, or shear, waves (e.g., Brekhovskikh 
and Lysanov, 1991) in addition to longitudinal, or compressional, waves. The 
seafloor has a roughness spectrum approximating that of a power law (Fox and 
Hayes, 1985; Goff and Jordan, 1989; Briggs, 1989). Acoustical energy incident 
on the seafloor will be scattered from this rough interface back into the water, 
and through this rough interface into the sediment. 

The main purpose of this work is to show the effect of roughness in the 
seafloor on scattering into the sediment. Although volume inhomogeneities in 
the sediment can be a significant source of scattering (Ivakin and Lysanov, 1981; 
Jackson, Winebrenner, and Ishimaru, 1986; Jackson and Briggs, 1992), volume 
scattering is beyond the scope of this work. However, the effect of vertical 
stratification of the sediment below the rough interface on the scattered field is 
discussed. Boyle and Chotiros (1992) and Chotiros (1995) neglect the effects of 
roughness in their interpretation of sediment penetration experiments, and 

v 

suggest that penetration of acoustical energy into a sandy sediment at shallow 
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grazing angles (below the compressional critical angle) can be explained by Biot 
wave propagation — provided nonstandard Biot parameters are used (Chotiros, 
1995). We show that their experimental results can be reproduced in a computer 
simulation that models the sandy sediment as a fluid and the interface between 
the water and the sediment as randomly rough. Ivakin and Jackson (in 
preparation) show that for the case of sandy sediments, the component of energy 
propagating as shear waves is small, and can therefore be neglected. 
Attenuation of acoustical energy is not negligible in the sediment for even very 
short distances, and the sediment is therefore modeled as a lossy fluid. 

Since it is not feasible to obtain an exact analytical description of even a 
relatively small portion of the seafloor, a probabilistic analysis is most practical. 
A statistical quantity of particular interest to rough surface scattering is the 
bistatic scattering cross section per unit area per unit solid angle (Ishimaru, 
1978b). For convenience, this quantity will be referred to as the differential cross 
section. It is useful to represent solutions for scattered acoustical energy in 
terms of differential cross section. However, exact analytical expressions for the 
differential cross section are not available. 

1.3 Rough interface scattering solution techniques 

Determining the field scattered from or through (penetrating) a rough 
interface given an arbitrary incident field is a problem with applications in several 
research areas, in addition to ocean acoustics (e.g., remote sensing, solid state 
physics, astronomy, optics). In the area of ocean acoustics, rough surface 
scattering is viewed as a contributor to propagation loss for sound propagating in 
shallow water (e.g., Kuperman and Ingenito, 1977) and a source of interfering 
reverberation in active sonar systems. In the acoustic detection of subsurface 
scatterers, the rough seafloor can be considered a source of interfering 
reverberation and signal loss (Nagy and Rose, 1993), and in medical imaging 
applications, rough interfaces in intervening tissue layers often lead to distorted 
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images (Berkhoff, Thijssen, and van den Berg, 1996). Although the rough 
interface scattering problem has been of interest for several decades, there is no 
exact solution method for the statistical problem. 

Solutions of the rough interface scattering problem are obtained by means of 
either numerical or approximate analytical methods. Although numerical methods 
can result in accurate calculations for the scattered field, statistical results require 
a numerically intensive “Monte Carlo” approach. In addition, numerical methods 
are very computer intensive for the three dimensional problem (Tsang, Chan, 
and Pak, 1993, 1994; Pak eta!., 1995). 

A classical analytical technique applied to acoustic scattering by 
Brekhovskikh (1952) and by Eckart (1953) is the Kirchhoff approximation, or 
tangent plane approximation. Although this approximation has a useful region of 
validity for Dirichlet boundary conditions (Thorsos, 1988a, 1988b, 1990; Thorsos 
and Jackson, 1989, 1991; Ishimaru and Chen, 1990), it is not as suitable for 
obtaining statistics of the scattered field for more general seafloor boundary 
conditions, such as fluid-fluid, where the reflection coefficient is a function of 
angle. However, as shown by Thorsos (1996b), the Kirchhoff approximation can 
be very accurate for scattering through a rough fluid-fluid interface near the 
Snell's law direction. 

More recent approximate methods include the phase perturbation method 
(Winebrenner, 1985; Winebrenner and Ishimaru, 1985a, 1985b), the small slope 
approximation (Voronovich, 1985), and the operator expansion method (Milder, 
1991). A summary of all of these approximate methods is given in Kaczkowski 
(1994). See also Thorsos and Broschat (1995) for a description of the small 
slope method. Although the phase perturbation method can provide an accurate 
approximation to the scattered field for Dirichlet and Neumann surfaces 
(Broschat et al, 1987; Broschat, Thorsos, and Ishimaru, 1989), it does not 
readily provide tractable solutions for the fluid-fluid seafloor scattering problem. 
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The operator expansion method is accurate in several scattering regimes for 
Dirichlet surfaces (Kaczkowski, 1994; Kaczkowski and Thorsos, 1994), but has 
not been developed for other boundary conditions and does not allow for formal 
averaging. A Monte Carlo approach is therefore required for obtaining field 
statistics. The small slope approximation can yield statistical solutions to the 
scattering problem with a fluid-solid boundary condition (Yang and Broschat, 
1994), using a separate supporting analytical technique in a region of known 
validity. The analytical technique commonly used with the small slope 
approximation is perturbation theory (Rayleigh, 1945). 

Perturbation theory, first used by Rice (1951) for randomly rough surfaces, is 
a straightforward technique for obtaining solutions to the scattering problem with 
complicated boundary conditions (Kuo, 1964; Dacol and Berman, 1988; 
Winebrenner, Farmer, and Joughin, 1995; Ivakin, 1994). The Rayleigh-Rice 
method (Rice, 1951) is a classical perturbative approach that uses the Rayleigh 
hypothesis as a starting point. In this approach, the incident and scattered fields 
are each expressed in terms of a Weyl, or a plane wave, representation. The 
Rayleigh hypothesis assumes that the scattered field can be represented 
everywhere in terms of a superposition of plane waves traveling away from the 
surface. Although this assumption is generally not valid close to the surface, 
Jackson, Winebrenner, and Ishimaru (1988) have established that the Rayleigh- 
Rice method and another more rigorous form of perturbation theory based on the 
extinction theorem (Nieto-Vesperinas and Garcia, 1981) yield identical solutions 
for the special case of a Dirichlet boundary condition through at least fifth order. 
They conjecture that the two perturbation methods yield identical results for all 
orders and all boundary conditions, and that the Rayleigh hypothesis is irrelevant 
to perturbation theory. 

The conditions for validity of the perturbation solution depend on the root 
mean square (RMS) height h of the rough surface as well as the wave number 
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k\ of the incident wave. It is generally accepted that for values of k\h much 
smaller than one, a first-order perturbation solution will usually, but not always, 
represent the scattered field with sufficient accuracy (Thorsos and Jackson, 
1989). For this reason, first-order perturbation theory is widely used for low 
frequency scattering from rough surfaces. Kuo (1964) has presented an 
expression for the differential cross section of a randomly rough interface 
separating two fluids based on first-order Rayleigh-Rice perturbation theory. This 
result has been frequently used to model the seafloor (Crowther, 1983; Jackson, 
Winebrenner, and Ishimaru, 1986; Stanic et al., 1988, 1989; Mourad and 
Jackson, 1989, 1993; Jackson and Briggs, 1992). Dacol and Berman (1988), 
Kuperman and Schmidt (1986), Liu, Schmidt, and Kuperman (1993), and Essen 
(1994) have applied perturbation theory to a rough shear-supporting sea-floor. 

Comparisons of first-order perturbation approximations with exact numerical 
results for the 2-D problem (1-D surface) have been made by Chen and Fung 
(1988) and Thorsos and Jackson (1989) for surfaces with a Gaussian roughness 
spectrum. These comparisons indicate that perturbation theory can be valid for 
values of k x h as high as 025. Thorsos (1990) shows by comparisons with exact 
integral equation results for scattering from a 1-D rough fluid-fluid interface, that 
first-order perturbation theory accurately determines the differential cross section 
for scattering from a "Pierson-Moskowitz" roughness spectrum with k\h=0A5. 
Using a truncated power-law roughness spectrum and sound-speed profile 
appropriate to the examples in this dissertation, Thorsos (1996a) shows by 
comparisons with exact integral equation results for scattering through a 1-D 
rough fluid-fluid interface, that first-order perturbation theory accurately 
determines the differential cross section over a wide range of scattering angles 
with k\h = 2. For scattering from a 1-D rough fluid-fluid interface, Thorsos shows 
that perturbation theory is valid for k\h- 035. Since a goal of this work is to 


7 


provide a tractable analytical model for seafloor scattering that is particularly 
accurate for scattering through a rough interface, perturbation theory is used. 

1.4 Overview 

Representing the ocean as a homogenous fluid, and the sediment as a lossy 
homogenous fluid separated by a rough interface, an analytical model based on 
perturbation theory is used to create computer simulations of the penetration 
experiments conducted by Boyle and Chotiros (1992) and Chotiros (1995) — 
even though these authors suggest that their results can be explained in terms of 
slow wave propagation in sediment that is a Biot medium. The general formalism 
for scattering of a pulse through a rough interface as well as from a rough 
interface is derived in chapter 2 of this dissertation. First-order fluid-fluid 
perturbation theory results used to approximate the rough interface scattering 
problem in the simulations are derived in chapter 3. In chapter 4, the results of 
the corresponding simulations of the experiments are given and compared with 
the experimental results of Boyle and Chotiros (1992) and Chotiros (1995). The 
simulations use the first-order perturbation model developed in chapters 2 and 3 
to calculate the incoherent scattered intensity at a receiver in the sediment, and 
use the flat surface, zero-order field to calculate the coherent intensity. Using the 
zero-order field to represent the coherent intensity is valid for small amounts of 
surface roughness; however, the surface roughness used in the ocean 
experiment simulations is k\h-\. The second-order component of the coherent 
intensity is derived in chapter 5, and compared with the zero-order component of 
the coherent intensity for the sediment roughness parameters used in the ocean 
experiment simulation. 

The simulation model also neglects gradients in sound speed in the 
sediment. In the approach taken in chapter 5, the medium below the interface is 
constrained to be a fluid only down to the lowest point of the rough interface. 
Below this point, the sediment can be an arbitrary vertically stratified medium, 
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characterized by its flat surface reflection coefficient. An application of the 
derivations in this chapter shows the effect of an upward refracting sound-speed 
gradient on the field scattered back into the water. This result is further 
generalized in chapter 6, where the sediment parameters below the lowest point 
on the rough interface are taken to be random. Although results showing the 
effect of gradients on scattering through a rough interface are not presented in 
chapter 6, this case can also be examined using the analytical approach in this 
chapter and chapter 5. Overviews of chapters 2 through 6 are presented in 
sections 1.4.1 through 1.4.5. 

1.4.1 Time-dependent mean square scattered field 

The total propagation time from source to receiver can be experimentally 
determined when the source signal is a pulse. For a known propagation path 
and known propagation speed, the distance of an object can be determined from 
the propagation time of a scattered pulse. Likewise, if the propagation speed is 
unknown, but the propagation path and locations of the source and receiver are 
known, the propagation speed can be determined from the pulse propagation 
time. In Boyle and Chotiros (1992), the propagation speed in the sediment is 
inferred from the arrival of signal energy. The main purpose of chapter 2 is to 
derive a near-field expression for the timG-dGpGndGnt mean square scattered 
field (average of square of field magnitude) due to a narrowband point source, 
including dispersion due to scattering and frequency dependence of wave speed. 
These results (Moe and Jackson 1996a, 1996b) are applied to the simulation of 
tank experiments (Boyle and Chotiros, 1992) and ocean experiments (Chotiros, 
1995); the results of these simulations are discussed in chapter 4. 

The time-dependent mean square scattered field can be expressed in terms 
of a two-dimensional integral over the two-frequency mutual coherence function, 
r((0!,(o 2 ) (see Ishimaru, 1978a, or Ishimaru ef a/., 1994a). Scattering of a pulse 
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from the ocean surface is addressed in Ziomek (1982a, 1982b), where r(coi,co 2 ) 
is presented in terms of a Kirchhoff representation. Ishimaru et al. (1994a) derive 
an analytical expression for r^coj,©^ for scattering from rough surfaces based 

on the Kirchhoff approximation, and Ishimaru et al. (1994b) use a second-order 
Kirchhoff approximation with shadowing corrections (Ishimaru and Chen, 1990, 
1991) to obtain an analytical expression for r(fo 1 ,to 2 ) that includes 
backscattering enhancement effects. In chapter 2, we obtain a general 
expression for r(co 1 ,o) 2 ) appropriate for rough interface scattering for a 

narrowband incident plane wave that does not depend on a particular scattering 
approximation. 

i 

In order to obtain simple analytic expressions for the time-dependent mean 
square scattered field, Ishimaru et al. (1994a, 1994b) assume the scattering 
channel is a WSSUS (wide-sense stationary uncorrelated scattering) channel 

(Ishimaru, 1978a) — the two-frequency mutual coherence function, r(co 1 ,o) 2 ) is 
only dependent on the frequency difference, coj = to 2 -fo 1 , and not on 
(©! +co 2 )/2. Because of frequency dependent attenuation in the sediment, we 
do not assume that r(©i,© 2 ) is only dependent on ©^. However, in chapter 2, 
we assume that the incoherent 7-matrix correlation function varies slowly over 
the source frequency range, and expand phase terms in a power series in 
frequency. Both the incoherent 7-matrix correlation function approximation and 
the power series expansion require that the source signal be narrowband. This 
approach is used to derive a general analytical expression for the time- 
dependent mean square incoherent field scattered from and through a 2-D fluid- 
fluid rough interface due to a narrowband incident plane wave which is 
expressed in terms of the second moment of the 7-matrix. At distances 
sufficiently far from the rough interface to neglect the evanescent waves, this 
expression for the time-dependent mean square incoherent field is given in terms 
of the differential cross section of the rough interface. This derivation leads to 
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the usual sonar equation in the limit as the narrowband signal approaches the 
CW (continuous wave) case. For the special case of a CW incident plane wave, 
this result becomes essentially exact for distances greater than about a 
wavelength from the interface. This result is surprising because the differential 
cross section is defined as a far-field entity (Ishimaru, 1978b; Winebrenner and 
Ishimaru, 1986). The expression for the time-dependent mean square 
incoherent field due to an incident plane wave is used to obtain the expression 
for the scattered field due to a narrowband point source heuristically. A simple 
analytical expression is obtained for the case of a Gaussian shaped source 
pulse, also in terms of the differential cross section. 

First-order perturbation calculations (first-order perturbation results are 
derived in chapter 3) for the case of a baseband Gaussian shaped source pulse 
illustrate narrowband pulse dispersion effects of the incoherent field due to 
forward scattering into a lossy sediment. For the case of incidence below the 
critical grazing angle, first-order perturbation computations also show that the 
incoherent field scattered through a rough surface can be much greater than the 
zero-order field transmitted below the corresponding flat surface depending on 
loss and receiver depth. These computations for the mean square incoherent 
field penetrating the rough interface are compared to the flat-surface result, for 
both plane wave and point sources. 

1.4.2 First-order perturbation theory 

Derivations for the first-order perturbation expressions used for the numerical 
computations in chapter 2 and in the computer simulations (chapter 4) are 
presented in chapter 3. Rayleigh-Rice perturbation theory is derived in chapter 3 
for the simple case of a rough interface separating two homogenous fluids. The 
water above the rough interface is modeled as a homogenous lossless fluid, and 
the sediment below the interface is modeled as a homogenous lossy fluid. The 
results of chapter 3 are equivalent to those published by Kuo (1964) and Mourad 
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and Jackson (1989), except that the field scattered into the second medium is 
also considered here. Although these results are derived for the two 
homogenous fluid case, this form of the first-order scattered field solution is 
shown to be more general in chapter 5. First-order perturbation theory solutions 
for this case are also used in the simulations of sediment penetration 
experiments in chapter 4. As in Moe and Jackson (1994c), a truncated or filtered 
power law is used for the roughness spectrum. This filter is shown to be 
equivalent to subtracting a moving average (Papoulis, 1984) from a 2-D surface 
with power-law roughness spectrum. 

1.4.3 Computer simulations of recent experiments 

A recent tank experiment (Boyle and Chotiros, 1992) and a recent ocean 
experiment (Chotiros, 1995) reveal acoustic penetration from water into sandy 
sediments at grazing angles below the compressional critical angle in relation to 
the mean surface. In their tank experiment (Boyle and Chotiros, 1992), a 
receiver buried in sand measures the intensity of a signal transmitted in water 
above the sand. The propagation time of the transmitted pulse is the same as 
that of slow refracted wave. Since only a single receiver is used in Boyle and 
Chotiros' tank experiments, no information about direction of the propagating 
energy in the sediment is available. In an attempt to determine both direction 
and speed of the energy propagating in the sediment, a buried array of receivers 
was used in an ocean experiment (Chotiros, 1995) that measures the intensity at 
an array of receivers in the sediment due to an acoustical source in the water. 
They use a simple algorithm to determine the speed and direction of the intensity 
front propagating in the sediment. The results of their experiments seem to 
suggest propagation of energy at speeds slower than the compressional speed 
of the sediment, in a direction concurring with Snell's law. 

Assuming a small amount of roughness at the water-sediment interface, and 
using the 3-D scattering model developed in chapter 2 with the first-order 
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perturbation theory results for 2-D fluid-fluid rough interfaces derived in 
chapter 3, we reproduce the results of the tank experiments, indicating that the 
acoustic penetration of the surface may be due to scattering (diffraction) from low 
levels of roughness (Moe et al., 1995) rather than slow-wave refraction. 
Assuming roughness parameters appropriate to the ocean experiment reported 
by Chotiros (1995), our goal is to reproduce his results with the model developed 
in this dissertation. An explanation by Thorsos (1995) of how diffraction can 
mimic a slow wave is given in chapter 4, along with discussions of simulation 
results (Moe et al., 1995). Some of the results presented in this chapter have 
been presented in Moe et al. (1995) and again in Jackson et al. (1996). 

1.4.4 Perturbation theory results for a fluid-fluid rough interface including 
gradients 

Moe and Jackson (1994a) follow the approach used by Winebrenner, 
Farmer, and Joughin (1995) in their treatment of electromagnetic scattering from 
rough dielectric surfaces, and constrain the sediment below the rough interface, 
down to the lowest point on the interface, to be a fluid. However, the medium 
below the lowest point on the interface is only constrained to be vertically 
stratified (e.g., sound speed can vary with depth, but not horizontal position) and 
is represented by its reflection coefficient. The derivation in chapter 5 takes this 
approach, allowing the sediment below the lowest point on the surface to be 
viscoelastic, or porous, supporting shear or Biot slow waves as well as gradients 
in physical properties. A first-order solution for scattering from a rough interface 
is also given in Ivakin (1994). Multiple layered rough interfaces in layered 
sediments are also considered in his presentation. 

First-order perturbation theory is valid for obtaining the incoherent field 
scattered into the sediment for the roughness parameters used in the simulations 
(Thorsos, 1996a), but generally, there is a need for higher order results to test 
validity and to correct lower orders. Also, sediments tend to be stratified, so the 
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assumption of spatial homogeneity of the second medium is sometimes 
questionable. Using Rayleigh-Rice perturbation theory, Nth-order recursion 
equations are derived in chapter 5 that provide a starting point for evaluating the 
Nth-order scattered field in terms of lower orders. The validity of using the zero- 
order result to represent the coherent intensity is shown in chapter 5, where an 
expression for the second-order coherent field is derived. 

1.4.5 Formal average with random sediment parameters 

In the simulation discussed in chapter 4, the sediment is taken to be a 
homogenous fluid. The perturbation derivations in chapter 5 allow for the effects 
of sound-speed gradients on the field scattered through and from the rough 
water-sediment interface. Strong gradients and discontinuities in density and 
sound speed on scales of several cm are often observed in shallow water 
environments (Briggs, Richardson, and Young, 1985; Richardson, 1986; Stanic 
etal., 1988, 1989; Jackson and Briggs, 1992; Lyons, Anderson, and Dwan, 1994; 
Jackson et a!., 1996), and gradients on scales of hundreds of meters are 
observed in deep ocean environments (Hamilton, 1979, 1980). Although the 
approach taken in chapter 5 includes the effects of layering and gradients in 
seafloor parameters below the rough interface, it does not allow for the inherent 
randomness of these parameters, which will sometimes "wash out" rapid angular 
variations in scattering caused by constructive and destructive interference. 
Perfect vertical stratification leads to interference effects which may not always 
appear. By allowing stratification parameters to be random, these interference 
effects can be averaged to reveal more robust features of the differential cross 
section. Note that this approach does not model volume scattering, since there 
is no horizontal dependence of sediment properties. 

In chapter 6, the first-order expression for the differential cross section 
derived in chapters 2 and 5 is generalized to account for randomness in the 
medium below the lowest point on the rough interface by expressing the 
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differential cross section in terms of moments of the stochastic flat surface 
reflection coefficient. Upward refraction or reflection due to large sound-speed 
gradients in the sediment results in increased scattering strength over the non¬ 
gradient case. This scattering gain occurs for various sediment types over a 
wide range of bistatic angles, but has a relative maximum when the incident and 
scattered grazing angles are equal (backscattering enhancement). 
Backscattering enhancement is observed and predicted in scattering from rough 
interfaces (e.g., Maradudin and Mendez, 1993; Ishimaruef a/., 1994b), but it is 
not manifest in first-order perturbation theory calculations that involve only 
scattering from a single rough interface with no additional specular or refracted 
path. Hanson and Zavorotny (1995) show that backscattering enhancement is 
apparent using first-order perturbation theory in a two-scale model for 
electromagnetic scattering from a rough ocean surface. In Hanson and 
Zavorotny, the energy path consists of a single scattering event and specular 
reflection; in the sound-speed gradient problem addressed in chapter 6, the 
energy path consists of a single scattering event and refraction. 





CHAPTER 2 


Near field scattering through and from a 2-D fluid-fluid rough 

interface 


2.1 Introduction to chapter 2 

In this chapter, a general analytical expression for the time-dependent mean 
square scattered field (average of square of incoherent field magnitude) due to 
an incident narrowband plane wave is derived and expressed in terms of the 
second order T-matrix statistics, as well as in terms of the bistatic scattering 
cross section per unit area per unit solid angle (differential cross section) of the 
rough interface. This novel approach yields an analytical expression that is not 
dependent on a particular scattering solution technique. In addition, the effect of 
loss in the medium below the interface is included in the result. Even though the 
differential cross section is typically considered a far-field entity (Ishimaru, 
1978a), the form of this result in terms of the differential cross section is valid for 
distances on the order of a wavelength from the interface. This rigorously 

derived result is used to derive the scattered field due to a narrowband point 

♦ 

source heuristically. This derivation leads to the usual sonar equation in the limit 
as the narrowband signal approaches the CW (continuous wave) case. The 
results of this chapter, along with the two-fluid perturbation theory results of 
chapter 3, are used in the simulation of sediment penetration experiments (Boyle 
and Chotiros, 1992; Chotiros, 1995), discussed in chapter 4. The results of this 
chapter are also presented in Moe and Jackson (1996a, 1996b). 

Figure 2.1 shows a 2-D slice of the 3-D ocean model. A time-dependent 
incident field in a homogenous fluid (water) above a rough water-sediment 
interface results in energy scattered from the interface back into the water and 
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through the rough interface into the sediment. The sediment is represented by a 
homogenous lossy fluid. The effect of depth-dependent sediment parameters on 
the time-dependent scattered field is not considered in this formalism. However, 
in the general perturbation theory derivation in chapter 5, as well as in the 
formally averaged sediment parameters derivation in chapter 6, the effect of 
sound-speed gradients on the intensity scattered back into the water for the CW 
case is discussed. General background for the scattering problem is given in the 

next section. 

The solution for a time-dependent scattered field requires obtaining the 
solution for the two-frequency mutual coherence function, r^G^) (Ishimaru, 
1978a). A general expression for r(a)i,co 2 ) f° r an incident plane wave is derived 

by assuming the 7-matrix correlation function is a slowly varying function of 
frequency. Dispersion effects due to frequency-dependent phase speed are 
treated by expanding phase terms in a power series about the center frequency. 
The validity of this expansion and the frequency-independence assumption both 
depend on the source signal bandwidth. The validity of the frequency- 
independence assumption for the 7-matrix correlation function is shown in 
chapter 3 to also depend on the distance from the interface. Using this 
approximate expression for r(coi,co 2 ), we obtain a near-field expression for the 
time-dependent mean square incoherent field (average of incoherent field 
magnitude squared) for narrowband source signals in terms of the second 
moment of the 7-matrix. Neglecting the evanescent field components in this 
result, we obtain a useful expression for the time-dependent mean square 
incoherent field in terms of the differential cross section. This result is 
remarkable in that a near field result is obtained in terms of what is generally 
considered to be a far field quantity. First-order perturbation results (derived in 
the following chapter) depict the contribution of the evanescent component to the 
incoherent scattered field for a CW plane wave source. The resulting profiles of 
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mean square incoherent field versus depth are compared with the flat-surface 
case, also illustrating the effect of roughness on the field penetrating a fluid-fluid 
interface. 

Using the narrowband time-dependent mean square incoherent field result 
for an incident plane wave, we heuristically obtain an analytical expression for 
the time-dependent mean square incoherent field due to a point source. First- 
order perturbation results illustrate dispersion of a pulse with Gaussian envelope 
due to forward scattering into a lossy sediment. These computations are carried 
out using a low-frequency cutoff for the bottom relief spectrum, and are 
compared to the flat-surface case. 

2.2 Scattering problem background 

The general scattering geometry is shown in Fig. 2.1. An arbitrary pressure 
field, \|/;(r), is incident on a 2-D rough interface separating the homogenous 

lossless fluid in medium 1 (water) from the lossy fluid in medium 2 (sediment). 
The zero-mean rough interface is defined by 

z = hf{ R), (2.1) 

where R is the transverse component of the three-dimensional position vector r, 
and h is the root mean square (RMS) height of the rough interface. The total 
field, \|/i( r )> above the rough interface is a sum of an arbitrary incident field, 
\|/-(r), and the resulting scattered pressure field, \| ff (r): 

¥ 1 (r) = Vi(r) + V|//(r). (2.2) 

The boundary conditions on the interface consist of the continuity of pressure, 

¥ 1 ( r )| z= /(R) = ¥2 ( r )l z =/ z /(R) ’ ( 2 ^ 


and the continuity of normal velocity, 
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:<24) 

where \|/ 2 is the pressure field in the sediment, h is the unit vector normal to the 
interface, pj is the density of the water, and p 2 is the density of the sediment The 
scattered field, yy(r), in the region above the highest point on the interface and 

the incident field are expressed in terms of the following (Weyl) plane wave 
expansions (Devaney and Sherman, 1973): 

\f f (r) = J d 2 K^ f { K)e'* lPl(K) V K ' R (2.5) 

and 

v .( r ) = J cP'K'Vi (K) e - '^ 1 ( K ) V K ' R , (2.6) 

where K is the two-dimensional transverse wave vector with magnitude K, k\ is 
the wave number in the medium above the interface, and 

p,(K) = Vl -K 2 lkl (2.7) 

The square root in Eq. (2.7) is chosen so that p] is either positive or positive 

imaginary. For the CW (continuous wave) case, a factor of e~ m is suppressed, 
but in general the above equations are the Fourier transforms of the time domain 

field, 

Vi(r) = Vi(r,co) = J<*\i/i(r,r)e l0W . (2.8) 

The field below the rough interface is expressed as a plane wave expansion of 
waves traveling in the negative z direction, 

\j/ 2 (r) = \|/ 2 -(r) = Jd 2 A'¥ 2 _(K) e ~ I * 2 ^ 2(K )V^‘ R , z<0, (2.9) 

where 


p 2 (K) = Vl -^ 2 /*2 


( 2 . 10 ) 
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and 

(HJ8)® (2.11) 

VC 1 

is the complex wave number in the lossy medium below the interface with 
frequency-independent loss parameter, 8. The speed ratio, v, is the ratio of the 
sound speed in the sediment, c 2 , to the sound speed in the water, q. Although 
the definition of k 2 in Eq. (2.11) violates causality (Aki and Richards, 1980; 
Wingham, 1985), we show these noncausal effects are negligible and consistent 
with the approximations made in deriving the analytical result. The square root in 
Eq. (2.10) is chosen such that p 2 is in the first quadrant of the complex plane. 
The field scattered back into medium 1 can be expressed in terms of the incident 


field and the T-matrix, 7j j: 

V f (K) = | d 2 K' 7] j (K, K ')% (K'). (2.12) 

Likewise, the downward component of the scattered field penetrating below the 
interface is expressed in terms of the T-matrix, T U - 

¥ 2 -(K) = |A / 7i 2 (K,K , )'l / t (K'). (2.13) 

The scattered field for the case of an incident plane wave is found by substituting 
'F ! (K) = 8(K-K i ) (2.14) 

into Eq. (2.12) and the result into Eq. (2.5), yielding 

y f (r) = | d 2 K 7] j (K, Kj) ( K ) z e ,KR . (2.15) 


The scattered field penetrating the interface for the incident plane wave case is 
found in the same way from Eqs. (2.9), (2.13), and (2.14): 

I|/ 2 _(r) = | rf 2 *rr 12 (K,K l -)e" tt 2 p2(K)z e*' R . 


(2.16) 
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For convenience, let \|/(r) represent both \|/y(r) and \K2-( r )- Define 

\j/(r) = \|/(r) — (2-17) 

where (\j/(r)) is the mean of the field averaged over the ensemble of rough 
interfaces. Similarly, define 


f(K,K i ) = r(K,K,)-(r(K,K i )) 


(2.18) 


to be the incoherent 7-matrix. Using the above definitions, \j/(r) is expressed as 


\j/(r) = JAf(K,K,)e ik 2 fo( K )z e ' KR . 


(2.19) 


2.3 General result for incident plane wave 

The purpose of this section is to derive an expression for the time-dependent 
mean square incoherent field at the receiver position, r, due to a narrowband 
plane wave source. Although the following derivation considers the more general 
case of scattering through a rough interface, it also applies to scattering from a 

rough interface. 

2.3.1 Time domain received signal 

A narrowband source signal, s(t), with units of pressure and with Fourier 
transform 

S(a) = J dts(,y°* (2-20) 

—oo 

can be represented as 

s(t ) = Rejw(r)e~ lC ° cf j, 


( 2 . 21 ) 
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where u{t ) is the corresponding complex baseband input signal with Fourier 
transform U( co), and (d c is the carrier frequency, or center frequency of the pulse 

(see, for example, Ishimaru, 1978a; Proakis, 1989). An incident plane wave 
pulse with direction denoted by the transverse unit vector, a,, and incident 

grazing angle, 6 ; , is expressed in terms of its 2-D Fourier transform as 

'F / (K,(o) = 5(co)6(K-K J )e'‘ to ^ sin0 ' /cl , O)>0, (2.22) 

where 


Kj = — cos0 ; d/ (2.23) 

<?i 

is the incident transverse wave vector, and z s is taken to be an arbitrary source 

height. Using Eqs. (2.9), (2.13), (2.17), (2.18), (2.20), (2.21), and (2.22) and 
chapter 5 of Ishimaru (1978a), we express the incoherent field xp bb (r,t) at a 

position r = (R,z) below the interface as 

Vbb M = ^ jd(aV((x>'p l2 {(0' + co c )e~ l(0 '', (2.24) 

where 


^ 12 (co' + co c ) = J d 2 K |f 12 (K,Kj,(o' + o) c )e , ‘ sine ' z ^ co ' +t0c)/ci 

xe -fKp 2 (K,co'+co c ) Z ((o'+co c )/c, £ ,kR 


(2.25) 


is the Fourier transfer of the impulse response of the incoherent field, generally 
known as the transfer function. The complex wave number ratio is given by 


K = &2 A 


(1 + /5) 

V 


(2.26) 


and 
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K/ = 


CO +C0, 


COS 9 ;&/ . 


(2.27) 


We have explicitly shown the dependence of (3 2 ( K ) on frequency and K. 


2.3.2 Time-dependent mean square incoherent field 

The time-dependent mean square incoherent field at position r can be 
expressed as (Ishimaru 1978a, chapter 4) 


I 2 (r,t) = (\y bb (r,0| 2 )- 


(2.28) 


Combining Eqs. (2.24) and (2.28) yields Ishimaru’s (1978a) Eq. (5-17): 


xr(co' + CO c ,CO" + Q) c )}, 


(2.29) 


where 


r(co' + C0 c , CO" + C0 c ) = (#12 K + C0 c )# 1*2 ((O" + 0) c )) 


(2.30) 


is defined by Ishimaru (1978a) as the two-frequency mutual coherence function 
Using the expression for #i 2 in Eq. (2.25), Eq. (2.30) is equivalent to 

r((o' + CD c ,C0 w + co c ) 

= j d 2 K" J d 2 K' (f 12 ( K',K;,co' + co^fj^K^K^co" + co c )) 

v0 ;K' R -;K" R isin B,z s co'/cj -/sin 0,z,co ”/c x 
XC C C 

_-*K|3 2 (K',CD'-HO c ) z(co'+co c )/c, ft «(KP 2 (K",C 0 "+C 0 c )) z(co"+CD c )/q ^ (2.31) 


where the transverse wave vectors for the incident field at the angular 
frequencies co" and co' are 
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K- / = ( ° c+CQ cos6,6; and K- = C ° c+(0 cosO/tx,-. 


q 


q 


(2.32) 


Define ^ ]2 (K',K'\K-,K",cfl c +co',(o c + ©") (see, for example, Zipfel and DeSanto, 
1972; Voronovich, 1995) such that 


(f 12 (K',K<,0)' + G> c ) 7 i* 2 ( + co c )) 

= eniK'X'XiXi&c +©> c + co")8 (K" - k'+k; - k;o . 


(2.33) 


We refer to Cn as the incoherent 7-matrix correlation function. In the above 
expression, the subscript is used to represent scattering into the sediment and 
the dependence of T n on co is included in the argument. Substituting Eq. (2.33) 
and 


K<'-K ;• = -——cos 0,a i = K d 

c\ 

into Eq. (2.31), and changing the integration variable K' to 
K = K' + Kj/2, 
yields 

r(o) , -i-co c ,co // + (o c ) 

= J^{^ 12 (K-K J /2,K + K J /2,K;,Kf,(0 c +(D',a) c +0)") 

x e -iK d R Q i sin 0 & (o'/q Q -i sin 0,-z, ©"/q 
x e -/Kp 2 (K-K d /2,G) c +co>co'/q 

x £ i(Kp 2 {K-K d / 2 ,© c +Co"))’ z(© c +©")/q 1 


(2.34) 


(2.35) 


(2.36) 


Equation (2.36) is now evaluated by making a few approximations. We assume 
is a slowly varying function of frequency, and for the case of a narrowband 

source signal, 
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e\ 2 (K - Krf /2, K + K d /2, k; , K£ Q) c + go", (O c + co') 

= ^ 12 (K,K,K,-,K/,ci) c ,( 0 c ) 

-^ 12 (K.Kf,co c ), (2.37) 

where the relation (Thorsos and Jackson, 1989; Berman, 1992) 

^ 2 (K,K,.0) c )8(K-K")»(3i 2 (K,K,.,m c )f 1 *2(K'',K i ,0) c )) (2.38) 

is a special case of Eq. (2.33), and K, is evaluated at the center frequency: 




OS 0yd/. 



(2.39) 


It is easy to show using first-order perturbation theory (see chapter 3) that 
the approximation in Eq. (2.37) is valid for K + |K^|/2 < , or 


K<k l - nf d cos0y/q, 


(2.40) 


where 2nf d = |co" - co'| can be taken to be the source signal bandwidth. For larger 
K, the phase of the quantity on the left side of Eq. (2.37) cannot be considered 
independent of frequency. The first-order 7-matrix representing scattering from 
the water into the sediment, T n , is derived in chapter 3 and given in Eqs. (3.41) 
and (3.45), and the first-order approximation of the incoherent T-matrix 
correlation function is given in Eq. (3.50). 

Define the vertical component of the wave vector in medium 2 as 

B 2 {K,(jd c + co',co c + co") = ((co c + co')/c,)kP 2 (K -K d /2,(a c + (o'). (2.41) 


Analytical expressions for dispersion in a waveguide can be found by expanding 
the propagation constant of the time-dependent field in a power series (for 
example, see Ishimaru, 1991, chap. 4). Similarly, with reference to Eq. (2.36), 
((o'jc\ )k(3 2 is expanded in a power series in co' and co". For the case of a 

narrowband signal, 
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#2 (/<f, CD c + (O', CO c + to") = #2c ( K ) + a), + 0^7 

’d 2 ft ,2 „ d 2 &2 

-^-co +2 ^ ■ - ■ * ■ 

3co /2 d(oda 


+• 


1 


, „ d^Bo „2 

“CD GO +- 4 -CO 

n -- //Z 


3 to' 


(2.42) 


where the derivatives are evaluated at co'= <o"= 0, and 
Blc(K) = k^ 1 P 2 (K,0) c ), 


(2.43) 


where the wave number, itj, in the above equation and all following equations is 
evaluated at the center frequency (k { = to c /q). Solving for the derivatives in Eq. 
(2.42) (see Appendix A), combining like powers of to' and co" in Eq. (2.36), and 
substituting the result into Eq. (2.29) yields 


h (r, t) = J d 2 K \e x 2 (K, K,, to c ) t 2klc Im(lcP2 (K ’® C ^ 

X j- j Jco" U* (to") e ,co, ^ _fl ~ h ) e i0) " 2p / 4 

x _L J </©' £/(©') ) e - w '“"V to ' 2 p/4 } 


(2.44) 


where 


R a,- „ z v . A , a,- -Kcos0,- 

r, =- L cos0,-+—sin0;+ z—- L Re 

C\ q qfcj 


A 


1 


\ 


V 


Kp 2 (K,to c ) 


(2.45) 


J 


and 


h - 


-z K /q 

p 2 (K,to c )' 


(2.46) 


The physical significance of these parameters will be discussed later. Note that 
t\ is real and r 2 is complex. The complex coefficient, 

P = {(2k 2 / E% c )(a r K C 0 c cos 0, / q 3 - ( tf/q f ) -i Im(l/52 c )(cos 0,- /q f f 
-i Im(l/ B% c )(a,- • K cos 0, /q f Jz 
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and the real coefficient, 

P = Imjjaj-Kcosej/Blc^^cAl 3 - a,- •Kcosejy^cfJJ-cos 2 ©,/^ 2 ^) -z, (2.48) 


can represent the contribution to pulse dispersion of the frequency-dependent 
phase term. This result is simplified by setting 

g\t) - 

xi fda>'U( 0 )')e- i “V i “ w 'V i<0 ' 2,> / 4 . (2-49) 

2 71 J 


The above definition leads to a simple expression for the time-dependent mean 
square incoherent field in medium 2: 




(2.50) 


In addition to being a function of the complex argument, t, g(t ) is a function of 
K, the incident field direction, and fluid-sediment parameters. Note that 


t d = t\+t 2 


(2.51) 


is also a complex function of K. For an arbitrarily small pulse bandwidth (CW), 
P, P, Im(r 2 ) are arbitrarily small, and g(t)->u(t). 


The definition of k in Eq. (2.25) violates causality (see Aki and Richards, 
1980, chap 5, and Wingham, 1985). When a corrected version of Eq. (8) from 
Wingham (1985) is used, causality will not be violated if 

K-25ln|(co c + co')/co c |/( 7 w) is substituted for K in (A2). Using this substitution, 

and following the derivation in Appendix A, results in causal versions of the 
parameters P,P,t\, and tj- Since 5«1 for ocean sediment, the difference 
between the causal representation of P,P-,t\, and t 2 and Eqs. (2.45) - (2.48) is 
negligible. 
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While Eq. (2.50) involves an integral over the transverse wave vector, it can 
be converted to an integral over the scattering interface through the change of 
variable, 

R' = R - vr d K/&!. (2.52) 


As shown in Fig. 2.2, R' is the location of a small scattering patch on the 
interface, and r d = \r d is the distance from this patch to the field point, r. 



R-R'| 2 +z 2 . 


(2.53) 


We constrain K to the range K<k l /\ in Eq. (2.52), that is, we neglect 
evanescent waves. This constraint is justified later when the contribution of 
evanescent waves on the mean square scattered incoherent field is illustrated 
using first-order perturbation computations for an incident plane wave. This 
constraint is also consistent with the frequency-independence assumption of the 
T-matrix correlation function. This assumption is invalid for arbitrarily large 
integration variable K in Eq. (2.50) (see Eq. 2.40). The Jacobian follows from 
Eq. (2.52). 



_ 

y d + z 2 x d y d 

{dR'J 

v 2 4 

x d y d x d ^ 


■l 


2 4 
v/j 



(2.54) 


where x d and y d are the transverse coordinates of r d . The geometric 
significance of these definitions is apparent when one assumes that the loss 
parameter is small (5 «1). Then 


*jKp 2 {K,<i> c )=f - a : 2 

= (itj sin 02 /v)f 1 + i§/:sin 2 02 j. 


(2.55) 


where 


sin 0 2 = -Jl-K 2 v 2 / k 2 = \z\/r d 


(2.56) 
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is the sine of the scattered grazing angle defined in Fig. 2.2. When terms of 
second order and higher in the loss parameter are dropped, the significance of q 

and t 2 becomes apparent, as 

tl = R' • a i cos 0,- fc\ + Zj sin 0,/q (2.57) 

is the time required for the incident plane wave front with direction 
(6c f cos9 ; ,sinQ,) to travel from the source reference point, (0,0,z s ), to the 

scattering patch, (R',0), and 


Re ( h) = r d/ c 2 ( 2 - 58 ) 

is the time required for a spherical wave scattered from the patch at (R',0) to 
travel through the sediment to the field point (R,z), z< 0. According to Eq. 
(2.55), the absorption exponent in medium 2 is given by 

Im(jt 1 KP 2 (^©c))N = “"N/ sine 2 = lm (h) r d- ( 2 - 59 ) 

Substituting Eqs. (2.54) and (2.59) into Eq. (2.50) results in 


I 2 (r,t) = jd 2 R' ^?M^ 12 (k,K f ,ffl c )e 2Im ^ 2)rrf g 2 (t~t d ). 
J v r d 


(2.60) 


Substituting the following quantities, 

<? 12 (a,d i ) * k ' Sm y^ e i2 {K,K i ,a c ), 

v 

a = K/K, and a t =KJKi 


(2.61) 

(2.62) 


into Eq. (2.60) yields 


/ 2 (r,0 = \d 2 R’ Cl2(& 4^e- 2Im t^)l r - r 'l g 2 (, - t d ) 

J r - rl 


(2.63) 
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The interpretation of o 12 as a differential cross section will be justified later in this 
section. Since K = * 1 (R-R , )/(vr J ), t\, t 2 , g, and therefore P and P are 
functions of |R-R'|. Note that the CW result follows from Eq. (2.63) by simply 
setting g(t-t d ) = 1. For this case of an arbitrarily narrowband source signal, the 

result in Eq. (2.50) for the mean square incoherent field is exact, and Eq. (2.63) 
approaches the exact result at distances at least a few wavelengths from the 
interface, where the evanescent waves contribute very little to the mean square 
incoherent field. 


Although the mean square incoherent field scattered into the water, /j (r,0, is 
similarly derived, the corresponding expression for I\(r,t ) is most easily obtained 
as a special case of Eq. (2.63) by substituting medium 1 parameters for medium 
2 parameters. For example, k\ is substituted for k 2 , |5j for k|3 2 , etc. For this 
special case, Eq. (2.63) simplifies to 


J r - r 


(2.64) 


where from Eq. (2.61) 

an («,&,) = kf sin 2 (0 2 \C\ i (K, K,, co c )| 2 . (2.65) 

In this case, 0 2 is also the scattered field grazing angle measured from the mean 
horizontal plane, \z\/r d =sin(0 2 ). Again, the quantity C\\ is found from the relation 
(Eq. 2.38) 

^ 1 (K.K„ro c )5(K-K") = (f„(K,K„m c )7i* 1 (K”,K i ,ffl c )). (2.66) 

There is no loss in this case: ^=0, and therefore, 






(2.67) 




30 



fifcW-*lPl(*-<»c)- < 269 > 

Note that t d is still given by Eq. (2.51), but the propagation time from the surface 
patch to the receiver/observation point above the interface is a real quantity. 

h =-p— v (2.70) 

As given in Voronovich (1994), the quantity Cufa,^) in Eq. (2.65) is actually the 
differential cross section defined in Ishimaru (1978b). 



(2.71) 


where \y, is a plane-wave field incident on a surface patch of area A A, and a,- is 
the unit vector in the direction of propagation of the incident field; vj/y is the 

scattered coherent field at a long distance |r-r '| 2 from the interface, in the 
direction denoted by the unit vector, a. Similarly, o J2 is the differential cross 
section relating the mean square incident field to the mean square scattered 
incoherent field in medium 2 (the sediment), defined as 


Oi 2 (&,&;) = 



2Im(k 2 )|r-r' 

» 


(2.72) 


where \j/ 2 _ is the incoherent field at the receiver in the sediment, a distance 
r-r'| from the surface patch, and the unit vector a represents the transverse 

direction of the scattered field in the lower medium. The result given in Eq. (2.63) 
(or Eq. 2.64) is remarkable, as the differential cross section is a far-field entity. 
Close inspection of Eq. (2.63) shows that it is the sonar equation when u(t) is 








substituted for g(t), in which the mean square scattered incoherent field is 
obtained by integrating the differential cross section over the interface with 
appropriate attenuation due to spreading and sediment loss. Equations (2.38) 
and (2.61) are a convenient way to find a 12 in theoretical developments; Eqs. 
(2.65) and (2.66) are a convenient way to find an- 


2.3.3 Dispersion of a Gaussian input pulse 


Two types of pulse dispersion occur in this problem. One type is simply due 
to scattering and is treated by the integral in Eq. (2.63) over the scattering 
surface. The other type of dispersion is due to the frequency dependence of the 
propagation constant and is embodied in Eq. (2.49). This subsection focuses on 
the latter form of dispersion. Consider the case of a pulse with a peak pressure, 
Ppeak- anc * Gaussian envelope, 

. <<)=/yak e" 2/ ' ? ■ < 2 - 73 > 

Its Fourier transform is given by 

U(a) = Ppeak^e - ^ 4 - (2-74) 


The parameter t s is chosen to vary the width of this input pulse, and therefore the 
pulse bandwidth and energy. Substituting Eq. (2.74) into Eq. (2.49) (see 
Appendix B) and integrating over cd' and co" yields 



(2.75) 


where 
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and 


P" = P * ~i(2P/q') 2 , (2.78) 

Substituting Eqs. (2.76) and (2.77) into Eq. (2.75), and simplifying, yields an 
expression equivalent to Eq. (2.75): 







exp 


-2Rej 


I'll 

V 

>s - >'f * 

2 o 

- (2f) J 


(2.79) 


As mentioned, Re(t d ) = Re(r 2 ) + h represents the propagation delay of the pulse. 

The imaginary part of t 2 is a result of modeling the sediment as lossy. Since 
higher frequency signals are attenuated more than lower frequency signals, the 
lower frequency components of a narrowband signal will be less attenuated than 
the center frequency. The imaginary component of t 2 , along with P, and the 
imaginary component of P compensate for typically excessive loss in the CW 
loss term — the exponential term in Eq. (2.50) or Eq. (2.63). In addition, P is a 
frequency dispersion term that results in pulse broadening of the signal below the 
rough interface. Note that for the zero loss case, the parameter P = 0. The 

effect of loss and the parameters P and P is seen in Fig. 2.3 for g 2 {t-t d ) from 

Eq. (2.79) evaluated at 


g‘ 


(t-td) 


t=Re{t d ) 


= 8 2 {~hi) • hi = lm ( f d ) = Im ( f 2 ) 


(2.80) 


and plotted as a function of scattered grazing angle. In Fig. 2.3, the pulse length 
parameter, t s , is equal to two periods of the center, or carrier, frequency for both 

cases shown. When the loss and dispersion are small, \g(-t 2 i )| approaches 

u(0) = 1. The smaller the scattered grazing angle, the greater the dispersion, and 
the smaller one would expect the pulse peak amplitude. However, the loss of the 

CW signal represented by the term e ~ 2 Im (*2)l r - r l is greater than the loss of the 
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narrowband signal, and the propagation distance becomes large for small 
scattered grazing angles (z is fixed). Even though there is dispersion, the loss 

for the narrowband case is significantly less than the CW case, and |g('-f 2 ,-)|> 1 • 
In the case of Fig. 2.3, with depth of 0.3 m, |g(-r 2 ,-)| is greater than one for very 
small grazing angles since the propagation distance becomes very large. An 
example with zero loss always has P"= P* and f 2 ; = 0, and it is always true that 

k(o)|<i. 

2.4 Applications and calculations 

Although the theoretical results in this dissertation are general and include 
scattering back into the water, the examples deal with scattering through a rough 
fluid-fluid interface (z< 0). Both plane-wave and point sources are considered, 
and the mean square scattered incoherent field is computed for both the flat 
surface case and the rough surface case. 

In the examples of the following section, the mean square scattered field in 
medium 2 is divided into coherent and incoherent parts. The incoherent 
component is treated using first-order perturbation theory, and the coherent 
component is treated in zeroth order — the flat surface solution is used. There is 
a slight inconsistency in this approach, in that the mean square incoherent field is 
proportional to the second power of k x h, while the square of the coherent field 
magnitude is only computed to an accuracy of zeroth order. As Rice (1951) has 
shown, the square of the coherent field magnitude to second order in kfi is 
obtained by subtracting the power carried by the mean square incoherent field. 
For small k\h, this correction is necessarily slight. In chapter 5, the second-order 
perturbation scattered field and second-order coherent field are derived. 



34 


2.4.1 Continuous plane wave source 

The contribution of the evanescent waves to the field scattering through a 
rough interface is shown in the following CW incident plane wave examples. In 
Fig. 2.4a, the first-order field penetrating through a rough interface, with 
roughness spectrum given in Eq. (3.66) and ^^ = 0.2545, is calculated using Eqs. 

(3.70) in Eq. (2.50) and setting g(t-t d ) = IPa. The solid line in the plot 

represents the first-order result including evanescent waves, and the dashed line 
is the first-order result excluding evanescent waves, which is the result in Eq. 
(2.63). Since the grazing angle in this example is below critical, the zero-order 
field (dotted line) is evanescent. Near the interface, the mean square incoherent 
field decays with increasing depth at a rate comparable to the square of the 
magnitude of the zero-order field, but decays at a rate dictated by the sediment 
loss farther from the interface. Including the evanescent waves (K > ^/v) does 
not make a significant difference at a distance greater than about a wavelength 
(7.5 cm) below the interface. The sediment parameters used in the plot in Fig. 
2.4b are the same as Fig. 2.4a, except the roughness parameters are chosen 
such that k\h = 1.0. Perturbation theory is shown valid for k x h = 1 by Thorsos 
(1996a) for scattering through a rough interface with truncated power-law 
roughness spectrum. In this case, including the evanescent waves of the integral 
(K> ki/v) does not make a significant difference at a distance greater than about 

12cm below the interface. 

These examples show that including the evanescent waves in the calculation 
of the first-order mean square incoherent field is only necessary close to the 
interface, and the contribution of evanescent waves to the mean square 
incoherent field is dependent on the interface roughness. It is important to note 
that in the region below the interface where the incoherent evanescent 
contribution is significant, the zero-order coherent contribution to the mean 
square field is typically much greater. 
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2.4.2 Mean square incoherent field due to a point source 


Starting with the plane wave result in Eq. (2.63) and expressing the 2-D 
integral in terms of a summation over the interface yields 


l 2 {r,t) = AA]r 

n 



hn *2 n) 


e -2Im (k 2 )r 2n 


where 


(2.81) 



n 


(2.82) 


is the distance from surface patch n at position r' n of size A A to the observed 
field point, r, and 



R-R' 

R-R' 


(2.83) 


The propagation time from this surface patch to the receiver, Re(r 2n ), is also 
given by 


Re(r 2n ) = r 2n /c 2 . 


(2.84) 


As the area of each surface patch, AA, approaches zero, the summation in Eq. 
(2.81) approaches the exact integral in Eq. (2.63). 

An expression for the first-order mean square field below the interface due to 
a point source follows heuristically from Eq. (2.81). Here, the surface patch (area 
AA) is chosen sufficiently small in relation to the distance, r ln , from the point 

source to surface patch n, to ensure that the field incident on the interface is a 
plane wave. However, the size of this patch cannot be chosen arbitrarily small, 
and the distance between the source and interface is therefore constrained to be 
greater than some minimum value. Limitations of this approximation remain an 
issue. The scattered or penetrating mean square incoherent field is given by 
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/ 2 (r,0£AAX — J2 —VsT-^e 2 " . 


n 


2 2 
r ln r 2n 


(2.85) 


where 




( 2 . 86 ) 


R 5 is the transverse coordinate vector of the source, g(t) is the pressure at 1 m 
from the source, and r 0 = 1 m. The total propagation time is given by the real part 
of 



(2.87) 


with 





( 2 . 88 ) 


When the approximation g{t) = u(t) is made, Eq. (2.85) is in the form of a 
convolution with reduced computation time. 

Using the expression for the first-order differential cross section given in Eq. 
(3.77), the approximate mean square incoherent field pulse in lossy sediment 
below a rough interface due to a Gaussian narrowband point source above the 
interface is found for two cases. In Fig. 2.5a, the time-dependent mean square 
incoherent field at a position below a rough interface in a lossy medium due to a 

point source is found using Eq. (2.79) for g 2 (t) in Eq. (2.85), and is compared to 
the result obtained setting P and P to zero — g(t) = u(t). Although the resulting 
pulse shapes are close in magnitude, as well as peak arrival time, using g(t) 
does result in an incoherent pulse that has a smaller peak magnitude. In this 
example, the parameter t s in Eq. (2.73) is equal to two cycles of the center 

frequency of f c = 20 kHz. This frequency results in k x h = 0.25, which is within the 
region of accuracy for first-order perturbation theory (Thorsos and Jackson, 
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1989; Thorsos, 1990; Thorsos, 1996a). In Fig. 2.5b, the point source is farther 
from the interface, but the incident angle is the same. The incident grazing angle 
is defined here to be the angle between the mean surface and the line containing 
the source and receiver points. 


2.4.3 Zero-order calculations 


Setting 

V 2 _(K) = T(K)^(K) 


(2.89) 


and 



ir 0 U c -iK.R, c ft,Pi(K)z, 
2tcA:iP 1 (K) 


(2.90) 


in Eq. (2.9) results in an expression for the field penetrating a flat fluid-fluid 
interface due to a point source at (R s ,z s ) with unit pressure magnitude at lm 

from the source. The quantity U has unit magnitude and dimensions of pressure 
x time. The resulting integral simplifies to 


\j/ 2 _(r) = fJ dK e f k* lPl(K) ^ 2p2(K) ) Jq(KR), 


(2.91) 


where J 0 is the Bessel function of order zero. For comparison purposes, the 
square of the time-dependent zero-order, or flat-surface coherent field 
magnitude, 


I(r,t) = 


J dc o U((a)r 0 e~ i<£>t J dK ^ ^ e^ lPl (K) ^ 2p2(K)) J 0 (KR) 


0 


(2.92) 


is plotted in Fig. 2.6a along with the time-dependent first-order mean square 
incoherent field using the same parameters as in Fig 2.5a. In this example, the 
incident field is below the critical angle. Although the zero-order component can 
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be approximated from Eq. (2.91) using the method of stationary phase to solve 
the resulting integral (see for example Brekhovskikh, 1980; Ishimaru, 1991) when 
the incident field is above the critical angle, a numerical method is more suitable 
when the incident grazing angle is close to critical (Westwood, 1989). Figure 2.6a 
shows an “exact” solution of the square of the flat-surface transmitted coherent 
field magnitude, along with the first-order mean square incoherent field from Fig. 
2.5a. For this example, the zero-order coherent pulse arrives earlier and with 
greater magnitude than the first-order incoherent pulse. In Fig. 2.6b, the source 
is moved farther from the surface than the previous two examples. Even though 
the sediment parameters, receiver point, and incident angle (the angle between 
the mean surface and the line containing the source and receiver points) are 
unchanged from Fig. 2.6a, the incoherent first-order pulse is greater in magnitude 
than the zero-order pulse. As the source is moved further from the surface, the 
stationary phase path for the zero-order pulse is longer, and the zero-order pulse 
amplitude is more attenuated in the lossy sediment. Since Fig. 2.4 is a plane 
wave CW version of this same example, one would expect that the first-order 
mean square incoherent field would be greater than the zero-order field squared 
for the source point far from the rough interface. 

2.5 Summary 

A general analytic expression for the time-dependent narrowband mean 
square incoherent field scattered from and through a fluid-fluid rough interface is 
expressed in terms of the second moment of the T-matrix, and in terms of the 
differential cross section. This result is valid for narrowband source signals when 
the T-matrix correlation function is a slowly varying function of frequency over the 
source signal frequencies. Using the exact CW incident plane wave result, along 
with first-order perturbation theory, and sediment parameters characteristic of a 
sandy seafloor, we show that the evanescent component of the scattered 
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incoherent field in the sediment becomes insignificant about a wavelength from 
the surface, depending on the surface roughness. These examples show that 
the expression in terms of the differential cross section is valid except within one 
or two (depending on the roughness spectrum) wavelengths of the surface. This 
result is surprising since the differential cross section is a far field entity. 

The above formalism applied to the scattering through a rough surface 
shows the effect of roughness on penetration through a surface. When the 
grazing angle of a plane wave incident field is below the critical angle in relation 
to the mean surface, the zero-order component of the transmitted field is 
evanescent, and does not penetrate deeply into the seafloor. Higher order 
components contain downward traveling waves, which can increase the depth of 
penetration of sound relative to the flat-surface case (Moe and Jackson, 1994b; 
Moe etal., 1995). 




Figure 2.1 : Scattering problem geometry for rough interface separating i 
lossless fluid in medium 1 (z > hf(R )) from a lossy fluid in medium 2 (z < hf{ R)) 

This diagram can be viewed as a slice through a two-dimensional surface. 















































































































































































































































































Figure 2.2 : Diagram describing transverse position variables 
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scattered grazing angle (degrees) 


Figure 2.3 : Illustration of effect of dispersion on peak mean square incoherent 
field. Peak of g 2 {t), normalized with respect to max(w 2 (/)) as a function of 

grazing angle. 8 = 0.0163, / = 20kHz, c, = 1500m s, v = 1.13, t s = lOOps. 
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First order and zero order field strength 



field strength (dB) 

Figure 2.4a : First-order and zero-order field strength as a function of depth. 
Field strength is given in dB with respect to incident pressure. Solid line 
includes evanescent component in first-order calculation, dashed line does not 

include evanescent component in first-order calculation. Incident angle is 20°. 

Critical angle is 27.75°, / = 20kHz, a = 0.1m, y = 3, w- 2 = 2xl0 _5 m, p = 2, 
8 = 0.019, cj = 1500 m/s, v = 1.13, k\h = 0.254. 
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First order and zero order field strength 



Figure 2.4b : First-order and zero-order field strength as a function of depth. 
Field strength is given in dB with respect to incident pressure. Solid line 
includes evanescent component in first-order calculation, dashed line does not 

include evanescent component in first-order calculation. Incident angle is 20 . 

Critical angle is 27.75°, / = 20kHz, a = 0.5m, y=3, w 2 = 6.2 x 10" 5 m, p = 1.96, 

8 = 0.0163, ci = 1500 m/s, v = 1.126, /cj/i = 1.00. 
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First-order mean square incoherent field 



Figure 2.5a : Time-dependent mean square incoherent field at depth 0.2 m below 
the mean interface. Point source height above interface is 0.528 m ; transverse 
distance between source and receiver: 2m. Source frequency = 20kHz, 
— i 00 p 5 . Solid line uses Eq. 2.85; dotted line is obtained from setting 

g[t) - u(t). Same sediment parameters and incident angle as Fig. 2.4a. 
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First-order mean square incoherent field 
xIO ' 5 



Figure 2.5b : Time-dependent mean square incoherent field at depth 0.2 m below 
the mean interface. Point source height above interface is 5.5 m; transverse 
distance between source and receiver: 2m. Source frequency = 20kHz, 
t s = 100pj. Solid line uses Eq. 2.85; dotted line is obtained from setting 

g(t) = u(t). Same sediment parameters and incident angle as Figs. 2.4a, 2.5a. 
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First-order mean square incoherent field and 
square of zero-order coherent field magnitude 



Figure 2.6a : Square of zero-order coherent field magnitude shown with dotted 
line, first-order mean square incoherent field shown with solid line. All 
parameters identical to Fig. 2.5a. 
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First-order mean square incoherent field and 
square of zero-order coherent field magnitude 


xIO ' 5 



Figure 2.6b : Square of zero-order coherent field magnitude shown with dotted 
line, first-order mean square incoherent field shown with solid line. All 
parameters identical to Fig. 2.5b. 





CHAPTER 3 


First-order perturbation solution for rough fluid-fluid interface 


3.1 Introduction to chapter 3 

In the numerical computations in chapter 2, and in the simulation in chapter 
4, the water above the seafloor is modeled as a homogenous fluid, and the 
sediment is modeled as a lossy homogenous fluid. The rough interface between 
the water and sediment is modeled as a 2-D random process with a truncated 
power-law roughness spectrum. The first-order fluid-fluid perturbation results 
used in chapters 2 and 4 are derived in this chapter. The power-law filter used in 
chapters 2 and 4 is also discussed in this chapter. 

Perturbation derivations are also given in chapter 5. However, the 
derivations in chapter 5 allow the sediment bulk properties below the lowest point 
on the interface to be arbitrarily stratified, and the approach allows evaluation of 
higher orders. The goal of this chapter is to obtain the zero-order coherent field 
and the first-order incoherent field for the fluid-fluid rough interface problem with 
homogenous half spaces. In addition to derivation of an expression for the field 
scattered back into the water as in Kuo (1964) and Mourad and Jackson (1989), 
a solution for the field scattered into the sediment is also derived. Although the 
results in this chapter are derived for the two homogenous fluid case, the form of 
the first-order scattered field solution derived in this chapter is shown to be more 
general in chapter 5. 

3.2 Scattering problem description 

A 2-D slice of the 3-D scattering problem is shown in Fig. 2.1. An arbitrary 
pressure field, \|/,(r), is incident on a 2-D rough interface, z = /i/(R), separating 



the homogenous lossless fluid in medium 1 (water) from a homogenous lossy 
fluid in medium 2 (sediment), resulting in the pressure field, i|/y(r), scattered 

back into the water and the pressure field, y 2 -( r ). scattered into the sediment. 
The mean of /(R) is taken to be zero for convenience and without loss of 
generality. Substituting the plane wave representations (Devaney and Sherman, 
1973) for the field above the interface (Eqs. 2.5 and 2.6) and the field below the 
interface (Eq. 2.9) into the equation for the continuity of pressure (Eq. 2.3) yields 

JA 4 / ,-(K)e“ iilPl(K)/(R) +'P/(K)e i ^ l(K)/(R) e lKR 

= (3.1) 

Substituting the plane wave representations for incident and scattered fields into 
the equation for continuity of normal velocity (Eq. 2.4) yields 

—J d 2 Kk$i(K) -^(K)e _I * lPl(W(R) +^/(K)e zfclPl(K)ft/(R) e ?KR 
J d 2 KK • V(/z/(R)) ^(K) e " /fc, Pi (K) ^ (R) +T / / (K)e^ lPl(W(R) e* KR 

= -J- 

P2 

- — f d 2 KK ■ V(hf( R))^2_(K)e~ iKk ^ 2(K)hm) e iKR 
P2 J (3.2) 

where the unit vector normal to the interface in Eq. (2.4) is given by 


I-W(R) 

” |£ - V/j/(R)| ’ 


(3.3) 


and V represents the gradient over the transverse coordinates. 
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Since 8«1, the expressions in Eqs. 3.1 and 3.2 are valid (Landau and Lifshits, 
1987). 

The plane wave representation for the field scattered from the interface, 
given in Eq. (2.5), represents the scattered field above the highest point of the 
interface. Likewise, the plane wave representation for the field scattered through 
the interface, given in Eq. (2.9), represents the scattered field below the lowest 
point of the interface. Using Eqs. (2.5) and (2.9) to represent the field on the 
interface, z = hf( R), in Eqs. (3.1) and (3.2) invokes the Rayleigh hypothesis (e.g., 

Jackson, Winebrenner, and Ishimaru, 1988). 

The goal of this chapter is to find the scattered field in terms of the incident 
field. Since 

'F f (K) = J d 2 K' T\ i (K, K')'F i - (K'), (3.5) 

and 

^ 2 _(K) = J /rr 12 (K,K')^(K'), (3.6) 

the scattering problem is solved when the T-matrices T u and T u are found from 
Eqs (3.1) and (3.2). 

3.3 First-order and zero-order perturbation theory solution for scattered 
field 

The perturbation theory approach assumes that the product of the wave 
number in the water and the RMS surface height, k\h, is small. The fields and 
T-matrices are expanded in a power series in k\h, and the problem is solved one 
order at a time. In this chapter, we are only interested in the first-order solution 
used in chapters 2 and 4, and therefore keep only the zero-order terms and the 
first-order terms of the scattered fields, 
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'F(K) = T (0) (K) + (fc 1 fi)'F (1) (K), (3.7) 

and the zero-order and first-order T-matrix components, 

r(K,K/) = r (0) (K,K,-) + (M)r (1) (K,K f ). (3.8) 

Likewise, exponential terms In Eqs. (3.1) and (3.2) are also expanded in the first 
two terms of a power series, 

e i*!MW(R) £ i + ifcjP^KMftR). (3-9) 

The Fourier transform of /(R) is given by 

F(K) = —Ur [d 2 Rf(R)e~ iKR . ( 3 -10) 

(2rc) 2 

Therefore /(R) and V/(R) can be expressed as 

/(R) = Jj 2 rF(K')e iK ' R , (3- 11 ) 

V/(R) = J A'/K'F(K')e' K,R . (3-12) 

Using Eqs. (3.9) and (3.11) in Eq. (3.1) yields an approximate expression for the 
continuity of pressure boundary condition appropriate for first-order and zero- 
order fields: 

jA[^(K) + 'F / (K)]e /KR 

+J d 2 Ki{k x h% (K)[^(K)(-1) + ^/(K)]J A'F(K')e‘ (K+K ' )R 

= j A^ 2 _(K)e /KR - jrf 2 ^i(K/t 1 ^ 2 (K)^ 2 _(K)p 2 ^T(K0e' (K+K,)R . (3.13) 

Multiplying both sides of the above expression by e -lK R , and integrating this 
equation over R by means of the operator J d 2 R , yields 
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fd 2 K ['F, (K) + (K)]5(K - K") 

+J d 2 K i(kih)$i (K)[-'F I - (K) + 'Fy (K)] J d 2 K'F(K')S(K + K' - K") 

= |A v F 2 _(K)8(K-K") 

-j J 2 /sr<K^^)P 2 (K)^ 2 _(K)J d 2 K'F(. K')5(K + K' - K"), (3 14) 

where 8(K) is the 2-D Dirac delta function. Integrating the above expression 
over K and over K' results in the following simplified expression for the 
continuity of pressure boundary condition appropriate for first-order and zero- 
order fields: 


V F J (K") + 'Fy (K") 

+} rf 2 Ki(* 1 /.)3 1 (K)[-'P,(K) + - K) 

= T 2 .(K") -1 <i 2 i:i(Kfc 1 *)|5 2 (K)'P 2 .(K)F(K" - K) 


Following the above procedure, we now use Eqs. (3.9), (3.11), and (3.12) in Eq. 
(3.2) and obtain an approximate expression for the continuity of pressure 
boundary condition appropriate for first-order and zero-order fields: 


J Apj (K)[-^,(K) + 'Fy(K)]e lKR 

+/J d 2 K(k l h)^ (K^K) + 'Fy(K)]J d 2 K'F(K')z i(K+K ’ >R 

-ijd 2 K(k l h)[ x ¥ i ( K) + 'F/fK)]J d 2 K' ^^F(K / )e , ’ (K+K ' ) ' R 

= --J J 2 /i:K(3 2 (K)'F 2 _(K)e i ' KR 
P J 



J d 2 K(k l h)K 2 p^K)^.^)! J 2 ^T , F(K / )e ,(K+K) ' R 
J d 2 K(kihy¥ 2 -( K)} d 2 K' ^-F(K')e i(K+K,) ' R . 


(3.16) 


Multiplying both sides of the above expression by e ,K R and integrating this 
equation over \d 2 R yields 
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J d 2 K$ j (K)[-^(K) + H / / (K)]5(K" - K) 

+/J d 2 K (K)[^(K) + ¥/( K)]J d 2 K'F{ K')5(K + K' - K") 

-if d 2 K fakpi (K) + Vf (K)]J d 2 K' ^F(K')6(K + K' -K") 

= d 2 K kP 2 (K)^- (K)5(K" - K) 

+i- \d 2 K{k x hy 2 ^ 2 2 {K)^2-{K)\ d 2 K'F( K')8(K + K' - K") 

- - J d 2 K (M)^- (K) J d 2 K' F(K')5(K + K' - K"). 


Integrating the above expression over K and K' results in the following 
simplified expression for the continuity of normal velocity boundary condition, 
appropriate for first-order and zero-order fields: 


PjCK'^-^K'O+^K")] 

+ij d 2 K(k l h)^ (K)[^(K) + 'F / (K)]F(K" - K) 

-ijd^fahpii K) + ^/(K)] K ( ^2~ ~F(K"-K) 

ACJ 

= _I K p 2 (K")^ 2 -( K ") 

9 

f d 2 K (fe 1 /?)K 2 P2 (K)^ / 2-(K)F(K" - K) 

P J 

- l -\d 2 K{k l hyY 1 -( K) K {K jl ~ K ~W-K). 


3.3.1 Zero-order field (flat surface reflected and transmitted field) 

Substituting ^(K) for ^-(K) and ^} 0) (K) for v Fy(K) into Eqs (3.15) and 

(3.18), and keeping only zero-order terms, yields the zero-order equations for 
continuity of pressure and continuity of normal velocity: 
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4 / I (K") + x pj 0) (K") = 'F^ ) (K") (3.19) 

%(K") - ^(K") = K q 2 ; ^ 2 -(K") ■ (3.20) 

J PPl(^ ) 

Solving Eqs. (3.19) and (3.20) for and 'f} 0) yields 


'F} 0) (K) = r(K)'F J -(K) 

and 

v F^ ) (K) = T(K) v f / ; -(K), 
where 

= 1-kP 2 (K)/pPi(K) 

V ; l + Kp 2 (K)/pp 1 (K) 


(3-21) 


(3.22) 


(3.23) 


is the flat surface reflection coefficient, and 


T(K) = 1 + T(K) = 


_ 2 _ 

l + Kp 2 (K)/pp 1 (K) 


(3.24) 


is the flat surface transmission coefficient. The zero-order field is therefore the 
field transmitted and reflected from the flat zero-mean surface. Comparing Eqs. 
(3.21) and (3.22) with Eqs. (3.5) and (3.6) shows that the reflection coefficient 

can be expressed in terms of the zero-order T-matrix, , 


r(K) = r 1 ( 1 0) (K',K)8(K'-K), (3.25) 

and the transmission coefficient can be expressed in terms of the zero-order 
T-matrix, r/ 2 0) , 

T(K) = r 1 ( 2 0) (K / ,K)6(K' - K). 


(3.26) 
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3.3.2 First-order scattered field 

Substituting Eq. (3.7) into Eqs (3.15) and (3.18), and keeping only first-order 
terms (terms with a k\h factor), yields the first-order equations for continuity of 

pressure, 

^(K") + J A/p^K^-^K) + T} 0 ) (K)]F(K" - K) 

= ^(K") - J A/Kp 2 (K)^°J(K)F(K" - K), (3.27) 

and continuity of normal velocity, 

pPi(K ) 

Kp 2 (K") f 

+-^-f Apf(K)k(K) + 'Fi 0 ) (K)lF(K' / -K) 

kP 2 (K") J L f (3.28) 

-£-f d 2 K k- (K) + ^ 0) (K)j K (K " F(K" - K) 

kP 2 (K") J ' f J *? 

= - 4 / | 1 J(K") 

+ ;- J A K 2 pi(K)^(K)^ - K) 

- i - f A^ 0) (K) K • (K " ~ F(K" - K). 

kP 2 (K") J 2 “' k} 

3.3.2.1 First-order field scattered from interface back into water 

Equations (3.27) and (3.28) are added together to eliminate Eqs. 
(3.21), (3.22), and (3.24) are used to express 4jS 0) and in terms of 
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r(l) 


f 


vpV'^K") 


1 + 




PMEl 

Kp 2 (K") 


\ 


J 


ip 


Kp 2 (K") 
zp 


Kp 2 (K") 


- f d 2 K$ (K)^ (K)(l + r(K))F(K" - K) 

r 

JA'PKKXi+HK)) 


K • (K" - K) 


F( K" - K) 


kP 2 (K") 


l —-j AK 2 pi(K)^(K)(l + r(K))F(K" - K) 


JA^(K)(l + r(K)) 


K • (K" - K) 


F(K" - K) 


Kp 2 (K") J ‘' ' " k\ 

+J A ip! (KWi (K)(l - r(K))F(K" - K) 

-J d^iKfciWi (K)(l + T(K))F(K" - K). 

Noting Eq. (3.5), and multiplying the above equation by Kp 2 (K")/p, we 
equivalent expression in terms of the first-order F-matrix: 


Pi(K") 


r 




kMIO 

PPi(K '0 


\ 


J 


7j§)(K", K) = \ i 


i(l-T(K)) K ^ {K \ (K) 

P 


-i(l + r(K)) 
+i(i+r(K)) 




1 „ 2 d 2 


Kp 2 (K") 


\ 


Pf(K )-—k P 2 (K) + 

P P 


kP 2 (K) 


) 


K • (K" - K) 


: l 


r 


i-i 


\ 


V 




J 


F(K"-K). 


For convenience, define H\\ such that 
r/J^K^K) = H U (K",K)F(K" - K). 
Substituting 


K-(K"-K) K K" 


1 + Pl (K), 


p ?( K )-( K p 2 ( K )) 2 =£* 2 

*1 


k 2 


K‘ 


= 1 — K' 


(3.29) 

obtain an 

(3.30) 

(3.31) 

(3.32) 

(3.33) 
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into Eq. (3.30), and using Eq. (3.31), yields 
P, (K")f 1 + , (K",K) = i(l - r(K)) 


v l p p P v ' 


*V 1 \ \r \r" 

+l(l + r(K)) 1-- —y- 

v U pj *? 


(3.34) 


With the use of Eq. (3.24), the left side of Eq. (3.34) is expressed in terms of one 
plus the reflection coefficient: 

Hl , (K",K)=-(i+r(K"))(i - r(K))^* 

(n-r(K”))(n-r(K)) r Kp,(K>p 2 (K) 11 1 ( w 

iPl(K") t p P P 1 'J 


(i + r(K"))(i + r(K)) ((i k • k" 

2*P 1 (K") Up J k? 


(3.35) 


The first factor on the right side of the above equation and the first two factors in 
the second term are converted from one plus the reflection coefficient to one 
minus the reflection coefficient through the following equality, obtained from Eq. 


(3.23): 

1 + T(K) pfo(K) 
l-r(K) Kp 2 (K)' 

The expression for H\ j follows from Eq. (3.35): 

H ' 1 (K "’ K)=+r ( K ">X! + HR)) 

+^K",K)(1 - r(K"))(i - r(K))}, 


(3.36) 


(3.37) 
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where 


a(K",K) = 1- 



+ 



\ 

-1 

) 


K" K 



(3.38) 


and 

b(K",K) = P!(K")Pi(K)(p -1). (3.39) 

The first-order incoherent field is given by k\h multiplied by the 2-D inverse 
Fourier transform of Eq. (3.5), where 7jy is given by Eq. (3.31) and H n is given 
by Eq. (3.37). 


3.3.2.2 First-order field scattered into sediment 

The field scattered into the sediment follows from rearranging the terms in 
Eq. (3.27), and expressing the reflected and transmitted fields in terms of the 
incident field by using Eqs. (3.21) and (3.22), 

^(K") = ^(K") - J d 2 K ip^KXl - T(K))F(K" - K) 

+J Ahcp 2 (K)(l + T(K))F(K" - K). ( 3 .40) 

As in Eq. (3.31), define in terms of the first-order F-matrix, T $, and 2-D 
Fourier transform of the interface profile function F, 

T${ K",K) = tf 12 (K",K)F(K"- K), (3.41) 

where 

H$( K",K) = H§{K",K) - fP^KXl - T(K)) + ncp 2 (K)(l + T(K)) (3.42) 

follows from Eqs. (3.40) and (3.41). Substituting Eq. (3.37) for H\\ into the 
above equation yields 
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(K " K)= 2iprj ^ (K " K)(l+r(K " ))(1+r(K)) 

+b(K", K)(l - r(K"))(l - r(K))} 

-ip, (K)(l - r(K)) + iKp 2 (K)(l + r(K)). (3.43) 

The following identity is shown by using Eq. (3.36): 

2/Pi (K")[-i‘Pi (K)(l - T(K)) + /k(3 2 (K)(1 + T(K))] 

+b( k",k)(i - r(K"))(i - r(K)) 

= -Pj (K)(3j (K")(p -1)(1 - r(K))(l + r(K")). ( 3 . 44 ) 

Substituting Eq. (3.44) into Eq. (3.43) yields 

H, 2 (K", K) = ^T[a(K",K)(l + r(K)) - i*K",K)(l - r(K))]. (3.45) 

where 7 is the flat surface transmission coefficient given in Eq. (3.24). The first- 
order approximation to the T-matrix for forward scattering into the sediment is 
obtained from Eqs. (3.6), (3.41), and (3.45). 

3.4 Time-dependent mean square incoherent field approximation 

Since the first-order T-matrices are proportional to the 2-D Fourier transform 
of the zero-mean interface, their expected value over the rough interface 
ensemble is zero. In addition, the zero-order component of each T-matrix is only 
dependent on the zeroth moment of the rough interface ensemble. Therefore 

ThfK.K,) = (*,/. jj<j)(K,K / ).(* 1 A)H 12 (K.K i )F(K-K i ) (3.46) 

and 

f, ,(K,K,) = (*,/ ! )jfl>(K,K i ) = (*,/.)», ,(K,K,)F(K - K,) (3.47) 

are the first-order approximations for the incoherent T-matrices given in Eq. 
(2.18). Using the first-order approximation to the T-matrix in Eq. (3.41), 
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(t I *) 2 (7<|)(K',K; > <D 1 )jg ), (K',Kr,<D 2 )) 

= (* 1 ft) 2 W 12 (K',K;,(fl 1 )H 1 * 2 (K",Kr,C 0 2 )(f|(K'-K;)F 1 (K"-K;')) (3.48) 

is the first-order approximation to the incoherent F-matrix correlation function, 
given in Eq. (2.33), for forward scattering through a fluid-fluid interface. The 
dependence on angular frequency is shown explicitly in the F-matrix and H\ 2 - 

Using the identity 

/j 2 (f(K" - K t )F(K' - K,)) = W(K' - K,)8(K" - K') (3.49) 

together with Eq. (2.33) and Eq. (3.48) yields 
^ 12 (k',k",k;,k;',o) c + co',o) c +©") 


= ^ 2 H 12 (K',k;,co c -Ko')H; 2 (K",Kr,co c + <o")w(K - (K< + Kfj/2) , (3.50) 


where Wis the roughness spectrum. For the case of a wide sense stationary 
process (WSS), 

W(K) = —Ur f d 2 R C(R)e -,K R , (3.51) 

(2tc) z j 


where 

C(R tf ) = (/(R + R J )/(R)) (3.52) 

is the surface autocorrelation function. For this special case, the identity in Eq. 
(3.49) is simply derived by making the substitutions 

Rj = Ri - R 2 , R a = (Ri + R 2)/ 2 (3.53) 
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and using Eqs. (3.51) and (3.52). 


3.4.1 Seafloor roughness spectrum 

The two dimensional Gaussian random process /(R) describing the seafloor 
surface is assumed to be an isotropic power-law (Fox and Hayes, 1985; Goff, 
and Jordan, 1989; Briggs, 1989) with a roughness spectrum 


W( K) = 



(3.55) 


In the notation of Mourad and Jackson (1989, 1993), the spectral strength and 
spectral exponent are represented by the parameters w 2 and y, respectively. For 
2 < y < 4, this surface random process is an example of a fractal (Mandelbrot, 
1983), and the perturbation theory requirement that k x h be small is not satisfied. 
In fact, the RMS height is infinite, and the random process is not WSS. However, 
since this random process is stationary in increments, Eq. (3.49) is valid 
(Ishimaru, 1978b), and the large scale roughness can be removed by a filtering 
operation (Jackson, Winebrenner, Ishimaru, 1986). 

The power-law spectrum is an estimate of the roughness spectrum over all 
spatial frequencies. The seafloor roughness spectrum on scales over one meter 
is poorly known, and therefore assuming a roughness spectrum that is a power- 
law over all spatial frequencies (Eq. 3.55) may not be valid. Large scale 
components of the roughness spectrum are set by the spatial scale of the 
experiment. In the case of a tank experiment, moderate to large scale roughness 
will not be present. Removal of the large scale roughness in Eq. (3.55) is 
therefore appropriate. 


Let /(R) be the surface with the power-law roughness spectrum given in Eq. 
(3.55), with 2-D Fourier transform, F(K). A filtered surface, /(R), is obtained by 
subtracting a moving average (Papoulis, 1984) of /(R) from itself: 
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27t a 


-4 vr 

/(R) = /(R)-— j JdfJ dR'f(R - R') . 


(3.56) 


Equation (3.56) can be expressed as a 2-D convolution of /(R) with a filter .s(R), 


/(R) = -l- ? |/(R-R'>(R'), 

{2%y J 


(3.57) 


./(R) = /(R)®j(R), 

where s(R) is of the form 


(3.58) 


s(R) = (2 tc) 2 5(R) - m(R), 


(3.59) 


and 


Mi( r) = ( 2 ^5W 

Tea 

is a moving average filter. The circ function is defined as (Goodman, 1968) 


(3.60) 


circ(/?) = 


1 R< 1 
0 otherwise 


(3.61) 


Filtering the original power-law random process with roughness spectrum given 
in Eq. (3.57) with the filter s(R) corresponds to subtracting from each point on the 

random surface, /(R), an average of the surface over a radius a about each 
point. The filter's Fourier transform follows from Eq. (3.59), 

S(K) = 1 - M(K) (3.62) 


where 


M(K) = 


2J\{Ka) 

~Ka 


(3.63) 
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The resulting roughness spectrum obtained by filtering the surface with the 
power-law roughness spectrum given in Eq. (3.55) using the filter in Eqs. (3.62) 
and (3.63) results in the 2-D surface random process given by 


W(K) = 


Wo 

F 


r 


i- 


V 


2 JjiKay f 
Ka j 


(3.64) 


Instead of subtracting a uniformly averaged area around each surface point 
from the surface, subtracting a Gaussian weighted average around each surface 
point from the surface results in a roughness spectrum without the Bessel 
function oscillations. In this case 


m 2 (R) = (2 ji) 


2 e 


-R 2 12a 2 


a-j2n 


(3.65) 


and the filtered, or truncated, roughness spectrum is given by 


W( K) = 



with RMS height (see Appendix C) 


(3.66) 


oo 



1 


2k j W(K)KdK = 


0 



(3.67) 


Here, T refers to the gamma function — not the reflection coefficient, and not the 
two-frequency mutual coherence function. 

The roughness spectrum given in Eq. (3.66) is the roughness spectrum used 
in this dissertation. The choice of the filter parameter, a, is discussed for each 
numerical calculation. The larger this parameter is, the closer the filtered 
roughness spectrum is to the original power-law (Eq. 3.55). Subtracting a 
moving average of the power-law random process from itself results in removal 
of low frequency components of the original power-law roughness spectrum. 
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Plots of this roughness spectrum as a function of K for two different values of the 
filter parameter, a, are given in Fig. (3.1). 

3.4.2 Incoherent r-matrix correlation function approximation in Eq. (2.37) 


In order to evaluate the approximation in Eq. (2.37) using first-order 
perturbation theory, K-K^/2 is substituted for K', and K + K^/2 is substituted 

for K" in Eq. (3.50) 

^ 12 (K-K^/2,K + Kj/2,K/,K",co c +co',o) c + 0)") 

= k?H u {K-K d /2,Ki,u> c +a')H* 2 (K + K d /2,K?,a> c +co") 


xW(K-(K' + Kf)/2), 

where K d is given in Eq. (2.34) and is repeated here for convenience. 


K // / 

j . IV j 


(£> c + 0 ) 


G)"-©' 


COS 0 yd/ - 


© c + © 


■cos 9, a / 


cos 0,6/, 


(3.68) 


(3.69) 


where 0, is the incident grazing angle and a,- is the unit transverse direction 
vector of the incident field. From Eq. (2.37), 0 n > n Eq. (3.68) is approximated by 

2 (K - Krf /2, K + Krf /2, k; , Kf, © c + ©', © c + ©") 

= <? 12 (K,K,-,© c ) 


s* 1 2 |ff 12 (K,K lS ©J 2 W(K-K f ), 


(3.70) 


where 


K/ — (© c /q)cos0,a,- 


(3.71) 
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The roughness spectrum, VF(K-(K; + K/)/2), is not a function of 
For a narrowband signal, it is a slowly varying function of 


(3.72) 


(o “7o\ 

CO a =£D c +-T-’ 

except for K= k\. However, scattered energy with K>k 2 is evanescent. From 
inspection of Eqs. (3.38), (3.39), and (3.45), the quantity 

HuiK-Kd/lXi^c + to'jHniK+Kd/l’K^c + to") depends on only in the 
last term of Eq. (3.38). In addition, the dependency on (o d is primarily in the 
factor 7/(3 outside the square brackets in Eq. (3.45). 


7(K + K d /2) __2p- (374) 

(3i(K + K d /2) k(3 2 (K + K rf /2) + pft(K + K d /2) 

Although the magnitude of Eq. (3.68) varies slowly with frequency, the phase 
becomes very dependent on a> d when either |K + K^/2| or |K-K^/2| is greater 

than both k k x and k x . For a narrowband signal, this situation occurs for 
evanescent scattered waves and possibly scattered waves at shallow grazing 
angles. However, as shown in chapter 2, including the evanescent scattered 
field is only necessary within a wavelength or two of the interface, depending on 
the roughness. 

Figure 3.2a is a plot of the magnitude of Eq. (3.68) (excluding the roughness 
spectrum, W) and Fig. 3.2b is a plot of its phase as a function of K = |K| and 

fd =®j/( 2 ^)- ( 3 - 75 ) 


The roughness spectrum is not included in these plots, because it is not a 
function of f d . The direction of K is the same as the incident transverse wave 

vector in this example. The center frequency is 20 kHz, and the sediment 
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parameters are the same as Fig 2.4b. Figure 3.2a shows that the magnitude 
becomes slightly dependent on f d only for shallow grazing angles. Figure 3.2b 
shows that the phase of Eq. (3.68) becomes dependent on large f d only for the 

evanescent scattered field and for scattering at shallow grazing angles — even 
for a relatively wideband signal. The simplifying approximation for e in Eq. 
(2.37) is valid for all values of K, except those values corresponding to 
evanescent waves and those corresponding to scattering at small grazing 
angles. Because of loss in the sediment, energy scattered at these small grazing 
angles is insignificant compared to energy scattered at larger grazing angles. 
Therefore, Eq. (3.68) can be considered to be independent of f d . 


3.4.3 First-order time-dependent mean square incoherent field 


Substituting £ 12 (K,K t ,co c ) in Eq. (3.70) into Eq. (2.50) yields the first-order 

expression for the time-dependent mean square incoherent field scattered 
through a rough fluid-fluid interface due to an incident plane wave: 


I 2 (r,t) s fd 2 Kk}\H l2 {lLK i ,(O c fw(K-K i )e 2kl M K MK*®c))z g 2 ^_^ (3-76) 


The differential cross section from first-order perturbation theory is found by 
substituting Eq. (3.70) into Eq. (2.61). 


o 12 (a,d / )s^^|ff 12 (l^K /> (D c f W(K-K f ). (3.77) 

Y 

Although the differential cross section in this approximation is actually second 
order in k x h (Thorsos and Jackson, 1989), we refer to it as the first-order 
differential cross section, and will refer to the corresponding mean square 
incoherent field in Eq. (3.76) (/ 2 (r,r)) as the first-order mean square incoherent 
field. Another first-order expression for / 2 (r,r) for an incident plane wave that 
neglects evanescent waves is obtained using Eq. (3.77) in Eq. (2.63). The effect 
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of neglecting evanescent waves is discussed in chapter 2. The first-order 
differential cross section is also substituted into Eq. (2.85) to obtain the first-order 
expression for / 2 (r,r) due to an incident field from a point source. This resulting 
expression is used in the simulations of experiments that are discussed in 
chapter 4. 

3.5 Summary 

First-order perturbation solutions of the scattered field from and through a 
rough interface between two homogenous fluids are derived in this chapter. 
These results are used to show the validity of approximating the incoherent 
r-matrix correlation function as only dependent on the center frequency of a 
narrowband source signal (Eq. 2.37). The results derived in this chapter are 
primarily used in chapters 2 and 4. The Gaussian weighted filtered power law 
(Eq. 3.66) is used for all numerical computations in this dissertation. The filter 
parameter, a , is chosen to create a spectrum appropriate for each particular 
example. For example, in the simulation of scattering into the sediment in a tank, 
discussed in chapter 4, a small value is used for this filter parameter. In ocean 
experiments, the value of this parameter is considerably larger. 
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Magnitude 



K (1/m) 


80 


Figure 3.2a : Magnitude of Eq. (3.68), incoherent T-matrix correlation function (divided 
by the roughness spectrum), incident angle =20 , f c = 20 kHz, p = 1.9608, 5 = 0.0163, 
v= 1.12565, cj =1536 m/s. 
















































































CHAPTER 4 


The effect of roughness on acoustic penetration of the seafloor as 
given by a fluid-fluid perturbation model and comparison with recent 

sediment penetration experiments 

4.1 Introduction to chapter 4 

Simulation results of a sediment penetration tank experiment by Boyle and 
Chotiros (1992) and an ocean penetration experiment by Chotiros (1995) are 
presented in this chapter. For the incoherent intensity calculations, the first-order 
perturbation theory results derived in chapter 3 are used in Eq. (2.85); the 
incoherent intensity in the sediment is obtained by dividing Eq. (2.85) by the 
characteristic impedance of the sediment (p 2 c 2 ). In the ocean experiment 
simulation, the time domain zero-order coherent intensity is calculated using Eq. 
(2.92) (divided by p 2 c 2 ). Note that the incoherent intensity calculated using first- 
order perturbation theory is proportional to the second power of k\h, and the 
coherent intensity is computed to an accuracy of only zero order. The 
contribution of the second-order coherent intensity contribution is derived and 
demonstrated in chapter 5. The effect of sound speed gradients in the vertical 
direction on the field scattered back into the water is considered in chapters 5 
and 6. Although the influence of sound speed gradients on the field scattered 
through the interface is beyond the scope of this dissertation, possible effects are 
mentioned in chapter 6. 

4.2 Background 

Figure 4.1 depicts the scattering model derived in chapter 2. The source 
pulse originates at a point in the water (medium 1) and propagates to the «th 
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surface patch, its mean square incoherent field magnitude is multiplied by the 
differential cross section of the patch, and then it propagates to the receiver in 
the sediment (medium 2). In addition to spherical spreading from the source to 
the patch, and from the patch to the receiver, there is a sediment loss term. The 
contributions of all patches are summed up and divided by the characteristic 
impedance of the sediment, yielding the total incoherent intensity below the 
rough interface. Figure 4.1 includes the term g(t) derived in chapter 2. We 
argue later that the baseband source signal, u(t), can be substituted for g(r) for 
both simulations. Substituting u{t) for g(t ) in Eq. (2.85) yields:*-- 


/jlrjJsAAVr,; 


2 g»K/; g "i) |, _ ^ f e -2 bn(*2 )r 2 „ 


n 


2 2 
r \n r 2n 


(4.1) 


where tj n is the sum of the propagation time from the source to the «th patch 
and the propagation time from the «th patch to the receiver. Note that the 
incoherent intensity is given by 7 2 (r,r)/(P2 c 2)> where p 2 c 2 is the characteristic 
impedance of medium 2. The first-order perturbation expression for scattering 
into the sediment, a„, is derived in chapter 3 and given in Eq. (3.77). The exact 
expression for the time domain coherent intensity is given by Eq. (2.92) divided 
by the characteristic impedance of medium 2 (also see Fig. 4.1). Since we use 
the zero-order field to represent the coherent field, the transmission coefficient, 7 
is substituted for j? 12 in this figure — which then becomes Eq. (2.92). Following 
the definition of the coherent reflection coefficient in Thorsos (1990), ^12 is 
defined in chapter 5 to be the coherent transmission coefficient. The power 
spectrum used for both simulations is that of a filtered power law (Eq. 3.66). 
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4.3 Tank experiment simulation 

This simulation determines the incoherent field scattered into the sand, 
assuming a small amount of roughness at the interface. A diagram of the tank 
experiment is shown in Fig 4.2. 

4.3.1 Simulation parameters 

A source in water just under one half meter above sand transmits a signal 
with Gaussian envelope and 60 kHz carrier frequency. The baseband source 
signal used in the simulation is 

2/2 

u(t) = constant x e _f ' ts , (4-2) 

with t s = 333 p.s. This signal approximates the time resolution achieved in the 
tank experiment. The water sound speed, ci = 1489m/s, and signal center 
frequency of 60 kHz correspond to a wavelength of X = 2.5cm in the water. In 
order to agree with the experiment described in Boyle and Chotiros (1992), the 
value for constant is chosen such that the field lm from the source is 191 pPa. 
The sound speed in sand is given by Boyle and Chotiros (1992) to be 
c 2 = 1 675 m/s . The loss parameter in the sand is not given, and is assumed to be 
a value appropriate for sand, 8 = 0.015, or 0.49dB/m/kHz . No roughness 
parameters for the experiment characterizing the rough interface for the sand are 
given. We assume only a small amount of roughness, and set the filter 
parameter in Eq. (3.66) to a = 0.01m. The spectral exponent and spectral 

strength in Eq. (3.66) are chosen to be y = 3 and w 2 = 2.17xl0 -5 , respectively, 
resulting in an RMS height (Eq. 3.67) of h = 1mm and k\h= 0.25. This roughness 
spectrum is well within the range of validity for first-order perturbation theory for 
scattering through a rough interface (Thorsos, 1996a). Of the three roughness 
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parameters (a,y,w 2 ), y = 3 was chosen arbitrarily, but it is consistent with 
observed seafloor roughness spectra. Values for w 2 and a in Eq. (3.66) are 
chosen to provide only a small amount of roughness while yielding scattering 
intensities close to those obtained by Boyle and Chotiros. The roughness 
spectrum is shown in Fig. 4.3. 

4.3.2 Simulation results 

As mentioned, the baseband source signal u(t ) is substituted for g(t) in Eq. 

(2.85) for the simulations (Eq. 4.1). The justification for using this simpler model 
in the tank experiment simulation is shown in Figs 4.4a and 4.4b, where the time- 
dependent mean square incoherent fields using u(t ) and g(t) are compared for 
two receiver depths, with 2 m transverse distance between source and receiver. 
At the lower depth, the two cases are almost indistinguishable; at the depth of 
0.08cm , the differences become more noticeable, but are still negligible. 

The incoherent intensity (first-order) as a function of time and horizontal 
distance between source and receiver is shown in Figs 4.5a, 4.5b, and 4.5c for 
receiver depths of 0.05 m, 0,33 m, and 0.614 m. Since the baseband source 
signal is Gaussian, the arrival time is the peak of the intensity pulse. The source 
level is given by Boyle and Chotiros (1992) to be at 191 dB referenced tolpPa at 
1 m (in the sediment). Figures 4.4a and 4.4b are vertical slices of Figs 4.5a and 

4.5b, respectively (in Pa 2 instead of dB). The simulation results have been 
presented in Figs. 4.5a, 4.5b, and 4.5c for the purpose of direct comparison with 
Figs. 5, 6, and 7 of Boyle and Chotiros (1992). By choosing the appropriate 
roughness parameters, these simulations match measured intensity space - time 
behavior as presented by Boyle and Chotiros (1992), not only in terms of pulse 
arrival times, but also in terms of the resulting field strength. Although the 
simulation only contains the incoherent intensity, the coherent contribution, 
dominates only for incident grazing angles above critical. The curve plotted in 
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Fig 4.5 is the propagation time for the pulse, assuming the existence of a 
1200m/s refracted "slow wave." The arrival time of the incoherent intensity pulse 
coincides with this curve. This agreement occurs for all three depths, as well as 
over a wide range of incident angles. This model does not allow for slow wave 
propagation, and the "apparent slow wave" propagation can be accounted for by 
considering the diffraction path for the signal. A result due to Thorsos (1995) 
showing the relationship between an apparent slow wave and downward 
diffracted energy is presented in the next section. 

4.4 How roughness can mimic a slow refracted wave 

The purpose of this section is to determine the apparent slow wave speed for 
energy diffracted in a given direction. This derivation follows Thorsos (1995), 
and allows for the diffracted energy in the sediment to be in directions ranging 
from straight down to the compressional refracted direction (determined by 
Snell's law). Figure 4.6 shows a plane wave pulse propagating with a sound 
speed q and at a sub-critical grazing angle incident on a water sediment 
interface. Two rays strike the interface; the second intersects the interface at a 
distance d from the first, and at a time x later than the first. Both rays are 
scattered by the rough interface, and propagate into the sediment at a grazing 
angle 0 2 and at a sediment sound speed c 2 . The path AB, at a grazing angle 

e^ = e 2 -e', (4.3) 

represents the apparent slow wave refracted direction corresponding to the slow 
wave speed c 2 . All points perpendicular to this path are therefore taken to have 

the same total propagation time. From Fig. 4.6, 


and 


^ CiT 

COS0]- — 

a 


(4.4) 
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cos 02 = (4.5) 

d 

Solving Eqs. (4.4) and (4.5) for x, and equating, results in Snell's law for the 
apparent refracted slow wave direction: 

COS0! _ COS 02 / 4 g) 

ci c' 2 

Using Eq. (4.3) and 

c 2 = c 2 cos0' (4.7) 

in Eq. (4.6), and solving for 0', yields 


0' = tan" 1 



COS01 -COS02 



(4.8) 


The apparent refracted grazing angle, 02, follows from the above result and Eq. 
(4.3). The apparent slow wave refracted speed is given by Eq. (4.7). 

For the same sound speed ratio used in the simulation, Fig. 4.7 shows the 
dependence of the apparent slow wave velocity on the scattered grazing angle, 
0 2 . Figure 4.7 indicates that the apparent slow wave propagation speed of 
1200m/s in the tank experiment simulations (see Fig. 4.5) corresponds to an 

"average" scattered grazing angle of about 80°. Because of loss in the sediment, 
the average diffraction path is at large grazing angles in a downward direction, as 
other directions suffer larger absorption loss. 


Taking the scattered grazing angle to be 90° (straight down), we plot the 
effect of sediment sound speed and incident grazing angle on the apparent slow 
wave speed and apparent refracted angle in Fig. 4.8a and Fig. 4.8b, respectively. 

For this case, 0 2 = 90°, the apparent refracted speed is equivalently expressed 
as 
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V 1 + ( c 2 M ) 2 cos2 0 1 

It is evident from this example that the apparent refracted slow wave speed and 
the apparent refracted grazing angle do not depend strongly on the sediment 
sound speed, even for these extreme sediment compressional sound speed 
examples. Although the apparent slow wave speed increases with increasing 
incident grazing angle, it is relatively constant for incidence below critical grazing. 
Below critical grazing is the region of interest, since the refracted component 
dominates above this angle. 

Consequently, roughness can mimic a slow refracted wave (=1200m/s), 
even though the acoustic sound speed is considerably larger (1675m/s in this 
example). This tank experiment simulation shows that the effects of roughness 
must be considered in interpretation of penetration experiments. 

4.5 Ocean experiment simulation 

The tank experiment determines the arrival time of the pulse, but not the 
direction that the wave propagates. In considering the propagation of the 
penetrating acoustic field, an ocean experiment using a buried array of receivers 
(Panama City experiment, Chotiros, 1991, 1995) measures the direction and 
speed of the intensity front. The results of computer simulations of this 
experiment using Eq. (4.1) for the incoherent intensity and Eq. (2.92) for the 
coherent intensity are presented in this section. 

4.5.1 Experiment simulation layout 

The receiver locations used in the ocean simulation are taken from Fig. 3.3.1 
of Chotiros (1991). This array layout is also shown in Fig. 3.4 of Chotiros (1995). 
As in their analysis of the experimental data, the receiver shown on the y axis 
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and the receiver farthest out on the x axis (see Fig 3.3.1 of Chotiros, 1991) are 
not used in the simulation. The layout used in the simulation is shown in Fig. 4.9. 
The source is at a height of 4m above the sediment, positioned along the x axis 
at several different distances from the array corresponding to incident angles with 
respect to the array center of 0,- = 79°, 40°, 24.9°, 12.7°, and 9.7°. 

As in Chotiros (1995), the center frequency for the ocean experiment 
simulations is 20 kHz. The water sound speed, sediment sound speed, and 
sediment loss parameter are given by Chotiros (1995) to be q = 1536m/s, 
c 2 = 1729m/s, and 5=0.0163 (,514dB/m/kHz), respectively, resulting in a critical 

incident grazing angle of 0, =27.3°. The spectral strength and spectral exponent 
are consistent with Stanic et al. (1988), which are from measurements in the 
vicinity of the experiment ( w 2 = 6.2xl0~ 5 m and y = 3). The filter truncation 
parameter (Eq. 3.66), a = 05m, is chosen to be consistent with the scale of the 
experiment, resulting in k\h = 0.98 from Eq. (3.67). Note that the transverse 
distance between all receivers is within 2a = lm. Thorsos (1996a) shows that 
perturbation theory is valid for scattering through the water-sediment interface 
with k\h as high as 2.0. 

4.5.2 First-order intensity model 

As in the tank simulation, the baseband source signal u(t) is also substituted 
for g(r) in Eq. (2.85) for the ocean simulations. The baseband source signal is 
taken to be 

u{t) = , (4.10) 

with t s = lOOpsec. This pulse width is comparable to that achieved in the actual 
experiment. As mentioned later, the signal at each receiver is normalized such 
that the signal peak is the inverse of the number of array elements. Therefore, 
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the magnitude of the source signal is irrelevant to the simulation. The 
justification for using this simpler model (Eq. 4.1) for the incoherent intensity is 
shown in Figs 4.10a and 4.10b, where the mean square incoherent field on the 
most shallow receiver (0.05 m depth) is plotted using both the expression in Eq. 
(2.85) and Eq. (4.1), that is, u(t) is substituted for g(t ) in Eq. (2.85). Although the 
mean square incoherent field magnitude using the simplification g{t) = u(t) is 22% 
greater than the result obtained from Eq. (2.85) at a depth of 5cm, and 17% 
greater than the result obtained from Eq. (2.85) at a depth of 18cm, the arrival 
times are not affected by the approximation, g(t) = u(t). Furthermore, the error in 
incoherent intensity resulting from this simplification does not affect the results 
presented in this chapter. 

4.5.3 Coherent intensity model 

Analytical approximations for the integral in Eq. (2.91) can be found (for 
example, Brekhovskikh, 1980; Ishimaru, 1991). When the incident field is above 
critical, steepest descent solutions agree closely with the exact numerical 
solution. Furthermore, for narrowband signals, the received zero-order signal 
amplitude is determined by the magnitude of the Green’s function evaluated at 
the pulse center frequency. Therefore, the time domain zero-order received 
pulse can be determined by evaluating the steepest descent solution at the pulse 
center frequency. Although the approximate steepest descent analytical 
solutions agree with the numerical solution of Eq. (2.91) very closely for 
incidence above the critical angle, the solution near critical can differ significantly 
(Westwood, 1989). In addition, the Green’s function is very frequency dependent 
when the incident field is at or below the critical angle. For this reason, the zero- 
order transfer function or Green’s function is found by evaluating Eq. (2.91) 
numerically for a discrete range of angular frequencies, to. In order to obtain the 
zero-order output pulse at a receiver, the discrete Fourier transform of the source 
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signal in Eq. (3.21) is multiplied by the corresponding Green’s function, and the 
result is transformed back to the time domain. The resulting exact zero-order 
time-dependent field magnitude is squared and added to the time-dependent 
first-order mean square incoherent field. The simpler steepest descent solutions 
have been used for comparison purposes. 

Although the zero-order field will be seen to be somewhat greater than the 
actual coherent signal (Rice, 1951; Thorsos, 1990) for the values of k x h used in 
these simulations, Chotiros' intensity velocity algorithm is affected only if 
coherent and incoherent intensities are close in magnitude. As shown in chapter 
5, using only the zero-order field to represent the coherent intensity results in 
coherent intensity calculations about 10% greater than calculations using both the 
zero-order and second-order contribution. Including the second-order component 
of the coherent intensity would not produce noticeably different simulation 
results. 

4.5.4 Algorithm for estimating intensity front velocity 

Since the receiver locations in Chotiros' ocean penetration experiment are 
not known precisely enough to use a coherent spatial spectral estimation 
technique, Chotiros (1995) proposed and used the following intensity speed and 
direction finding algorithm. This velocity finding algorithm assumes a direction 
and speed for a plane front of propagating energy, determines the corresponding 
delays for each receiver with respect to a reference receiver, and then adds 
normalized and delayed versions of all receiver intensities. If the position of the 
wth receiver is denoted by r m and the reference receiver by rj, the relative signal 

delay at the mth receiver with respect to the reference receiver is given by 

A 

^m = ^-2 ' r m ! c 2 


(4.11) 
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A 

for a plane wave traveling in a direction k 2 and speed c 2 . For a given speed and 
propagation direction, the normalized received signals are delayed by Eq. (4.11), 
and then summed. The resulting magnitude is plotted as a function of grazing 
angle (depression angle) and wave speed values. A peak in the resulting plot 
occurs when the assumed direction and speed coincide with the actual intensity 
front. As can be deduced from section 4.4, this type of incoherent processing 
cannot unambiguously resolve wave velocity. In addition, the intensity front 
measurement will be off if the flat surface coherent component of the signal 
dominates at the reference receiver but incoherent diffracted energy dominates 
at the other receivers. 

In order to determine the direction of the intensity front, an array is generally 
required to span 3-D space (see Appendix D). However, the direction of the 
intensity front is assumed to have no y component, and the receivers are 

therefore only required to span the xz plane. 

4.5.5 Results 

Speed-angle distributions for the five source locations for 0,- = 79°, 40°, 

24.9°, 127°, and 9.7° are shown in Figs. 4.11a, 4.11b, 4.11c, 4.1 Id, and 4.11e, 
respectively. Chotiros (1991) states that the receiver at z=-0.05m is usually 
taken to be the reference receiver. These plots are therefore obtained using the 
receiver at z = -0.05m for the reference. Figures 4.1 la and 4.11b correspond to 
incidence above the critical grazing angle, and show energy propagating 
approximately at the sediment sound speed and at angles given approximately 
by Snell’s Law. The curve in the figures obeys Snell's Law. Since the incident 
grazing angle is above critical in Figs. 4.11a and 4.11b, almost all of the energy 
penetrating into the sediment is due to the zero-order refracted energy. The 
coherent intensity due to the flat surface solution is also greater than the 
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incoherent intensity in Fig. 4.11c, but the first-order incoherent intensity is of 
comparable magnitude. At sub-critical angles (Figs 4.1 Id, 4.11e), the energy 
appears to be propagating at about 1400m/s, at angles that also agree with 
Snell’s Law for this apparent slow wave speed. This model does not include 
slow wave propagation; this apparent slow wave propagation is explained by 
diffractive scattering. Equation (4.6) shows that the diffracted intensity front will 
always appear to obey Snell's Law. 

Although Chotiros reports apparent slow wave propagation at speeds 
between 1 lOOm/s and about 1350m/s (see Fig. 6, Chotiros, 1995), we observe 
an apparent slow wave speed of 1400m/s in Figs 4.11c, 4.1 Id, and 4.11e. 
However, different choices of roughness spectrum will result in different apparent 
slow wave speeds. Although the roughness spectrum used in the simulation is 
suitable for the ocean experiment location, burying receivers into the sediment 
will likely alter the roughness spectrum at the interface above the receiver. 
Therefore, it is possible that the roughness spectrum used in the simulation is not 
an accurate representation of the seafloor at the experiment site. The effect of 
decreasing the roughness by changing the filter parameter in Eq. 3.66 is shown 
in section 4.5.7. In the next section, another factor that can affect the apparent 
slow wave speed is discussed. 

4.5.6 Significance of reference receiver choice 

The apparent slow wave speed as determined by Chotiros' sound speed and 
direction finding algorithm can be affected by reference receiver choice. 
Refraction is dominant when the incident angle is above critical, and diffraction 
can be dominant at subcritical incidence, depending on the receiver depth and 
incident grazing angle. Since the refracted path is much shorter for the receivers 
closest to the interface, the refracted coherent intensity will sometimes dominate 
over the incoherent intensity only at shallow receivers. In fact, the zero-order 
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field is negligible for the transverse distance of 19.19m between source and 

receiver (incident grazing = 12.7°) at all receiver depths except for the depth of 
0.05 m — the reference receiver. Figures 4.12a, and 4.12b show the mean 
square incoherent field and square of coherent field magnitude (1 Pa at 1 m from 
source in water) for all receivers in the simulation as a function of time for the 

12.7° incidence grazing example. For all receivers except the reference receiver, 
the mean square incoherent field is much greater than the square of the coherent 
field magnitude; at the reference receiver, the coherent field is much greater than 
the incoherent field. Since the refracted path is the fastest, the peak in the 
incoherent intensity pulse will always arrive after the zero-order component of the 
coherent pulse. As a result, the relative time delay is greater than would be the 
case if another receiver were chosen as the reference receiver and the receiver 
signal at z - -0.05m not included in the intensity velocity algorithm. The apparent 
slow wave speed using this reference receiver is therefore slightly slower than 
the apparent slow wave speed using any other receiver. For both incidence at 

12.7° and 9.7°, the first-order incoherent intensity is less than the zero-order 
coherent intensity for the receiver at depth 0.05m; for all other receivers, the 

opposite is true. The outputs for 12.7° incidence and 9.7° incidence are plotted in 
Figs. 4.13a and 4.13b using the receiver at z = -0.18m as the reference receiver 
and removing the receiver at z = — 0.05m. These simulation results show slightly 
faster apparent slow wave propagation (1500m/s) than the results in Figs. 4.1 Id 
and 4.11e. Note that this speed is significantly smaller than the sediment sound 
speed of 1729m/s. 

Note that if the apparent slow wave in the experiments is also due to 
diffraction, then using any other receiver as the reference than the one at 
z = -0.05m in the experiment would also result in observation of an apparent slow 
wave that is faster than would otherwise be observed. It is important to note that 
the slower the observed apparent slow wave, the greater is this discrepancy. 
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4.5.7 Effect of the filter parameter 

The purpose of this section is to show the effect of the filter parameter a (see 
Eq. 3.66) on the apparent slow wave speed using Chotiros' speed and direction 
finding algorithm. In Figs 4.11 and 4.13, the filter parameter « = 05m, resulting 
in kih = 05m, was chosen. Figures 4.14a and 4.14b show speed-angle 

distributions for the two source locations corresponding to 0 ; - =12.7° and 9.7°, 
respectively, for a filter parameter value of a = 2.0 m — k\h = 1.96. As shown by 
Thorsos (1996a), perturbation theory is valid for this case. As in Fig. 4.13, the 
receiver at the depth of 0.18 m is used as the reference receiver, and the receiver 
at a depth of 0.05m is discarded. The differences between Figs. 4.14a and 4.13a 
and Figs. 4.14b and 4.13b are negligible. Increasing the filter parameter beyond 
a - 05 m to a value more than twice as large as the scale of the receiving array 
has negligible effect on the speed-angle distribution output. 

Figs. 4.15a and 4.15b use a filter parameter a = 0.02m, resulting in an 
apparent slow wave speed of about 1200 m/s — the same as the tank experiment 
simulation results. Since Chotiros' sound speed and direction finding algorithm 
measures the propagation of the intensity front, it does not determine the actual 
direction of propagation. The reason that this algorithm yields a direction and 
speed corresponding to an apparent refracted slow wave follows from the 
discussion in section 4.4 (Thorsos, 1995). 

4.6 Summary 

By modeling sand as a fluid, and including a small amount of roughness 
h rms = 1, * he tank simulation results match the acoustical penetration 

experimental results obtained by Boyle and Chotiros (1992), both in magnitude 
and in arrival time of the intensity pulse. The simulation results of this tank 
experiment show that roughness must be included in a sediment penetration 
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model, and that even a small amount of roughness can account for energy 
penetrating the interface measured in their experiments. 

Although the measured results in the ocean array experiments (Chotiros, 
1995) yield an apparent slow wave that is between 50 m/s and 300 m/s slower 
than the results obtained from the simulation, the assumed roughness spectrum 
may not be an appropriate representation of the experimental site after the 
receivers were buried. However, the results of this chapter show that diffraction 
results in apparent slow wave propagation that appears to obey Snell's Law. 
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Figure 4.1 : Scattering model diagram 
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Figure 4.3 : Roughness spectrum used in tank simulation, a = 0.01m, y = 3.0 
w 2 = 2.17 x 10 -5 m , h = 1mm. 
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First-order mean square incoherent field 

xio -4 



arrival time (ms) 

use g(t) -> u(t) 
use g(t) 


Figure 4.4a : Time-dependent mean square incoherent field determined in tank 
experiment simulation; solid line uses Eq. (2.85), dotted line uses Eq. (4.1). 
Roughness parameters are given in Fig. 4.3, frequency = 60kHz, p = 2.0, 6 = 0.015, 
v= 1.125, c\ - 1489m/s, = 33.33|j.s. Source height = 0.49m, receiver depth 

= 0.08m, transverse distance between source, receiver = 2m. 
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First-order mean square incoherent field 



1 2 3 

arrival time (ms) 


result with u(t 
result with g(t 


Figure 4.4b : Time-dependent mean square incoherent field determined in tank 
experiment simulation; solid line uses Eq. (2.85), dotted line uses Eq. (4.1). Receiver 
depth = 0.33 m, otherwise same parameters as Fig. 4.4a. 
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Receiver Depth = 0,08 m 
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Figure 4.5a : First-order incoherent intensity as a function of time and transverse 
distance between source and receiver. Arrival time is at pulse peak; black curve 
corresponds to i200m/s arrival time of apparent slow wave. Same parameters as Figs 
4.3, and 4.4 -- frequency = 60kHz, p = 2.0, 5 = 0.015, v = 1.125, Cj= 1489 m/s, 
t s = 33.33 ps. Source height = 0.49 m, receiver depth = 0.08 m. 
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Receiver Depth = 0.33 m 



horizontal distance from source (m) 


Figure 4.5b : First-order incoherent intensity as a function of time and transverse 
distance between source and receiver. Receiver depth = 0.33 m, otherwise same 
parameters as Fig. 4.5a. 
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Receiver Depth = 0,614 m 




Figure 4.5c : First-order incoherent intensity as a function of time and transverse 
distance between source and receiver. Receiver depth = 0.614m, otherwise same 
parameters as Fig' 4.5a. 
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Apparent slow wave speed as a 
function of scattered grazing angle 



Figure 4.7 : Dependence of the apparent slow wave velocity on scattered grazing angle 
for incident grazing angle, 0, = 10°. 
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Apparent refracted slow wave 
speed vs. incident grazing angle 



incident grazing angle (degrees) 

c 2 = compressional wave speed in sediment 


Figure 4.8a : Dependence of the apparent slow wave velocity on incident grazing angle 
assuming average scattered field is straight down. 
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Apparent refracted grazing 
angle vs. incident grazing angle 



Incident grazing angle e 1 (degrees) 


Figure 4.8b : Dependence of the apparent refracted grazing angle on incident grazing 
angle assuming average scattered field is straight down. 
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Figure 4.9 : Ocean experiment simulation layout. 
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First-order mean square incoherent 
pressure at 5 cm below interface 
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Figure 4.10a : Time-dependent mean square incoherent pressure in ocean experiment 
simulation for receiver depth = 0.05 m. Solid line uses Eq. (2.85), dotted line uses Eq. 
(4.1). Transverse distance between source and receiver = 19.19m, a = 0.5m, y = 3.0, 

w 2 = 6.2 x 10 -5 m, h = 1.2cm, frequency = 20kHz, p = 1.9608, 8 = 0.0163, v = 1.1257, 
C\ = 1536m/s, t s = 100 (is. 
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First-order mean square incoherent 
pressure at 18 cm below interface 



Eq. (4.1) 
Eg- (2.85) 
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Figure 4.10b : Time-dependent mean square incoherent pressure in ocean experiment 
simulation for receiver depth = 0.18 m, otherwise same parameters as Fig. 4.10a. 
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Incident Grazing Angle = 79.0 (degrees) 



sound speed (m/s) 


Figure 4.11a : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Velocity distribution for 0 i =79° (0.43 m 

transverse distance between source and reference receiver), same parameters as Fig. 
4.10. White curve corresponds to energy propagating at speeds and directions 
coinciding with Snell's law. 
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Incident Grazing Angle = 40.0 (degrees) 



sound speed (m/s) 


Figure 4.11b : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Velocity distribution for 9/ =40° (4.91m 

transverse distance between source and reference receiver), otherwise same 
parameters as Fig. 4.11a. White curve corresponds to energy propagating at speeds 
and directions coinciding with Snell's law. 
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Figure 4.11c : Normalized output intensity of ocean simulation using intensity velocity 

finding algorithm given in Chotiros (1995). Velocity distribution for 0 y =24.9° (9.22m 
transverse distance between source and reference receiver), otherwise same 
parameters as Fig. 4.1 la. White curve corresponds to energy propagating at speeds 
and directions coinciding with Snell's law. 
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Figure 4.1 Id : Normalized output intensity of ocean simulation using intensity velocity 

finding algorithm given in Chotiros (1995). Velocity distribution for 0 Z - =12.7° (19.19 

transverse distance between source and reference receiver), otherwise same 
parameters as Fig. 4.11a. White curve corresponds to energy propagating at speeds 
and directions coinciding with Snell's law. 
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Incident Grazing Angle = 9.7 (degrees) 



sound speed (m/s) 


Figure 4.11e : Normalized output intensity of ocean simulation using intensity velocity 

finding algorithm given in Chotiros (1995). Velocity distribution for 0,- =9.7° (25.83 m 
transverse distance between source and reference receiver), otherwise same 
parameters as Fig. 4.11a. White curve corresponds to energy propagating at speeds 
and directions coinciding with Snell's law. 
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Figure 4.12a : Time dependent mean square incoherent pressure in ocean simulation 
at each receiver for the 12.7° incidence grazing (Fig. 4.1 Id example). 
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Square of zero-order coherent pressure magnitude 
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Figure 4.12b : Square of coherent pressure magnitude at each receiver in ocean 
simulation for the 2.7° incidence grazing (Fig. 4.11e example). 
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Figure 4.13a : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Same parameters as Fig 4.1 Id except 
receiver at z- -0.05m is removed, and receiver at 0.18 m depth (see Fig. 4.9) is used 
as reference receiver. 
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Figure 4.13b : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Same parameters as Fig 4.11e except 
receiver at z = -0.05m is removed, and receiver at 0.18 m depth (see Fig. 4.9) is used 
as reference receiver. 
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Figure 4.14a : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Same parameters as Fig 4.13a except 
roughness spectrum filter parameter a - 2.0 m, rather than a = 0.5m. 
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Figure 4.14b : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Same parameters as Fig 4.13b except 
roughness spectrum filter parameter a - 2.0 m, rather than a — 0.5 m. 







Incident Grazing Angle = 12.7 (degrees) 


item: 


• -■M&M 
,*XyiwT 
■m^ms 






sound speed (m/s) 


Figure 4.15a : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Same parameters as Fig 4.13a except 
roughness spectrum filter parameter a = 0.02 m, rather than a = 0.5 m. 
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Figure 4.15b : Normalized output intensity of ocean simulation using intensity velocity 
finding algorithm given in Chotiros (1995). Same parameters as Fig 4.13b except 
roughness spectrum filter parameter a = 0.02 m, rather than a = 0.5 m. 









CHAPTER 5 


Perturbation derivation including the effect of gradients 


5.1 Introduction to chapter 5 

Solution to the first-order scattered field from and through a rough interface 
separating two homogenous fluids is derived in chapter 3. The simulations of the 
tank and ocean experiments in chapter 4 use first-order perturbation theory to 
calculate the incoherent scattered intensity, and the flat-surface (zero-order) 
results are used to calculate the coherent intensity. Although first-order 
perturbation theory results allow solution for the incoherent intensity proportional 
to the second power of k x h, the zero-order solution results in computation of the 

coherent intensity only to {k x hf. In order to calculate the coherent intensity 

rs 

component proportional to (k^h) , the second-order scattered field is required. 

The second-order field and second-order coherent field are calculated in this 
chapter, and compared with the zero-order component of the coherent field. 

Since the sediment in both the tank and ocean experiment simulations is 
modeled as a homogenous fluid, the effect of sound-speed gradients on the 
scattered field is not included. In this chapter, the sound-speed gradients are 
included in the model. Moe and Jackson (1994a) take an approach suggested 
by Winebrenner et al. (1995) that allows for layering below the rough interface. 
Although the medium just below the rough interface is constrained to be a 
homogenous lossy fluid, the medium below the lowest point on the interface is 
constrained only to be vertically stratified (see Fig. 5.1), and its complexity is 
contained in the reflection coefficient for the corresponding flat-surface problem. 
This approach is used to derive recursion equations that allow for calculation of 
the Nth-order scattered field. The zero-order field and first-order scattered field 
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are derived and compared with the results in chapter 3. The second-order field is 
derived in order to obtain an expression for the second-order component of the 
coherent field. The second-order coherent field is compared to the flat surface 
transmitted field. The effect of including sound-speed gradients as well as a 
shear supporting medium is also discussed. 


5.2 Scattering problem 


The seafloor is represented as a two dimensional random rough fluid-fluid 
interface. The geometry of the scattering problem is shown in Fig. 5.1. A 
pressure field, \|/,(r), is incident on a 2-D rough interface separating the 

homogenous lossless fluid in medium 1 (water) from a homogenous lossy fluid in 
medium 2 (sediment). Unlike the problem in chapter 3, the region below the 
lowest point on the rough interface (z = -d) is arbitrary. As in chapter 3, the total 
field above the highest point of the rough interface in medium 1 is a sum of an 
arbitrary incident field, \|/,(r), and the resulting scattered pressure field, V/(r), 

and is expressed as 


\j/j(r) = J J 2 A:k(K)e" l7:iPl(K)z + ^(IQe 


i^iPi(K)zl /K R 


(5.1) 


A factor of e iat is suppressed. Because of the flat interface between medium 2 
and medium 3, the field in medium 2, \|/ 2 , is expressed as a sum of plane waves 
traveling in the positive z direction, as well as a sum of a plane waves traveling 
in the negative z direction: 


V2 ( r ) = J d 2 K 4V- (K) e _I * 2p2 (K) z + ^ (K) e '* 2p2 (K) Z 


iKR 


(5.2) 


Note that in the absence of upward refracted or reflected energy, 'P 2+ = 0. The 
derivation proceeds similarly to the derivation in chapter 3. In Eq. (5.1), 'Fy is 

unknown; in Eq. (5.2), both ¥ 2 - and ¥ 2 + are unknown. Three equations are 
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required to solve for these three unknowns. Two of the boundary conditions are 
the continuity of pressure on the interface, given in Eq. (2.2), and the continuity 
of normal velocity on the interface, given in Eq. (2.4). These boundary conditions 
are used in the scattered field derivation for the fluid-fluid model in chapter 3. An 
additional equation is supplied by a technique originated by Winebrenner et al. 
(1995) in their treatment of electromagnetic scattering from rough dielectric 
surfaces. This third equation uses the general boundary condition that relates 
v f / 2 + to at z--d in terms of the reflection coefficient r 23 at z = -d for a wave 

incident from medium 2 to medium 3 (Bucker, 1970): 

TWK) = ¥2±W e -2ik 2 MK)d' (5.3) 

23V ' ^ 2 _(K) 

This equation is key to Winebrenner's approach. All of the complexity of medium 
3 is contained in r 23 . Note that media 2 and 3 can both be lossy, with a complex 

wave number, but the sediment in region 3 is not constrained to be a fluid, e.g., it 
could be viscoelastic or porous. Furthermore, arbitrary gradients in density and 
wave number are allowed in medium 3 provided variations are only in z. 

The goal of this chapter is to find the scattered field in terms of the incident 
field. Since 

47 (K) = JA'7i 1 (K,K') v T ( (K / ) (5.4) 

and 


^.(K) = J At 12 (K,K')'F 1 (K') (5.5) 

(Eq. 3.5 and Eq. 3.6), the scattering problem is solved when the T-matrices T u , 
and T n are found. Perturbation theory solves this problem by expanding the 
fields and each T-matrix in terms of a power series representation, 
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to= y!M!!tI")(K) 

iS n! 

T(K,K,') = Y®2.r < ">(K,K / ), 

£S n! 


(5.6) 

(5.7) 


and solving one order at a time. Both the factorial and the factor ( k\h)" do not 
need to be displayed, but this representation emphasizes that the perturbation 
expansion can be viewed as a power series expansion. 


5.3 Derivation of Nth-order recursion relation 

Although only the zero, first, and second-order fields are derived in this 
chapter, the results of this section can be used as a starting point to derive 
higher-order results. An immediate application of higher-order results would be 
to use the second-order field component as a correction to the first-order 
coherent field. Higher-order results serve as both a correction and measure of 
accuracy of lower-order results. 

5.3.1 Background 

Define the Fourier transform of the wth power of the surface, 

F„(K) = —U- f d 2 Kf n (R)e~ iKR , (5.8) 

nV ’ (2tc) 2J 

so that F„(K) is obtained from a convolution of F„_,(K) with F( K), 

F„(K) = F„_ 1 (K)0F(K). (5-9) 

Also, 

Fq(K) = 8(K), (5-10) 


and 
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F,(K)sF(K), (5.11) 

is the Fourier transform of the surface profile function. Taking the gradient (Eq. 
3.4) of 

/ n+1 (R) = (AF„ +1 (K)e' M (5.12) 


(KR 


and then dividing by n+1 yields 

=- i - J d 2 K KF n+1 (K)e 
n +1 n +1 J 

The following identity, 


/"(R)V/(R) = -i-f d 2 K KF„ +1 (K)e 

n +1 J 

follows from Eq. (5.13) and 

V/" +1 (R)=(« + l)/"(R)V/(R). 


(KR 


(5.13) 

(5.14) 

(5.15) 


5.3.2 Continuity of pressure 

Invoking the Rayleigh hypothesis, we use Eqs. (5.1) and (5.2) in Eq. (2.3) to 
obtain the equation of continuity of pressure: 

= |Ak.(K)e-^ < K >« R >+ f 2+ (K)e“* < W<*>] e ® * (5 A 6) 


Expanding the exponential terms in a power series, 

e rt 1 P, (K)A/(R) . 2(, Pl(K)/(R) )»MF, 


(5.17) 


n =0 


and using Eqs. (5.3), (5.12), and (5.13) yields 
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J d 2 K £(iP 1 (K))"^^-['P / (K)(-l) B + 'F / (K)]J A' F„(KV (K+K ' )R 


n =0 

= fd 2 K £(/p 2 (K))"^^-^ 2 -(K){[(-l) B + r 23 e /2 * 2p2(K)d 

TV~0 

x J d 2 K'F n (K') e /(K+K ' ) ' R }. 


(5.18) 


After using the operator — \d 2 R 

(2tc) 2 J 

over R, and over K, to obtain 


2 n e -/K R 


, the above expression is integrated 


f d 2 K X(iP,(K))"^-[ , P j (K)(-l)" +l' / (K)]jv,(K"-K) 

n=0 "• 


= JA £(iP 2 (K))"M->P 2 .(K)[(-l)'' 

n =0 


+ r 23 e«* 


F n (K"-K)' (5 . 19) 


Expanding fields in Eq. (5.19) in a power series (Eq. 5.6) yields 


jd 2 K X(ipi(K))" fe?i->F i (K)(-l)' , F„(K"- K) 
+\d 2 K f, (>P 1 (K)rM-£^-'P}' , (K)F m (K" 


m =0 

oo 


-K) 


=J A £ (ip 2 <K))4M 

m=0 { 


/=0 

m 


+ r 23 e 


*2* 2 p 2 (K)rf 


) 


F m (K" - K) 


(5.20) 


Combining powers of k\h results in the following expression for the continuity of 
pressure on the interface: 
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^ ) (K")(y /, +i) 

= jd 2 K(-l) N {iV l (K)) N F N (K"-KW i (K) 

+jd 2 K X (m) (*'Pi(K)) m F m (K"-K)^" m) (K) 

m -1 

~\d 2 K £ (^)(iKP 2 <K))"F w (K”-K)[Y + (-l)’"]T<!'-'” ) (K) 
m =1 

hM^K"), 

7 (5.21) 

where 

Y = r 23 (K)e i2 ^ 2(KM . (5.22) 


5.3.3 Continuity of normal velocity 


An expression for the equation for continuity of normal velocity on the 
interface is obtained by substituting Eq. (3.3), 


n = 


z - V¥(R) 

£ - Vhf( R)| ’ 


into Eq. (2.4), 


!*.V V i(r)-—*-Vy 2 (r)| 

Pi P2 Z J[ ) 


Pl p 2 --'"z-JW 

yielding an expression for the continuity of normal velocity, 


(5.23) 


(5.24) 


-AVy, -“i 

Pi Z =/(R) Pi W(R)L -U/(R) 


(5.25) 


Using Eqs. (5.1) and (5.2) in Eq. (5.25) results in the following expression: 
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-1 \d 2 K k x p l( K)[-^(K)e-^P> ( W (R ) + 

_ JL Jj 2 £ k • V(/j/(R)) Y^K)e -, * lPl(K) ^ (R) + ¥ / (K)e'* l|3,(K)/l/(R) e* K R 

= — jd 2 K /t 2 p 2 (K)[-^ 2 -(K)e" , ' <:2p2(K)/!/(R) + 'F 2+ (K)e^ 2P2(K)¥(R) ]e i ' K ’ R 

- — [ d 2 K K • V(/t/(R))[ V P 2 _(K)e -1 * 2 ^ 2(K) ^ (R) + 'F 2+ (K)e^ 2P2(K)/i/(R) |e I ' K ' R , 

P2 L (5.26) 

which becomes 

(5.27) 

Pi J n= 0 „ 

| d 2 K K • V(/*/(R)) £(/p!(K)/(R)) n ['F i (K)(-l) w + ^(KjJe®* 

Pi n=0 " - 

= — f ^ 2 AT fc 2 p 2 (K) 2(iP2(K)/(R))" ^-[-'F 2 _(K)(-l) n + 4 / 2+ (K)]e fKR 

P2 n=0 "• 

f d 2 K K • V(/j/(R)) Y (ip 2 (K)/(R)) B fef-[T 2 _(K)(-l) n + 'P 2+ (K)]e iK ‘ R 

P2 J n=0 K ’ 

with the use of Eqs. (5.6) and (5.17). Using Eqs. (5.3), (5.12), and (5.14) in the 
above expression yields 
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— J d 2 K kfaiK) ^(/PjCK))" M-[-(-l) n ^(K) + 'Fy(K)] (5.28) 

Pi n =0 

xfd 2 K'F„( K')e''< K+K ’> R 

J A). X(<-U.(K)r^-h(K)(-l)» + Vj (K)] 

Pi n=0 "• 

x — f J 2 /*:' K • K' F n+l (K') e z ( K+K 0 R 
n + l 3 

= ± J A* 2 p 2 (K)£(<P 2 (K))" Ml>P 2 _(K)[-(-l)” + r 23 e i2t *(KV 

P2 n=0 ”• 

x J A'F„(K')e'< K+K '> R 

f A* Y(ip 2 (K)Y , &^>F 2 _(K)[(-l)" +r 23 e i2 '* <K >‘'' 

p2 J n= Q w! 

x— f A'K • K' F n+ \ (K') e^ K+K ^ R . 
n + l 3 

Equation (5.28) is integrated over R using —f d 2 Rc~ lK R , and over K. The 

(2jc) j 

substitution in Eq. (5.22) is made, yielding the equivalent expression, 

- 1 d 2 K k^(K) XO'Pi(K))" M-[-(-l)”^(K) + ^(K)]F n (K" - K) (5.29) 

Pi n=0 ”• 

--LJ d 2 Kh X(iP i(K)) n ^-[^(K)(-l) n + V f (K)] K ' { ^~ K) F n+l (K" - K) 

Pi n=0 n ‘ 

= —\d 2 K i 2 p 2 (K) £(ip 2 (K))“ ^->F 2 _(K)[-(-l)” + y]f„(K" - K) 

P2 n =0 "• 

~j d 2 Kh £(iP 2 (K)f ^->F 2 .(K)[(-1)" + Y] K( ^ K) WK~ - K). 


Dividing the above expression by k h multiplying by p lf and expanding in a power 
series — Eq. (5.6) — yields 
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JAp,(K)X(iPl(K))"^- 

n=0 n - 


-(-1)"'P I (K) + 

1=0 L 


oo 


-ijd 2 K XO'Pl(K)) 


n (k\ h ) 


«+l 


AZ—0 


(« + 1 )! 


x 


K • (K" - K) 


'P i (K)(-l)" + £^'P}' ) (K) 

/=o /- 

= ij Akp 2 (K) £(ip 2 (K)) n ML 

P az=0 

<0(K)[-(-l) n +y]f„(K"-K) 

1 A t • 


F W+1 (K"-K) 


/=0 


oo 


n (ty) ^ \p(/) ^ 


x[( _ ir+Y ]WE) 


F„ +1 (K"-K). 


4(K 


-K) 


(5.30) 


Separating terms, renumbering indices (setting n = l+m), and equating like 
powers of k\h yields 
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- J A Pi (K) ( 1} (K)F n (K" - K) 

J Ni 

+J d 2 K p! (K) X . ^~ m) (K)F w (K" - K) 


ra!(N — m)! 

+ fd 2 X (- 1 ) A 'OPl(K)) A ' y (K) K ■ < K " - K > f „(K" - K) 

J iV! Pl(K )*, 2 

_[A y (WK))\ w-,») K.(K"-K) 

J £;«**-»)■' P.W 2 m( ’ 

= - J A Kp 2 (K) X ^ <K ^7 T^-'” ) (K)[-(-l)'” + 7]F m (K" - K) 
P m -0 

-±ja x ( '^ <K)> > g'" ) < K )H- i ) m + y] k r ( ^"~ 2 ) ^( k " 

P j ^m\(N-m)\ 1 J kP 2 (K)^ 

Corresponding to the definition in Eq. (5.22), define 


—^F,v(K" - K) 

Pi (K)fcf 


-K). 


y" = r 23 (K")e 


i2k$ 2 (K")d 


The following definitions, 

AN J l hS^L . FN(K ''. Kh 

/p 1 (K)F(K"-K) 

M^(K'-K) 

2! 2 

A(K",K) = 


(5.31) 


(5.32) 


(5.33) 


(5.34) 


0‘P i(K))‘ 


Fm (K" - K) 
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/Kp 2 (K)F(K" - K) 
(*k| 3 2 (K)) ^ K // _ K -j 


B + (K", K) = 


(iK(3 2 (K)) i 


fy(K"-K) 


-/k| 3 2 (K)F(K" - K) 
(~ngp 2 (K))_ _ 

2 ! 


B_ (K", K) = 


(~/kP 2 (K)) 


Fm (K" - K) 


UN-l) y(N 2) 

Yf(K) = — - --—- ... 'i’f 

f{ ; (JV-1)! (N- 2)! f 


and 


Yl_(K) = 


—1) W (N 2) 
(N - 1)! (AT-2)! 


... 


result in more concise forms of Eq. (5.21), 

_2-= (v +1)—-- 

Nl y ’ N\ 


J d 2 K [yBl (K", K) + B‘_ (K", K)]y 2 _(K) 

Jrf 2 F(-l) A, A Ar (K // ,K)^(K) 

fd 2 KA t (K",K)Y / (K), 


(5.35) 


(5.36) 


(5.37) 


(5.38) 


(5.39) 


and Eq. (5.31), 
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Pi cm 


'py v) (K") 

TV! 


= J d 2 K Pj (K)(-1)' N A n (K'\ K) ^ (K) 
-u(N-\)\d 2 K K ' (K "^ (-l^A^CK^K^-CK) 

Pi (K)&i 

+ f d 2 K K ’ (K "~ K) - A r (K",K)Yf(K) 

J Pi (K)fcf 7 

-J A PjCk^ck'.kjy/CK) 


i J A Kp 2 (K)[^(K",K) - BL(K",K)]y 2 _(K) 
V(t"-1)kE 2 (K") 2~ n ) 


- f d l K K (K " !P f'yB+(K ,, ,K) -BL(K",K)1 y 2 _(K) 

p J Kp 2 (K)fcf 1 J 


(5.40) 


Adding Eq. (5.39) multiplied by (y"-1)kP 2 (K"), and Eq. (5.40) multiplied by 
p(y" + l), eliminates and results in an expression for the Nth-order 

scattered field: 
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vpW/'K"') 

-i-T-- -[p(t" + l)Pl (K") - (y "- l)Kp 2 (K")] 

N\ 

= p(y" + 1)J A P 1 (K)(-1) A 'a a ,(K",K)'F i (K) 

-p(y" + 1)J d 2 K p 1 (K)V(K",K)Y / (K) 

-u{N - l)p(y" +1) f A K • (K "-f l(-lf AN (K",K)%(K) 
PV P,(K)jt, 2 

+p(Y" +1) f A K-(K// ~ 2 K U f (K // ,K)Y / (K) 

V JJ pj(K )kf J 

+( Y " + 1)J d 2 K k|3 2 (K)[^(K",K) - BL(K",K)]Y 2 _(K) 

-(y"+ 1)f A K (K " ~ [yb+(K",k) - bL(K",k)1y 2 _(K) 

v kP 2 (K.)*? *• J 

+(y" - i)kP 2 (K")| A (-if a n (K",k) %(K) 

+(y" - 1)kP 2 (K")| A A f (K",K)Yy(K) 

- (y" - 1)kP 2 (K") J d 2 K (K", K) + Bi (K", K)] Y 2 _ (K). 


(5.41) 


Define 



pP^kVi+y'^ 

Kp 2 (K")ll-Y"J 


i + r(K w ) 
l - r(K") ’ 


(5.42) 


where r in second equality in Eq. (5.42) is the total reflection coefficient of the 
corresponding flat surface problem at the zero-mean rough interface, obtained 
from the well known reflection coefficient expression (see for example Sec. 3.3 of 
Brekhovskikh and Lysanov). 



r 12 (K) + r 23 (K)e‘ 2 * 2p2(K)t/ 

1 + TuiKfoiiK)^ 2 ^ 2 ^ ’ 


(5.43) 


where 

r rT n PPi(K)-Kp 2 (K) 
l2{ } P Pi(K) + Kp 2 (K) 


(5.44) 
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is the reflection coefficient corresponding to a flat surface fluid-fluid interface 
between medium 1 and medium 2 in the limit as the width of medium 2 , d , 
becomes arbitrarily large, or medium 3 parameters are identical to medium 2 
parameters. The term r 2 3 is a general reflection coefficient corresponding to the 
medium 2 - medium 3 flat interface. Although medium 2 is constrained to be a 
fluid, medium 3 is only constrained to vary in the z direction. Dividing Eq. (5.41) 
by (y"-1)kP 2 and using Eq. (5.42) yields 


T (W)(K") 

te"+1] 


N\ 


+jd 2 K 


(~l) N A N (K",K)'¥ i (K) 


jTPl(K) l ^"K • (K" - K) 

Lpkk") fronPi (K)^ 2 j 

_- 14 . ' (K . . _ 1 1 A f (K",K)Yy(K) 


(5.45) 


Pl(K") 

^"kP 2 (K) 


l"K ■ (K" - K) 


( d 2 K _ 

J [ pPKK") P p 1 (K")Kp 2 (K)fc 1 2 

j J 2 tf[yBt(K",K) + BL(K",K)]y 2 _(K). 


[^(K",K)-BL(K",K)]y 2 _(K) 


Dividing the above expression by k\ , multiplying by pj, substituting Eq (5.33) 
through Eq. (5.38) into the above expression, using the binomial coefficient 
notation, 


(N\ 


V 



N\ 

(N — m)\m\ ’ 


(5.46) 


and collecting terms in the same power in k\h yields 
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(5.47) 


Apply Eq. (5.42) to Eq. (5.47), 
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'py v) (K") 


= }d 2 K 


i+r(K") pi(K) i - r(K") i+r(K") k -(K"-k) 


2 Pi(K") 

/ *p 


P^K'OPitK)^ 


x(-l f (^(K))" F n (K" - K)^(K) 
+\d 2 K 


1 + r(K") p^K) 1 - T(K") ( l + r(K") K • (K" - K) 


Pi (K") 


2 

x X (^)^K)) m F m (K" - K)^- w) (K) 

m -1 


2 p^k'OP^k^j _ 


+jd 2 K 
N 


1 + r(K") Kp 2 (K) 1 + r(K") K • (K" - K) 


pPi(K") 


(5.48) 


x 


2 pp 1 (K")Kp 2 (K)^J 

Z (m)'KP2(K)rF m (K" - K)[y - (-irW"-”) 


(K) 


m=1 


TV 


+ J A X 


aA i - r(K "),. 


(iKp 2 (K)) m F m (K"-K)[y + (-l) 


m 


'F 


(W-m) 


2 - 


(K). 


/n=l 


Equation (5.48) can be further simplified by using the identities 


K-(K"-K) K K" 


*? 


-l+Pl(K) 


(5.49) 


and 


(Pi(K)) -(kP 2 (K)) = 


2 _kf — K‘ 


kj-K 2 

*? 


= 1-K 


(5.50) 


to give 
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^(K") 

4 |A( 1 + r ( o) pi(K „ )Pi(K) 

x(-l) A '(iP,(K)) A, F A ,(K"-K)'P i (K) 




( 1 +r(K")) 


4 J-* +r(K,,) )p4® [4r - 1 - m <K) Si. 

N / \ /.. v 


XI f m)(>Pl (K)f F m < K " - K)^"'~ m) (K) 


m-\ 


+ 


-Jj 2 ^(i+r(K"))- 


i 


P p 1 (K")kP 2 (K) _ 


2 K-K" 

K - T~ 

*? 




xz(«)^(K)r^-K)[T-(-ir^m 
m=i ^ 

+i J A (1 - r(K")) £ (^(ilcp 2 (K)) m F m (K'' - K)[y + (-1) 

m=l 


m 


^ m) (K). 


(5.51) 


The Nth-order fields scattered away from the interface and through the interface 
are found from the recursion equations, Eq. (5.21) and Eq. (5.51). Each order 
field is found in terms of all lower-order fields. Although reciprocity is not obvious 
in Eqs. (5.21) and (5.51), nevertheless, it should be possible to show they are 
reciprocal when each order is calculated. 


5.4 Zero-order 


The flat surface reflected field follows from Eq. (5.51) by setting N = 0. 

4f <K") - i J AW - KW® j&S 


X 


—5— + 1 _ Pi (K^Pi (K) 
kf 


Pi(K")Pi(K) 

(l-rgn) 

(1 + T(K")) 


(5.52) 


¥P\k") = J d 2 KS( K" - K)r(K") v P ( (K) 


(5.53) 


or 
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^ 0) (K) = r(K)^-(K). (5.54) 

This result is the same as the zero-order result obtained in chapter 3, except that 
the reflection coefficient in this result, given by Eq. (5.43), includes medium 3 in 
the partial reflection coefficient, r 23 , at the flat interface between medium 2 and 

medium 3. The zero-order T-matrix, 1§\ is the flat surface reflection coefficient 
multiplied by a delta function; 

T§\k",K) = r(K)5(K" - K), (5.55) 

where r is given in Eq. (5.43), representing the total reflection coefficient at z = 0 
when the rough interface is replaced by its mean. The zero-order field 
transmitted through the interface follows from Eq. (5.21) by setting N = 0, 

^ ) (K")(y" +1) = J d 2 Kb( K" - K)^(K) + ¥f\K") (5.56) 

or 

^°}(K) = T(K) x ¥ i (K), (5.57) 

where 

T(k) = ^ + (5.58) 

v ’ 1 + y 

is the flat surface transmission coefficient. Note that this transmission coefficient 
relates the down-going transmitted field in medium 2 to the incident field. For the 
case of two homogenous fluids — medium 2 parameters are equal to medium 3 
parameters, or d -»°o, then y = 0, and 

T(K) = 1 + T(K), (5-59) 


which agrees with Eq. (3.24). 
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5.5 First-order 

The first-order incoherent field scattered back into medium 1 and the 
downward traveling field scattered into medium 2 are derived in this section and 
compared with the first-order fluid-fluid results of chapter 3. 


5.5.1 Derivation of T-matrix 


For the first-order field, set N = \ in Eq. (5.51): 


'pj 1) (K") 
- A 


+ 


(1 + T(K")) 

KK" 

Pl(K") 

k \ 

(1 + T(K")) 

K K" 

Pl(K") 

. k l 


+ 




i f.2.1 (1 + r(K")) 


1 ~ Pi (K")Pi (K) 


(i-ran) 

(l + T(K")) 

(l - T(K")) 
(1 + T(K")) 




P p,(K")Kp 2 (K) 


2 K K" 

K -2— 


|F(K"-K) x Fj 0) (K) 

( 0 ), 


kP 2 (K)F(K" - K)[y +1]^2— (K) 


+-J d 2 K{ 1 - T(K"))kP 2 (K)F(K" - K)[y - 1]^?(K). 


(5.60) 


The first-order F-matrix then follows from Eq. (5.4), the zero-order quantities 
given in Eqs. (5.54) and (5.57), and from the above equation: 
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11(1+ran) \K■ K" R ™ R rin (1 - r (K"))"L F „ F , 

= 5 p l( K") "P - " 1 * 1 * 1 ® )Pl<K) (I7fan) F(K _K) 


1 ir(K)(i +ran) k-k" 

2 Pi OH k\ 


p!(K")Pi(K) 


(i+r(K"))J 

( ™ (i - ran) 


(5.61) 


(i+ran) 


ww - K > 


+ 11 ; ( 1 + r (K // )) k 2 _™ K p 2 (K)F(K"-K)(l + r(K)) 

2pp 1 (K")Kp 2 (K)[_ *? J V ” 

+|(1 - r(K"))lKp 2 (K)F(K" -K) ^jj(l + r(K)). 

By applying the second equality in Eq. (5.42), the above equation simplifies to 


(K", K) = H n (K", K)F(K" - K), 


(5.62) 


where 


H t 1 (K", K) = — l — {6(1 - r(K"))(l - r(K)) + a(l + r(K"))(l + r(K))}, (5.63) 

Z«Pl(lv ) 


, k 2 ( 1 . | K" • K 

a = 1-+ --1 - 

p Ip J kf 


(5.64) 


and 

b = p 1 (K")Pi(K)(p-1). (5.65) 

Reciprocity is manifest here. Also note that this result is identical to the fluid- 
fluid first-order result in chapter 3. In this result, the complexity of medium 3 is 
included in the reflection coefficient given in Eq. (5.43). Although it appears as 
though Eq. (5.63) follows from the two-homogenous fluid result in Eq. (3.37), one 
would not typically be able to generalize the two fluid result to Eq. (5.63) simply 
by substituting the general expression for the reflection coefficient given in Eq. 
(5.43) for the two-fluid reflection coefficient given in Eq. (3.23). Using Eq. (3.36) 
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in Eq. (3.35) to eliminate the icp 2 was necessary to obtain the expression 
identical to Eq. (5.63). If this substitution were not made, then the first-order 
results in chapter 3 would not generalize to the results in this chapter — simply 
by using the general reflection coefficient given in Eq. (5.43). 

The first-order field scattered through the surface is found by setting N = 1 in 
Eq. (5.21): 


^ 1 _ ) (K")(y" + 1) = ■ -J d 2 K i’Pj (K)F(K" - K)'Pj(K) 

+J d 2 K$ x (K)F(K" - K)^} 0) (K) 

-J AiKp 2 (K)F(K" - K)[y - 1]^(K) 

+4 /( f 1) (K"). 

7 (5.66) 

Using the zero-order results in Eqs. (5.54), (5.57), (5.58), and (5.62) in the above 
equation results in the following expression for the first-order field scattering 
through the rough interface: 


+ 1) = J A/Pi(K)F(K" - K)r(K) v F I -(K) 

-J d 2 K ipj (K)F(K" - KfViiK) 

+J d 2 Kif>fo (K)F(K" - K)(l - T(K))'F i (K) 

+ f d 2 KH x !(K",K)F(K" - 

J (5.67) 


The first-order F-matrix 7$ follows from the above equation and Eq. (5.5): 


fW(K",K)(y" +1) = /(p - l)Pi(K)(l - r(K))F(K" - K) 

+tf u (K",K)F(K"-K). 


(5.68) 


Similar to Eq. (5.62), define 

rW(K",K) S tf 12 (K",K)F(K" - K). 


(5.69) 
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If we use the representation for the transmission coefficient in Eq. (5.58), and the 
definitions of a, and b in Eqs. (5.64) and (5.65), H n follows from Eq. (5.68) and 
its definition in Eq. (5.69): 

tf, 2 (K".K) = + r < K » - *0 - r(K))}. (5.70) 

This result is a generalization of the result in Eq. (3.45), where the transmission 
coefficient relates the downward traveling field below the zero-mean rough 
surface to the incident field on the rough surface. For scattering back into 
medium 1, the complexity of medium 3 is contained in the reflection coefficient; 
for scattering through the sediment, the complexity of medium 3 is contained in 
both the transmission and reflection coefficients. 

5.5.2 The effect of shear waves on the differential cross section 

As a check of the formalism developed in this chapter, we reproduce a result 
from Essen (1994) that includes shear supporting sediment below a rough fluid- 
fluid interface. This particular example is an appropriate check, because the 
medium immediately below the rough surface is a fluid and Essen uses 
perturbation theory to calculate the differential cross section. 

The differential cross section obtained from first-order perturbation theory for 
scattering back into the water (medium 1) follows from Eq. (3.77) by replacing 
H \2 from Eq. (5.70) with H n from Eq. (5.63) and setting the speed ratio, v, to 

one. 

<T„(a / .a,) = ij 4 sin 2 (e 2 )|Hii(K / ,K ( -)| J w(K / -K,). (5.71) 

Here, 0 2 is the scattered grazing angle measured from the mean surface, and 
K r and K, are the transverse components of the scattered field and incident 
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field wave vectors, respectively. This example is concerned with the 
backscattering strength, lOlog^n), where Ky = -K,. 


In Eq. (5.63), the complexity of medium 3 is characterized by the flat surface 
reflection coefficient, T. Eq. (5.43) is used to find the reflection coefficient from 
the partial reflection coefficients, T n and r 23 . The partial reflection coefficient 

r 12 is the fluid-fluid reflection coefficient at the interface between medium 2 and 
medium 3 (z = 0). The partial reflection coefficient r 23 is the fluid-solid reflection 

coefficient at the interface between medium 2 and medium 3 at z = -d. When 
medium 3 is shear supporting, but gradient free, the partial reflection coefficient 
r 23 is given by (Brekhovskikh, 1980; Brekhovskikh and Lysanov, 1991) 



Zm ~ Z 2 

Z in + Z 2 ’ 


(5.72) 


where 

Z in = Z M cos 2 (2y 3 ) ■+ Z 3t sin 2 (2y 3 ), 


(5.73) 



(5.74) 



(5.75) 


and 



(5.76) 


The parameter K/ is the ratio of the compressional wave number of medium 3 to 
k h and the parameter k ( is the ratio of the shear wave number of medium 3 to 
kj. The complex angle y 3 is determined from 
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cos(y 3 ) = MK), (5-77) 

where 

P3U, = | - *V( W ■ < 5 ' 78 > 

In Fig. 5.2, the parameters given in Fig. 5c of Essen (1994) are used to plot the 

backscattered field strength, 101ogio(on(K/- K / = -K,, as a function of 

grazing angle. In this example, medium 2 consists of a very fine sand and has 
thickness d = X. The wave number and density ratios are taken to be 
k = 1/(1.12 - 0.004/) and p =1.85. The sedimentary rock in medium 3 supports 
both shear and compressional waves with wave number ratios 
k ( = l/(2.3 - 0.004i) and k, =1/(1.3-0.110- The ratio of the density in the 

sedimentary rock (medium 3) to that of the fluid in medium 2 is equal to 2.5. This 
example also uses the power-law roughness spectrum given in Eq. (3.55) or Eq. 
(3.66) with a = °°. Using a finite value for the filter parameter, a, would only 
affect the result at larger grazing angles than those plotted. Since y = 4.0, the 
backscattering strength does not depend on frequency. The spectral strength is 
given to be w 2 - 0.04/(2rc). This result is identical to the result obtained by 

Essen. 

5.5.3 The effect of sound-speed gradient on incoherent backscattered field 

In order to illustrate that gradients in sound speed can have a significant 
effect on the intensity scattered back into the water, the backscattered field 

strength, 101ogio(on(K/ K f =-K h is computed as a function of grazing 

angle for a sediment with a sound-speed gradient in medium 3, and is contrasted 
with the no sound-speed gradient case— all medium parameters in medium 3 
are identical to those of medium 2. In this example, we have followed Mourad 
and Jackson (1993) in assuming that the square of the wave number in the 
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medium below the rough seafloor decreases linearly with depth. The choice of 
linear decrease in squared wave number with depth, 


k%{Q = {\-aQkfK 2 , 


(5.79) 


where 

t, = -(z + d) 


(5.80) 


yields the well known solution for the field in medium 3 proportional to the Airy 
function, Ai(x), and solution for the reflection coefficient in terms of Airy functions 

(Brekhovskikh, 1980): 


ibhQQHAKtri-Aqta) 
231 j ik 2 $ 2 (K)HAi(t 0 ) + Ai'(t 0 ) 


(5.81) 


with 

ff = (* 2 2 a)' 1/3 , (5.82) 

t 0 = H 2 (K 2 -k}), (5.83) 

and 

x = H 2 a&% +H 2 (K 2 -kiy (5.84) 

The parameter a is expressed in terms of the sound-speed gradient, g, the 
sound-speed ratio, v, and the sound velocity in medium 1, q: 

a = 2g/{yc\), (5.85) 

where 



(5.86) 


The parameter g is taken to be the sound-speed gradient at the interface 
between medium 2 and medium 3. An example that shows the effect of a sound- 
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speed gradient on scattering back into the water is given below. In this example, 
medium 2 is taken to be arbitrarily thin; the reflection coefficient r can be found 
from substituting Eq. (5.81) into Eq. (5.43), and setting d = 0. 

PCIQ _ ikMK)pHAi(t 0 ) - Ai (fp) /g gy% 

ikM^)pHAi(t 0 ) + Ai'(t 0 )' 

Since the medium above z--d must be medium 2, z = -d must be at or below 
the lowest point on the rough interface. However the sound-speed gradient is 
small enough in this example that setting d = 0 results in a negligible difference 
in scattering strength. Since the backscattered field strength is plotted only to a 

grazing angle of 70°, values of the filter parameter a>lJm in the power-law 
roughness spectrum given in Eq. (3.66) result in the same scattered field 
strength. 

Figure (5.3) compares computations of seafloor backscattered field strength 
made with and without a sound-speed gradient for a frequency of 300Hz and a 

water sound speed of 1545 m -1 , with parameters appropriate to a soft silty or 

clayey bottom (p 21 = 1.548, v 2 i=056, and 8 = 0.1445xl0 -3 ) and an upward 

refracting sound-speed gradient at the interface of 1.3 s -1 . This example shows 
that the sound-speed gradient can cause an appreciable change in the seafloor 
scattering strength. The gradient in this example is high enough, and the loss is 
low enough, so that much of the energy that is transmitted through the interface 
is refracted upward to the interface to strike it again, resulting in constructive or 
destructive interference. The significance of this interference is discussed further 
in chapter 6. While this example illustrates the effect of large-scale gradients on 
low frequency scattering, very steep gradients on centimeter scales can result 
from biological reworking of sediments (Richardson, 1986) and may influence 
high frequency scattering. Gradients in porosity are equivalent to gradients in 
density (Hamilton, 1974), which can have strong effects on the reflection 
coefficient (Rutherford and Hawker, 1978; Holthusen and Vidmar, 1982). 
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Since the sound-speed profile at the experimental site is not known, the 
effect of gradients in sound speed on the incoherent and coherent intensity 
scattered through a rough interface is not included in the model used in the 
simulations (chapter 4). 

5.6 Second-order field 

For the ocean experiment simulation (Chotiros, 1995) we use a zero-order 
approximation in the coherent intensity calculation. There is a slight 
inconsistency in this approach, in that the incoherent intensity is proportional to 
the second power of k x h, while the coherent intensity is only computed to an 
accuracy of k x h to the zeroth power. The second-order field is found in the next 
section. Using this result, we obtain an expression for the second-order coherent 
field. 

5.6.1 Derivation of T-matrix 

The procedure for finding the second-order field is similar to finding the first- 
order field. However, both first-order and zero-order results are required in order 
to obtain the second-order field. Setting N=2 in Eq. (5.51), substituting K' for K 
in the F x terms, and using Eqs. (5.54), (5.57), (5.58), and (5.69) and 

F 2 ( K" - K) = J d 2 K'F( K" - K')F(K' - K) (5.88) 


from Eq. (5.9) yields 
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m>f\ K") = J d 2 K' J d 2 KF(K" - K')F(K' - K)^(K) 

..L , rfV ^ Pi(K) [k-k" r (l-r(K"» 

x { (1+r(K tell” W ^ K (i+r (K ")) 

lh , rrTr ,^ Pi(K) F k-k" r ¥ ( i - r(K")) 

+(1+r(K )} p^n "1T +1 Pl( )p i (K) (i + r(K")^ ( } 


(i+r(K")) 


(l-Y) Kp 2 (K)f 2 k-k 


-(i - r(K"))(i + r(K))(Kp 2 (K)) 2 } 

+2.ij d 2 K'\d 2 KF( K" - K')F(K' - K)^(K) 


x H u (K',K)(l + r(K")) 


K' K" 


(i - r(K")) 


^ Pi( k")L -a? 

+H n ( K',K)(l + r(K"))(l + Y')^ 
-H 12 (K',K)(1 - r(K"))Kp 2 (K')(l - Y')} • 




2 K' • K" 

K-*— 


(5.89) 


Define G n (K",K',K), such that 

4 / j- 2) (K") = J d 2 K'jd 2 KF(K" - K')F(K' - K)'P l (K)G 1 \ (K",K',K). (5.90) 

With much work and use of the second identity in Eq. (5.42), the above 
expression is equivalent to 
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Gn(K",K',K) 




(i+r(K"))(i - r(K»|^- - (1 - r(K"»(i+r®) 


“h 


iH n (K',K) 

Pi OH 


( 


K' K" 
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-1 
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: 1 


Pl(K") 


(1 +r(K")) - Pi (K")Pi (K')(i - r(K")) 


+/H 12 (K',K)(l + r(K")) 


) 

(1+y') 

Pi(K") 


x 


K' 


P P k 


i k' ■ k" 0 „ n (i - r(K")) i - r(K') 

—pP 1 (K")Pi(K') v ——— 


(i+ r(K")) i+ r(K') 


In order to find the second-order field scattering through the interface, set 
in Eq. (5.21). 

yW (K")(l + y") = - J d 2 K (pi (K)) 2 F 2 ( K" - K Wi (K) 

+J d 2 K 2ipj (K)F(K" - K)'FP(K) 

-J d 2 K F 2 (K" - Kf¥f\K) 

+J AkP 2 (K)F(K" - K)[i - yH-( k ) 

+J d 2 K (kP 2 (K)) 2 F 2 (K" - K)[l + yM? (K) 

+'fJ 2) (K"). 

The expressions 

¥^(K) = jd 2 K'H n (K, K')F(K - K')^ (K') 
and 

'pW(K) = J d 2 K'H l2 (K,K')F(K - K')^-(K') 


(5.91) 


N = 2 


(5.92) 


(5.93) 


(5.94) 


are obtained by combining Eqs. (5.4), (5.62) and Eqs. (5.5), (5.69). Combining 
Eqs. (5.57) and (5.58) yields 


145 


^(O)^) = 1+£(K) (K) 

(i+y) 

Using Eqs (5.88), (5.93), (5.94), and (5.95) in Eq. (5.92) yields 

ip( 2 ) ( K")(l + Y") 

= J d 2 K J d 2 K'F(K" - K')F(K' - K)^-(K) 

x{-(3? (K)(l + r(K)) + (kP 2 (K)) 2 (1 + r(K)) + 2iH\ j (K', K)^ (K') 

+2iff 12 (K',K)icP 2 (K')(l - Y0 + G n (K",K',K)}. 

Using the following definition for G 12 (K",K',K), 

4/(2) (K ") = J d 2 K$d 2 K'F(K" - K')F(K' - K)^(K')G 12 (K",K',K), 

in Eq. (5.96) yields 

G 12 (K",K',K)(l + Y") = {-pf (K)(l + r(K)) + (Kp 2 (K)) 2 (l + r(K)) 

+2iH{ j (K', K)p! (K') 

+2/H 12 (K / ,K)kP 2 (K / )(1 - YO + G„( K",K',K)} 
By noting the following identities, 

(p 1 (K)) 2 - (kp 2 (K)) 2 = k i^~ - = 1 - k 2 , 

Eq. (5.98) simplifies to 

(1 + y")Gi 2 (K",K',K) = (k 2 - l)(l + r(K)) + 2iH u {K',K)fo(K') 

+2iH x 2 (K', K)Kp 2 (K')(l - Y) + Gi l (K", K', K). 


(5.95) 

(5.96) 

(5.97) 

(5.98) 

(5.99) 

(5.100) 
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5.6.2 Second-order coherent field 

In order to check the zero-order approximation for the coherent intensity, the 
goal in this section is to find the coherent component of the second-order T- 

matrix, (r^ 2 ^(K",K)\. If we allow G to represent both G\\ in Eq. (5.91) and G \2 


in Eq. (5.98), it can be expressed as the sum: 

G(K", K, K ) = G a (K", K) + G b (K", K', K). (5.101) 

The 7-matrix is found from Eq. (5.90) or Eq. (5.97): 

T (2) (K",K) = J d 2 K'F{ K" - K')F{K' - K)G a (K",K) 

+| d 2 K'F(K" - K')F(K' - K)G fc (K",K',K). (5.102) 

Substituting 

f l (R) = f(R) (5-103) 

and 

/ 2 (R) = /(R)e iK ' R (5.104) 

into the identity 

—Ly f i 2 Rh (R)/ 2 (R)e- iK R = J d 2 K'Fi (K-K')i^(K') (5.105) 

(2k) z j j 

yields 

1 - T f d 2 Rf 2 { R)e _/(K,, “ K)R = f d 2 K'F{K"-K')F{K'-K). (5.106) 

(2k) 2 j j 


Note that / is a WSS process with unit variance. Therefore, taking the expected 
value of Eq. (5.106) over the ensemble of interfaces results in the following 
identity: 
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8(K" - K) = (J d 2 K’F x ( K" - K')F 2 (K'- K)}. (5.107) 

Consequently, the expected value of the first term on the right side of (5.81) 
results in the following expression for the coherent second-order F-matrix: 

/J d 2 K' F(K" - K')F(K' - K f}G a (K",K) = 8(K" - K )G a (K,K). (5.108) 

Define Int2 to be the second term on the right side of Eq. (102). 

Int2 = J d 2 K'F(K" - K')F(K' - K)G b {K",K',K). (5.109) 


Substituting 


F(K)= (2^^ 2/i:/(R)e ’ ,KR 


(5.110) 


for both F(K" - K') and F(K'-K) in Eq. (5.109), and taking the expected value, 
yields 


(Int2) = JA'G i (K /, ,K , ,K)|jj 2 Fe /(K ' K)R 

T , j2 p , „-«(K"-K 
(In)' 


x— ^ J d 2 R' e -'( K "- K, ) R/ c(R' - R )/h 2 ■, 


where 

C(R'-R)^ 2 (/(R')/(R)) 
is the surface autocorrelation function. By noting 

-Xr f d 2 R’ e-( K "- K ') *' C(R' - R) = e-'< K '- K '> R 
(2n) 2 1 

x ‘ C(U) 

(2tc) 2 J 

= e -«(K"-K')R w(K ,_ K0/ft 2 


(5.111) 


(5.112) 


(5.113) 


and 
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1 f d 2 R e - ! '( K '- K > R e-‘( K "- K ')- R w( K" - K') = 8(K - K")W(K" - K'), (5.114) 

( 2 n) 2 J 


Eq. (5.111) simplifies to 


(Int2) = 8(K - K")J d 2 K'G b { K",K\K) W(K" - K')/k 


(5.115) 


By using Eqs. (5.107) and (5.115), the expected value of Eq. (5.102) simplifies to 


= 5(K" - K){g o (K, K) + J d 2 K'G b (K, K', K) W{ K - K') A 2 } - (5-116) 

By making an approximation suggested by Thorsos and Jackson (private 
communication), Eq. (5.116) is approximately 

/r (2 >(K",K)} = 6(K" -K){G a (K,K) + G fc (K,K,K)J d 2 K' W( K - K')/ h 2 ] 

= 8(K" - K){G fl (K,K) + G fc (K,K,K)}. (5.117) 

To second-order, the expected value of the T-matrix is given by 

(r(K",K)> = ^7’ (0) (K",K)^ + ^(fc 1 /i) 2 ^r (2) (K // ,K)y (5.118) 


The intensity reflection or transmission coefficient, £ 7 , (Thorsos, 1990) follows 
from Eq. (5.118). 


S, = |S(K)| 2 s |S 0 (Kf + (hhf Re{^(K)^ 2 (K)}, 

(5.119) 

where 


^■(K) = ^ 0 (K) + |-(*lA) 2 *-2(K), 

(5.120) 

(7-(K",K)) = 8(K" - K)iE(K), 

(5.121) 

and (see also Berman, 1992) 
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/j( 2 )(K",K)^ = 5(K" -K)je 2 (K). (5.122) 

For the case of scattering from an interface, substitute the flat-surface reflection 
coefficient r(K) for je 0 (K) and *n(K) for ^( K )i for the case of scattering 
through an interface, substitute the flat surface transmission coefficient 7(K) for 
*o(K), and *i 2 (K) for *(K). 

For scattering from the interface, Eq. (5.91) along with the definition in Eq. 
(5.101) yields 


GnaCK'.K) 




(i+r(K"»(i - - (i - r(K"»(i+r<K» 

Pl(A 1 


(5.123) 


and 


G llfc (K",K',K) 


iH n (K',K) 

Pl(K") 


x 


r 


K' K" 


\ 


-1 


V 


■l 


(i+r(K")) - p^K'OPitK'Xi - r(K")) 


r K 2 lK'-K'^ 


V 


P P k 


l J 


r 


(1 + nm) - pp,(K")Pi(K')(i - r(K")) 


v 


i - hkq 

i+ r(K') 


\ 


) 


(5.124) 


Likewise, for scattering through the interface, using Eq. (5.100) and the definition 
in Eq. (5.101) yields 


G 12 o (K",K) = -J-(k 2 - l)(l + r(K)) + -J-Gn^K'-.K) 


(5.125) 


and 
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G m ( K",K',K)(l + y") 

= 2iH\ j (K', K)P! (K') + 2iH} 2 (K', K) K p 2 (K')(l - y') 
+G nb { K",K',K). 


5.6.3 Second-order coherent field calculations for incident plane wave 

In Eq. (5.119), £7 consists of the zero-order term and the second-order term; 
the first-order T-matrix is zero mean. In order to compare the second-order term 
to that of the zero-order term, the second term on the right side of Eq. (5.119) 
divided by the first term on the right side the equation is plotted in Fig. 5.4 as a 

function of the RMS (root mean square) surface height. The grazing angle is 20 ° 

in Fig. 5.4a, and it is 40° in Fig. 5.4b. The solid curve results from using Eq. 
(5.116), the dashed curve is due to Eq. (5.117). 

5.6.4 Modification of coherent field calculation 

Taking the expected value of both sides of Eq. (5.5) and using Eq. (5.121) yields 
('P 2 .(K)) = ^ 12 (K)>P,(K). (5.127) 

For the case in which medium 3 parameters are identical to medium 2 
parameters, 

\j/ 2 _(r) = |A 4 / 2 _(K)e"' 7 :iPl(K) V' K ' R . (5.128) 

Taking the expected value of the above equation and using Eq. (5.127) results in 
an expression for the coherent field penetrating the rough interface. 

( ¥ 2 .(r)> = j4 2 Ar? 12 (K)>P / (K)e- , '*' | 5 > <K)z e iK R . (5.129) 

Likewise, the coherent field scattered from the interface is given by 

( ¥/ (r)) = JAs„(K)>P,(K)e it ' |i > (K, V KR . (5.130) 
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Note that % 12 (K) in Eq. (5.129) is a coherent transmission coefficient, but ^(K) 
in Eq. (5.130) is a coherent reflection coefficient. Referring to Eq. (2.91), the 
exact expression for the zero-order field is easily modified by substituting j£ 12 (K) 
for the transmission coefficient. 

<V 2 _(r)) = ,j<iA:4^l^le i(z <‘ 1? ' (K) ' z * 2|i2(K)) J 0 (KK). (5.131) 


The square of the time-dependent coherent field magnitude scattered through 
the rough interface due to an incident point source with baseband source signal 
U( cd) is therefore 


h cM = 


oo 


J^cot/(co) 


-ICO? 


f r 0^12( K ) c >tz.^lPi(K)-^ 2 P2(K)) J f KR \ 

{ *lPl(K) 


(5.132) 


5.7 Summary 

Nth-order perturbation recursion equations are derived in this chapter that 
allow for simple computation of the first and second-order fields for scattering 
through and from a fluid-fluid rough interface. An expression for the second- 
order coherent field is derived that illustrates when it is necessary to include the 
second-order field in the coherent intensity calculations for scattering through 
rough interfaces. In the ocean experiment simulations, k\h=l, and including 
the second-order component reduces the coherent field calculated with the zero- 
order result by almost ten percent. Since the peak incoherent intensity and 
coherent intensity typically differ by over an order of magnitude for all incident 
angles, and at all receivers, the difference in the calculated coherent field by 
excluding the second-order correction does not affect the intensity velocity 
algorithm output. 
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The results in this chapter are general in that the medium below the lowest 
point on the rough interface is allowed to be vertically stratified, viscoelastic, or 
porous, and is contained in the reflection coefficient of the corresponding flat 
surface problem. For the case of increasing sound speed with depth, these 
results show that an interference effect can occur when enough energy is 
redirected back toward the interface. Since the sediment parameters such as 
sound speed, sound-speed gradient, and density are not everywhere constant, 
the effect of treating the sediment parameters as random is discussed in the next 

chapter. 
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V|/ i : incident field 


lj if : scattered field 




Medium 1 




Figure 5.1 : Definition of the scattering problem for an arbitrary two-dimensional rough 
interface. Medium 1 (z > hf( R)) and medium 2 (-d <z< hf^R.)) are homogenous 

fluids; the medium parameters can be a function of z for (z < -d ). This diagram can be 
viewed as a slice through a two-dimensional interface. 
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grazing angle (degrees) 


Figure 5.2 : Medium 2 is very fine sand with d = X, k = 1 /(1.12 - i • 0.004), and p = 1 .85. 
Medium 3 is sedimentary rock with k ( = 1 /(2.3 — / - 0.004), K t = 1/(1.3 — * - 0.11), and 
p 3 /pj =2.5. Surface: w 2 = 0.04 /(2jc), y = 4.0. The solid curve uses these 
parameters. For the dashed curve, k,. = 1/1.13. For the dotted curve, there is no 
layering, i.e., d —» oo, 
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Upward refracting sound-speed gradient 
0 _- 


CD 



XJ 


-100 f ■ | i-1 ■ | < | ' | ' 

0 10 20 30 40 50 60 

grazing angle (degrees) 


Figure 5.3 : Bottom backscattering strength with and without a sound-speed gradient 
for soft, silty, or clayey bottom: y = 3.25, w 2 = 1.636 x 10 -4 m 3 / 4 , frequency = 300 kHz, 
p = 1.548, 5 = 1.445 x 10 -4 , v=0.96, q =1545m/s, gradient = 1.3 s" 1 . 
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Coherent intensity ratio 





Figure 5.4a Ratio of second-order coherent term to zero-order term, 
(k\h) 2 Rej^o(K)*?2(K) J /\%o (K)| 2 with frequency = 20 kHz, p = 2.0, 5 = 0.0163, 

y — 3,0, w 2 =6.2xl0~ s m, v=1.13, q = 1536 m/s, and incident grazing angle, 

e, = 20 °. 
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co 




Figure 5.4b Ratio of second-order coherent term to zero-order term: 
{k x hf Re{^o(K)^ 2 (K)}/|^o(K)| 2 with 0 y =40°, otherwise same parameters as Fig. 

5.4a (frequency = 20 kHz, p = 2.0, 8 = 0.0163, y = 3.0, m >2 = 6.2x10 v=1.13, 

C\ -1536 m/s). 







CHAPTER 6 


First-order perturbation solution for scattering from a rough seafloor 
including the effects of random subsurface stratification 


6.1 Introduction to chapter 6 

The example in Fig. 5.3 shows that upward refraction due to a sound-speed 
gradient in the sediment results in constructive and destructive interference in the 
field scattered back into the water. This example assumes the sound-speed 
gradient and other seafloor parameters are not random. Since seafloor 
parameters such as sound-speed gradient and density actually vary with 
position, a more suitable model may be one that takes into account the inherent 
randomness of the seafloor parameters. In the last chapter, a first-order 
expression for the scattering cross section per unit area per unit solid angle 
(differential cross section) is derived that includes the complexity of the medium 
below the lowest point on the rough interface (medium 3) in the flat surface 
reflection coefficient. In this chapter, the reflection coefficient, r 2 3 , is modeled as 
a random variable. The expected value of the differential cross section over the 
random sediment parameters is then expressed analytically in terms of reflection 
coefficient moments. Analytical expressions for the reflection coefficient 
moments are found using ray theory. These analytical results are compared with 
modeling the sound-speed gradient as a random variable with Gaussian 
probability distribution (Moe and Jackson, 1994b), where the expected value of 
the differential cross section is calculated by means of a numerical integration. 
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6.2 Averaged scattering strength with random sound-speed gradient 


The first-order perturbation solution for the differential cross section for 
scattering from the rough interface can be found by combining Eqs. (5.71) and 
(5.63). 



where 



and 

^(K/’K/) = Pi(K/ )P 1 (K,)(p -1) 


( 6 . 2 ) 


(6.3) 


are given by Eqs. (5.64) and (5.65). The vectors K f and K, are the transverse 

components of the scattered field and incident field wave vectors, respectively. 
The mean (flat) surface reflection coefficient, T, is written explicitly in terms of the 
reflection coefficients r 12 and r 23 in Eq. (5.43), where r 12 represents the mean 

surface reflection coefficient at z = 0 as c/ and r 23 is the reflection 

coefficient corresponding to the interface between medium 2 and medium 3 
(z = -d). All of the complexity of the lower medium is contained in r 23 . As 

shown in Fig. 5.3, the lower medium, characterized in Eq. (5.43) by r 23- can 
have a notable effect on the scattered field strength if it redirects a significant 
amount of the downward propagating energy in the sediment back up towards 
the rough interface. Two first-order effects as depicted in Fig. 6.1 can result in 
scattering energy away from the interface, resulting in destructive or constructive 
interference. 
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Since sediment parameters are typically not constant, the interference 
effects in Fig. 5.3 may not actually be observable. By treating the gradient as a 
random variable, the expression for the differential cross section given in Eq. 
(6.1), weighted by a Gaussian probability distribution for the sound-speed 
gradient (linear increase in squared wave number with depth) with variance 

(\g ~(#)| 2 ^ = 0.1 sec -2 and mean (g) = 1.3sec _1 , is numerically integrated. The 

resulting averaged backscattering strength is plotted in Fig. 6.2 as a function of 
grazing angle and compared with the backscattered field strength due only to the 
rough seafloor surface with no gradient, r 23 = 0. In this example, a sound-speed 

gradient in the sediment results in a substantial increase in the average 
backscattered field strength. This scattering gain occurs at all grazing angles, 
but is greater at smaller grazing angles. Large grazing angles result in long 
refracted paths in the lossy sediment, and therefore contribute less to this gain 
than small grazing angles with shorter paths. 

Using the same parameters as in the above example, we plot the bistatic 
scattering strength in Fig. 6.3 as a function of scattered grazing angle for an 

incident grazing angle 9, =15°. The scattered direction in all plots is chosen 
such that the incident transverse wave vector, K,-, and the scattered transverse 
wave vector, K f , are in the same line, and the scattered grazing angle, Qf, is 

measured from the horizontal in the backward direction. In this case, 

backscattering (K^ =-K,) corresponds to an angle of =15°. Scattering gain 

due to the upwardly refracting sediment occurs in all scattered directions, 
including the backscattered direction, where a local scattering gain maximum, or 
backscattering enhancement, is apparent. The reason for this backscattered 
peak, as illustrated in Fig. 6.4, is that the coherent sum of the two paths 
corresponding to the incident and scattered direction occurs only in the 
backscattered direction, because only in this direction are the paths ABCD and 
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D'CBA' of equal length. This effect is illustrated for similar scattering problems 
by Lang and Khadr (1992) in their analysis of remote sensing of a vegetation 
canopy over a flat ground, and in Hanson and Zavorotny (1995) in their analysis 
of backscattering from an ocean-like surface. The width of the peak depends on 

the sound-speed gradient variance, ^|g-(g)| 2 y The larger the sound-speed 

gradient variance, the narrower this peak. Another oscillation that depends on 
the sound-speed gradient can be observed just below a scattered grazing angle 
of 10°. At these angles, the oscillations in the scattering strength are not quite as 
rapid as those at larger grazing angles, and the variance chosen for the gradient 
is not sufficient to completely wash out the dependency of scattering strength on 
scattered grazing angle. 

6.3 Formally averaged differential cross section with random sediment 
parameters 

Rather than using numerical integration as in the last example, using a 
general approach to this problem leads to a formally averaged analytic 
expression for the differential cross section per unit area that includes the 
randomness of seafloor parameters below /(R), in addition to the randomness of 
the interface. All of the complexity of the lower medium in Eq. (6.1) is contained 
in the reflection coefficient, T. Consequently, the inherent randomness of the 
sediment is equivalently represented in the random reflection coefficient, r. This 
point of view implies that the expression in Eq. (6.1) is actually a conditional 
expected value over the rough interface given deterministic sediment parameters 
below the interface, or a deterministic reflection coefficient, and the goal is to find 
the expected value of Eq. (6.1) over the random fluctuations in the seafloor 
parameters below z = /(R), 
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(6.4) 


The approach is to express Eq. (6.4) in terms of moments of the flat surface 
reflection coefficients. The notation ry = T(K y) and T, = r(K,) is used for the 

reflection coefficient in Eq. (6.1). Expanding the square in Eq. (6.4) yields 






+ 2frRe( a (i+1>)(1+r,)(i - r;)(i - r*)) 


(6.5) 


For monostatic scattering, the random variables Ty and T, are equal 
(K f = -K/), and Eq. (6.5) becomes 



+^Re(fl)[2-8(|r| 2 ) + 2(|T| 4 ) 

+ 8 i>Im(«)[lm(r|r| 2 ) - im(r)]J. 


+ 4Re(T 



( 6 . 6 ) 


Small changes in medium parameters result in rapid changes in phase of the 
reflection coefficient. The reflection coefficients, iy and ry are therefore 

assumed to be independent random variables for bistatic scattering. However, 
for scattering in directions close to the backscattered direction, the reflection 
coefficients corresponding to the incident and scattered directions are neither 
equal or independent (see Fig. 6.3). For independent reflection coefficients, Eq. 
(6.5) becomes 
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k * 

d n =^W(K f -Ki) 

x||aj 2 i + 2Re(r / )+^|r / | 2 ^ l+2Re(r / )+^|r i | 2 ^ 

+b 2 l-2Re(r / ) + /|r / | 2 \ l-2Re(r / ) + (|r i | 2 ^ 

+2bRc(a) f 1 - - (|T f f}) - 4Im(r / )lm(r i ) 

-4Z>im(a) im(r^i - (|r / | 2 }] + i m (r / )(i - (|r*| 2 )) >• 

In order to find expressions for the moments (r), (r 2 ^, (|r| 2 ^, (r|r| 2 J, and (|r| 4 ^, 

an equivalent form of Eq. (5.43) is used to express the reflection coefficient, r, in 
terms of the partial reflection coefficients r 12 , and r 2 3 as a geometric sum, 


r 2 3 e 2ik2M X 

n =C 

Small variations in medium parameters result in small deviations of the 
magnitude of the reflection coefficient r 23 from its mean, but result in large 
fluctuations of its phase. The partial reflection coefficient r 23 is assumed to have 
the form 


r = r 12 + (i-r 2 2 ) 


{-ri 2 r 23 e 2012 ^) . 


( 6 . 8 ) 


r 23 (K) e 2ik2 ^ 2 ^ =ti(K)e^, 


(6.9) 


where 0 is a uniformly distributed random variable between 0 and 2ti, and 
r| 2 (K) is approximated as the mean square reflection coefficient /|r 23 (K)| 2 y 
With much work, Eqs. (6.8) and (6.9) yield 


<r> = r 12 , 


( 6 . 10 ) 
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( 6 . 12 ) 


r|rf) = r 12 r 12 + 2 r 12 i-r/ 2 
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-rr 2 (i-ri 2 2 ) 
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T) 


(6.13) 


(6.14) 


The derivations for Eqs. (6.12), (6.13), and (6.14) are given in Appendix E; Eqs. 
(6.10) and (6.11) follow immediately from Eqs. (6.8) and (6.9). 


6.4 Effect of sound-speed gradient 


The above formalism is now applied to the case of a linear decrease in 
squared wave number with depth introduced earlier (Figs. 5.3, 6.2, and 6.3). In 
order to use the differential cross section expressions in Eqs. (6.6) and (6.7), the 


value for r| 2 (K) = /|r 23 (K)| 2 ^ (d 2 = 0) is computed and used in Eqs. (6.10) 

through (6.14). These equations are then substituted in Eq. (6.6) to obtain the 
monostatic differential cross section, or in Eq. (6.7) to obtain the bistatic 
differential cross section for scattered directions not close to the backscattered 


and specular directions. Rather than use Eq. (5.81) to determine t| 2 (K) exactly, 
a computationally efficient, yet accurate, ray theory expression (Appendix F) is 

used to obtain the following analytical expression for ^|r 23 (K)' 2 
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T| 2 (K) = exp 


| sin 3 (0 2 ) + 2 cos 2 (0 2 ) sin(0 2 ) 

\ 


l <L 

_3 

J 


(6.15) 


where 0 2 is the grazing angle just below the mean surface, 


cos(0 2 ) = 


K 

Re(*2)‘ 


(6.16) 


The gain parameter, q, is the sound-speed gradient, g = dc/dz evaluated at 
z = 0, divided by the product of the frequency, /, and the loss parameter, 6: 


g _ 1 dX 
fb 8 dz 


(6.17) 


This parameter will be shown to determine the average gain of the scattered field 
strength due to gradients and discontinuities below the rough water-sediment 
interface as compared to the scattered field strength due solely to the rough 
water-sediment interface. For example, a small value of q due to a relatively 
large loss at a high frequency will result in the attenuation of much of the signal 
energy before it strikes the rough interface again from below, resulting in little 
scattering gain over the homogenous fluid-fluid case. Also, a small sound-speed 
gradient results in large path lengths, and attenuation of much of the signal 
energy before it strikes the rough interface for the second time. When q is small, 
only the sediment properties at the water-sediment interface are of importance in 
the scattered field strength calculation . 

An expression for r| 2 (K) for the linear increase in sound speed with depth 
case is also derived in Appendix F. 
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6.4.1 Effect of sound-speed gradient on backscattered field strength 

As a check of the analytical results, the backscattering strength using Eq. 
(6.6) is plotted in Fig. 6.5 along with the numerically obtained result from Fig. 6.2, 
and the bistatic scattering strength obtained from Eq. (6.7) is plotted in Fig. 6.6. 
It is evident that the scattered field strength obtained by numerically averaging 
the expression in Eq. (6.1) over a Gaussian probability distribution function for 
the gradient is indistinguishable from the result due to the formally averaged 
differential cross section in the backscattered direction and for bistatic scattering 
at least a few degrees from the backscattered direction. In this example, the 
gradient gain parameter has the relatively large value # = 30. Equation (6.6) is 
used to find the point on the plot corresponding to backscattering. 

When Kf =-K ; , then r(K / ) = r(K l ), resulting in a peak in the 

backscattered direction. Depending on the value of the sound-speed gradient 
variance in the numerical calculation, the reflection coefficients T(Ky) and r(K,) 

are correlated when K f is close to K it and therefore the bistatic scattered field 

strength is dependent on the gradient variance in the region close to the 
backscattered direction. However, even for a small variances in seafloor 
parameters, this region is very small. 

The moments obtained in Eqs. (6.10) through (6.14) using the representation 
for the random reflection coefficient in Eq. (6.9) also agree very closely with the 
corresponding reflection coefficient moments obtained numerically with a 
Gaussian probability distribution function for the sound-speed gradient. 

Figure 6.7a illustrates the dependency of the scattering gain, 

101ogio(dn/cii), (Ki = -Kf), where Gji is given by Eq. (5.71), on the mean 

density ratio and mean sound-speed ratio for q = 30. This value for q is the 
same as that of Figs. 6.2, 6.3, 6.5, and 6.6. Fig. 6.7b.is an identical plot, but with 
a lower gain parameter, q= 10, resulting in correspondingly lower values for the 
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scattering gain. For a density ratio close unity, the scattering gain can be very 
large. However, the density ratio must always be greater than one for the 
seafloor. 

The dependence of the scattering gain on q and grazing angle is shown in 
Fig. 6.7c. As can be seen in Fig. 6.3 and Fig. 6.7c, the scattering gain is strongly 
dependent on the grazing angle. A low grazing angle results in a significantly 
higher scattering gain because the resulting ray path length in the lower sediment 
is much shorter and therefore attenuated much less than the large grazing angle 
case. The dependency of the scattering gain on q is greater for larger grazing 

angles. 

6.4.2 Bistatic scattering and backscattering enhancement 

Figures 6.8a, 6.8b, 6.8c, and 6.8d use Eqs. (6.6) and (6.7) to show 
backscattering enhancement corresponding to the same sediment parameters in 
Figs. 6.7a, 6.7b, 6.7c, and 6.7d. Note that these results do not depend on 
interface roughness. Since the two first-order paths through the sediment add 
coherently only for monostatic scattering (see Fig. 6.4), the field strength 
scattered back in the direction of the transmitter is greater than the field strength 
scattered in other nearby directions. As is clear from Fig 6.8c, the dependency of 
backscattering enhancement on q is negligible. 


6.5 Conclusion 

A formally averaged expression for the scattered field that includes the 
random stratification below the rough seafloor surface as well as the randomness 
of the rough seafloor is derived in this chapter. Although this formalism can be 
applied to a rough interface above any arbitrary vertically stratified sediment, ray 
theory is put to use to obtain formally averaged closed form analytical 
expressions for the special case of linear decrease of squared wave number with 
depth sound-speed gradient. This result shows that the scattering gain due to 
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upward refracting sound-speed gradients is determined by the quantity 
q = g/(fo), the sediment parameters (speed ratio, and density ratio), and also the 

incident angle. 

Gradients and layering below the seafloor can strongly affect the scattering 
from the interface if a significant amount of the energy that penetrates the 
interface is redirected back up towards the interface. The effect of gradients in 
sound-speed on scattering can be significant even for incident energy at large 
grazing angles (for situations with large q — large gradient, low frequency, low 
loss), but the effect of gradients on the scattered field strength is most significant 
at small grazing angles. Backscattering enhancement is also shown to occur in 
sediments that exhibit a relatively large scattering gain. By using numerical 
solutions for the flat surface reflection coefficient, the monostatic and bistatic 
differential cross sections of a rough interface above an arbitrarily complicated 
vertically stratified sediment that includes the randomness of the sediment 
parameters below the interface, in addition to the randomness of the interface, 
can be solved using the formalism in this chapter. 

The results of this chapter can also be used to obtain expressions for the 
field scattered through the rough interface when the sound speed is depth 
dependent. These results could then be incorporated in the model used in the 
ocean experiment simulation discussed in chapter 4. Except for downward 
propagating energy, sound-speed gradients would increase the path length of 
energy scattering into the sediment. Consequently, sound-speed gradients 
would seem to favor downward propagating energy, resulting in a slower 
apparent slow wave propagation than that predicted by the model used in 
chapter 4. The effect of sound-speed gradients on the scattered field through the 
interface is, however, beyond the scope of this dissertation, and a subject for 
future investigation. 
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Constructive or destructive interference 


Increasing 
sound spee 




Figure 6.1 : Two first-order effects that can result in destructive or constructive 


interference 
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Average effect of gradients 



grazing angle (degrees) 


Figure 6.2 : Averaged backscattering strength. Linear decrease in squared wave 
number with depth sound-speed gradient, g, with Gaussian probability function 


= i 3 /|g| 2 - ^ = 0.1 s' 2 ). All other parameters same as Fig. 5.3. 





bistatic scattering strength (dB) 
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Numerically averaged bistatic scattered field strength 
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incident angle = 15 degrees 
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Figure 6.3 : Numerically averaged bistatic scattering strength; same parameters as Fig. 
6.2 and Fig. 5.3. 
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Backscattering enhancement 





Figure 6.4 : Backscattering enhancement. If paths ABCD and D'CBA' are equal 
(A = A' and D = D'), then backscattering enhancement due to the sound-speed 

gradient occurs. 
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Comparison of analytical and numerical results 



Figure 6.5 : Average backscattering strength: comparison of numerical average 
obtained with random sound-speed gradient (Fig. 6.2) with analytical formal average 
approximation. Same parameters as Fig. 6.3. 
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Comparison of analytical and numerical results 



scattered grazing angle (degrees) 


Figure 6.6 : Average bistatic scattering strength: comparison of numerical average 
obtained with random sound-speed gradient (Fig. 6.3) with analytical formal average 
approximation. Same parameters as Fig. 6.3. 
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Scattering Gain (dB), q = 30 



speed ratio 


Figure 6.7a : Dependency of scattering gain, 101og 10 (d/a), (k,- = -Ky), on the mean 

density ratio and mean sound speed ratio for q = 30 (frequency = 1 kHz, g = 30 s \ 
8 = 0.001, incident angle = 45°) 
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Figure 6.7b : Same parameters as Fig. 6.7a except g - 10 s \ resulting in q - 10 
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Backscattering Enhancement (dB), q = 30 
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Figure 6.8a : Backscattering enhancement as a function of p and v using same 
parameters as Fig. 6.7a (q = 30). 
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Figure 6.8b : Ba 














CHAPTER 7 


Dissertation Summary and Proposed Future Work 


7.1 Introduction to chapter 7 

By means of analytical models for seafloor scattering derived in this 
dissertation, roughness is shown to be an important mechanism for penetration 
of acoustical energy into the sediment at subcritical grazing angles. 
Experimental results (Boyle and Chotiros, 1992; Chotiros, 1995) showing 
acoustic penetration from water into sandy sediments at grazing angles below 
the compressional critical angle in relation to the mean surface are explained in 
terms of diffraction of an ordinary longitudinal wave, rather than refraction of a 
slow wave. 

Since the transmitted signal in the experiments is a pulse, formalism for 
scattering of a pulse through a rough interface as well as from a rough interface 
is derived in terms of the scattering cross section per unit area per unit solid 
angle (differential cross section). This analytical result for the time-dependent 
scattered field intensity is used to show that the CW approximation is appropriate 
for simulations of both experiments. First-order perturbation theory is used to 
obtain an approximate expression for the differential cross section. Since using 
the zero-order field to represent the coherent intensity is only valid for small 
amounts of interface roughness, the second-order component of the coherent 
field intensity is derived, and compared with the zero-order component of the 
coherent intensity for the sediment roughness parameters used in the ocean 
experiment simulation. These analytical results justify using the flat surface field 
to represent the coherent component of the field. Although the main problem is 
acoustic scattering through the water-sediment interface for an incident field in 
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the water, analytical results derived in this work also include scattering from the 
interface. Furthermore, these analytical results are general, and can be applied 
to other rough interface scattering problems. 

Gradients in sediment sound speed are also allowed in the perturbation 
results. A novel approach for finding the differential cross section, given 
vertically stratified random sediment parameters, yields analytical results that 
show the average effect of sound speed gradients in the sediment on the field 
scattered back into the water. Although gradients in sound speed may affect the 
apparent slow wave speed, this effect is not calculated in this dissertation. 
However, the formalism in this dissertation can be generalized to include the 
effect of sound-speed gradients on the field penetrating the surface. 

7.2 Time dependent scattered intensity — chapter 2 

An original expression for the two-frequency mutual coherence function for 
scattering from and through a 2-D fluid-fluid rough interface due to a narrowband 
incident plane wave is derived and expressed in terms of the second moment of 
the r-matrix, resulting in a general analytical expression for the time dependent 
incoherent intensity due to a narrowband incident plane wave. 

Two approximations are made to the two-frequency mutual coherence 
function when obtaining the results of chapter 2. 

• The T-matrix correlation function is assumed to vary slowly over the source 
frequency range. This assumption is shown to be valid in chapter 3 for a 
narrowband source signal. 

• The phase terms are expanded in a power series in frequency. 

An equivalent expression, omitting evanescent waves, is obtained in terms of the 
differential cross section. It has been demonstrated that the evanescent 
component of the scattered incoherent field is insignificant compared to the 
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propagating component of the incoherent field for distances as small as one or 
two wavelengths from the interface, depending on the surface roughness 
spectrum. In the limit as the narrowband signal approaches the CW (continuous 
wave) case, the result in terms of the T-matrix correlation function is exact, and 
the result in terms of the differential cross section approaches the exact solution 
for distances sufficiently far from the interface. This result is interesting, because 
the differential cross section is defined as a far-field entity. 

This expression for the incoherent time dependent scattered field intensity in 
terms of the differential cross section is used to heuristically obtain the 
corresponding expression for the scattered field due to a narrowband point 
source. A simple analytical expression is obtained for the case of a Gaussian 

shaped source pulse, also in terms of the differential cross section. The 

* 

analytical results and calculations in this chapter show that the sonar equation is 
appropriate for the simulations of the experiments. 

7.3 First-order perturbation theory — chapter 3 

Derivations for the first-order perturbation expressions used for the numerical 
computations in chapter 2 and in the computer simulations (chapter 4) for 
scattering both from and through a rough fluid-fluid interface separating two 
homogenous fluids are presented in chapter 3. These first-order perturbation 
theory results are also used to show that the T-matrix correlation function varies 
sufficiently slowly over the source frequency range so that its dependence on 
frequency can be ignored in the time dependent scattered field derivation. 

7.4 Computer simulations of recent experiments — chapter 4 

The results of computer simulations of recent experiments, presented in 
chapter 4, show the effect of roughness on acoustic penetration of the seafloor 
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and are compared with recent sediment penetration experiments. These results 
show that the incoherent intensity scattered through the rough seafloor can 
mimic a slow refracted wave. 

By assuming a small amount of roughness at the water-sediment interface, 
and using the 3-D scattering model developed in chapter 2 and the first-order 
perturbation theory results for 2-D fluid-fluid interfaces derived in chapter 3, it is 
possible to reproduce the results of the tank experiments of Boyle and Chotiros 
(1992), both in magnitude and arrival time of the incoherent intensity pulse 
(apparent slow wave). These results indicate that the acoustic penetration of the 
surface may be due to scattering (diffraction) from low levels of roughness rather 
than slow-wave refraction. 

Simulations of the ocean sediment penetration experiments (Chotiros, 1995) 
are presented and compared with Chotiros' experimental results. Assuming 
roughness parameters appropriate to the experiment location and using the 
intensity propagation speed and direction finding algorithm used by Chotiros, 
these simulations also show propagation of energy at speeds slower than the 
compressional speed of the sediment, in a direction concurring with Snell s law. 
Although the simulations yield an apparent slow wave that is between 50m/s and 
300m/s faster than the results presented by Chotiros, the apparent slow wave 
speed is dependent on the roughness spectrum. Because of the potential for 
disturbance of the sediment above the receivers, the roughness spectrum 
measured in the area of the experiment is not the same as that of the actual 
roughness spectrum at the experiment site. 

7.5 Perturbation theory results for a fluid-fluid interface including 
gradients — chapter 5 

Using Rayleigh-Rice perturbation theory, Nth order recursion equations are 
derived in chapter 5 that provide a starting point for evaluating the Nth order 
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scattered field in terms of lower orders. The sediment below the lowest point on 
the interface can be viscoelastic, or porous, supporting shear or Biot slow waves 
as well as gradients in physical properties (e.g., sound speed can vary with 
depth), and is represented by the corresponding flat surface reflection coefficient. 
For the case of increasing sound speed with depth, these results (first-order) 
show that constructive and destructive interference in the field strength can occur 
for scattering back into the water when a significant amount of energy is 
redirected back toward the surface. 

An expression for the second-order coherent field is derived from the Nth 
order recursion equations that illustrates the necessity of including the second- 
order component in the coherent intensity calculations for k\h >1. 

7.6 Formal average of sediment parameters — chapter 6 

Since the sediment parameters such as sound speed, sound speed gradient, 
and density generally show variability in lateral direction, the effect of treating the 
sediment parameters as random is shown in chapter 6. The first-order 
expression for the differential cross section derived in chapters 2 and 5 is 
generalized to account for randomness in the medium below the lowest point on 
the rough interface by expressing the differential cross section in terms of 
moments of the stochastic flat surface reflection coefficient. This approach 
results in a formally averaged expression for the scattered field that includes the 
random stratification in the sediment, in addition to the randomness of the rough 
water-sediment interface. 

Although this formalism can be applied to a rough surface above any 
arbitrary vertically stratified sediment, ray theory is used to obtain formally 
averaged closed form analytical expressions for a few special cases that include 
sound speed gradients. This result shows that the scattering gain due to upward 
refracting sound speed gradients is determined by the quantity q = g/{f 5), the 
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sediment parameters (speed ratio and density ratio), and also the incident angle. 
Upward refraction or reflection due to large sound speed gradients in the 
sediment results in increased scattering strength over the non-gradient case. 
Backscattering enhancement is demonstrated in this first-order result, since the 
scattering gain has a relative maximum when the incident and scattered grazing 
angles are equal. 

Summary and ideas for future work 

The tank and ocean experiment simulation results presented in this 
experiment show that surface roughness is a mechanism for subcritical 
penetration of sediment. Although the issue of the dominant mechanism for 
subcritical penetration of the seafloor can only be resolved by experiment, further 
improvements on the scattering model can be made. For example, the effect of 
sound speed gradients on the scattered field through the interface is a subject for 
future investigation. The results of chapters 5 and 6 can be used to obtain 
analytical expressions for the field scattered through the rough interface when 
the sound speed is depth dependent. These results could then be incorporated 
in the model used in the ocean experiment simulation. 

The Nth order recursion equations provide a starting point to solve higher 
orders. Solving for the third-order scattered field would provide a means of 

obtaining the differential cross section to (kfi) 4 . This result could then be used to 
determine the accuracy of the first-order results for scattering through and from 
2-D surfaces. 

The results in chapter 5 provide a straightforward approach for obtaining the 
first-order scattered field from multiple rough interfaces in terms of partial 
reflection coefficients. The effect of multiple rough interfaces could then be 
compared with the results of chapter 6. 
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APPENDIX A 


Derivation of narrowband dispersion coefficients 


Equations (2.38) - (2.41) (expressions for t h t 2 , P, and Q) are derived in this 
section. For convenience, define 

co; = co c + G)' (A1) 

©"=© c + ©". 

From Eqs. (2.10), (2.36), and (A1), 

s^,©;,©;^ (©; M )kP 2 (1k-Kj/2|,©;) 


= K 


-K 2 + K-K d --K d -K d 


Substituting Eq. (2.34) for in the above expression yields 


Mk, ©;,©r) 


K 2 + K ——— cos 0,cos<t) 
v c i y 


co t- Q>r 

c i 


cos 0,- 


Equation (2.42) is included below for convenience: 


+ m") = /i2 C (A') + |^(o'+^-o)" 

2 a © /2 a© 9 ©" a ©" 2 


The derivatives in (A4) are found from (A3), 


a©; 2 B 2 c\ q 1 , c , 2 JV 2 


(A2) 


. (A3) 


(A4) 


(A5) 


and 
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as2 = _L 

ao>;' 2 B 2 


. — COS 6; -2 
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V'-co'^ 




Cl 


( COS 0; 
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(A6) 


Note that 


a(Bz ) 1 _ 1 a^2 

aco^ S 2 
Using (A5) and (A7), 
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or equivalently, 
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Also, 


d 2 B2 _ 1 
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Evaluating the derivatives at co'= co"= 0 or 


to' = (o'/= co c 


(A10) 


(All) 


yields 


B2c(K) = B 2 {K,G) c + (S)',G) c +(O")\ (o , =m „ =0 = k& 1c P 2 (K, (O c ), 
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B 2 (K,(x) c + co',co c +0)") = B 2c (A r ) 
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(A19) 


Subtracting Eq. (A18) from the conjugate of (A19) yields an approximation for 
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Combining the co' terms in (A20), 
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leads to Eq. (2.45). 


Starting with the to' 2 terms in (A22), define 



(A24) 


Equation (2.47) is obtained by rearranging terms in (A24), and Eq. (2.48) is 
obtained from the cross term, g)"co'. 
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APPENDIX B 

Dispersion of a Gaussian pulse 


2/2 

Beginning with the baseband source term in Eq. (2.73), «(r) = e -r /ts , and the 
definition of g 2 (t) in Eq. (2.42), Eq. (2.68) is derived. Using the well known 
identity 

1 ? c -* 2 /2a 2 e /(ax _ e -a 2 co 2 /2 (B1) 

oV 2k J 

—oo 


in the second integral of Eq. (2.49) yields 

]dc = ke< x -‘ Pa ') 2/ 1 2 , 
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where 

£/(«') = 
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and 

q = -Jt 2 + iP. 

Using (B2), the expression for g 2 (t) becomes 
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(B5) 


where 

P" = P* -i(2P/q') 2 . ( B6 ) 

Using the identify in (B1) to evaluate (B5) yields Eq. (2.75). 
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APPENDIX C 


RMS height: Subtracted Gaussian weighted moving average 


The RMS height of the surface is the surface correlation function evaluated at 


R = 0. 


h 2 = j d 2 K W(K). 

Substituting Eq. (3.66) for the roughness spectrum into (Cl) yields 
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The integral I 2 is solved first. Using integration by parts (let w = e K a / 2 ) t this 


integral simplifies to 
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Next, 73 from (C6) is evaluated using the solution for 72. Let b 2 =2a 2 ; then 
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The first three terms in (CIO) sum to zero in the limit; the last two terms can be 


evaluated using the identity 
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or equivalently, 
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Using Eqs. (C11) and (C12) in (CIO) and (C2) yields 
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h 2 = 2iw 2 a \ . - r(2-Y/2)[2 2 ~' y/2 -1 . 

The RMS surface height is therefore 



(Cl 3) 


(Cl 4) 
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APPENDIX D 

Number and placement of array elements 


Using the intensity front velocity algorithm (Chotiros, 1995), determine the 
number of receivers that are required to obtain a point solution in speed- 
depression angle space. 

Define the reference sensor r 0 to be at the origin 

r 0 =0. (D1) 

The signal arrives at the y'th receiver located at ry = (xj,yj,zjj at time tj. 



a-r, 



(D2) 


where a = a x x + o^y + a z z is the unit magnitude direction vector of the plane 
wave, and c is the propagation speed. Equation (D2) for receivers r 0 ,r 1 ,...,r„ 
can be written as 


~*1 

y\ 

zf 


V 

x 2 

yi 

Z2 



h 

• 

* 

« 

a y 

= c 

• 

♦ 

♦ 

♦ 

_ 


• 

_ x rt 

y n 

V 


Jn_ 


(D3) 


There are n +1 receivers including the reference receiver r 0 . Note that all of the 
receiver locations r!,r 2 ,...r n and the receiver times t h t 2 ,...t„ are assumed to be 

known. The vector a and the propagation speed c' are unknown — three 
unknowns. In order to have a unique solution, the receiver position matrix r is 
required to have rank 3, requiring a minimum of three receivers (r^,^) that 

span 3-dimensional space. Including the reference receiver, a minimum of four 
receivers are required that are not in the same plane. 
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However, Chotiros’ speed-depression angle solution space need only be a 
slice of the 3-D solution space since the plane wave direction has no y 
component (a,, =0). Therefore, the receiver matrix r is required only to have 

rank 2, spanning only the x,z plane. Including the reference receiver, a minimum 
of three receivers (n = 2) are required. 

The following example shows that a linear array results on a curve in the 
speed-depression angle solution space. Assume a four element linear array, 
such that 


ri = (0,0,zi) 

r 2 = (0,0,z 2 ) 
r 3 =( 0 ,0,z 3 ) 


The reference receiver is tq = (0,0,0). Equation (D3) results in 


/ 0 / *2 / h 

a z =c— = c— -c — 
Z] Z 2 Z3 


a 2 x = l 


f A 2 
c'± 

V z lJ 


-a y = l- 


f t ^ 2 
c'4- 

V 4 ) 


(D4) 


(D5) 

(D6) 


because of the assumption a y = 0. A linear array therefore results in a curve 
(2-D solution) in speed-depression angle (velocity) solution space. 
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APPENDIX E 

Calculation of reflection coefficient moments 


The moments of the reflection coefficient (T) are obtained by multiplying Eq. 
(6.8) by the conjugate of itself. To simplify Eq. (6.8), let 


a = T u , 

(El) 

<n£J 

1 

T—H 

11 

(E2) 

c = r 23 e i2k $ 2d , 

(E3) 

so that 


r=a+YCX(-“0"- 

(E4) 


n -0 


Taking the expected value of (E4) multiplied by its conjugate yields 

(m 2 )=(w 2 )+/w 2 x + ( 2 £<-««" }y • < E5 > 

Noting that 

^(C) m )=n% m <E6> 

and summing the geometric series results in 

Equation (6.12) is obtained by substituting (El) - (E3) into (E7). Multiply (E4) by 
(E5) and take the expected value to obtain the moment 
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or 
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+a|y| 2 | r i | 2 


vIyPmV 1 

a (l-|a| 2 r| 2 j 
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1-N 2 ri 2 ' 


(E10) 


Substituting (E7) into (E10) yields 

(Ell) 

Using the relations in (E1)-(E3), the above equation is equivalent to Eq. (6.13). 
Taking the magnitude of (E4) to the fourth power gives an expression for the 


(nrf)=atf - ^Yl 2 - a*Yhr| 2 


moment 
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We use (E7) for the first term on the right side and use the relation (E6) for the 
second, fourth, and sixth terms. The third term evaluates to zero, and the fifth 


term uses the identity 


Re 2 «r 

ln=0 


^ oo 

n =0 


(El 3) 


With much work, the fourth order moment of the reflection coefficient can be 


shown to be 


(|rj 4 ^ = jot | 4 + 4|a| : 
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Using the relations in (E1)-(E3), the above equation is equivalents Eq. (6.14). 
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APPENDIX F 

Ray approximation calculation of reflection coefficient magnitude 


F.1 Linear decrease in squared wave number with depth 


The ray paths for propagation into the lower medium are parabolic. If the 
lower medium were lossless, then all of the energy propagating into the lower 
medium would be redirected upward, and the magnitude of the reflection 

t 

coefficient, r 23 , would be unity. However, loss in the lower medium results in a 
magnitude squared reflection coefficient equal to the energy that reaches the 
surface. The square of the wave number in medium 3 is given by 

k&) = k 2 {\-aQ, (FI) 

where k 2 is the constant wave number in the arbitrarily thin layer below the mean 
surface, and the depth coordinate 

C, = -z (F2) 

increases with depth. By defining the sound-speed gradient as the change in 
sound speed with respect to depth at a point just below the mean surface, a can 
be expressed in terms of the sound-speed gradient g at C, = 0, the speed ratio v, 
and the sound speed in medium 1 , c,. 



(F3) 


Since energy propagating through the sediment will be attenuated by an amount 
determined by the line integral over the ray path, we can approximate the 
expected value of the magnitude squared of the partial reflection coefficient r 23 


as 
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where the path of integration in over the ray path, 0(0 is the grazing angle of the 
ray as it propagates through the lower medium, and k 2 is the wave number just 
below the surface (medium 2). In this example, the mean depth of medium 2, d 2 , 
is allowed to be arbitrarily small. The turning point of the ray, Z max , is found from 
Snell’s law noting that the grazing angle is zero at this point: 

cos(9) _ cos(0 2 ) /pg\ 

< 5(0 c 2 

where 0 2 is the grazing angle of the field just below the mean surface, related to 
the incident field grazing angle 0j, and speed ratio v 2 i is by Snell’s law 


£S5<M = v 21 . (F6) 

cos(0i) 


Evaluating the integral in (F4) yields Eq. (6.15). 


F.2 Linear increase in sound speed velocity with depth 

The sound speed in medium 3 is set equal to the sound speed in medium 2 
at the medium 2 - medium 3 interface; it increases linearly with depth, and is 
given by 

< 5 ( 0=0 +« 0 ^ 2 - ( F7 ) 

The lower medium reflection coefficient follows from the following line integral: 

/ 9 \ f C=Zmax 1 dL 

(r 23 ) = exp -2ImU 2 \ - -- . , 

\ / r=o l + flCsm(0(O) 

By using Snell’s law, the ray turning point, Z max for this constant sound-speed 
gradient profile is easily shown to be 
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max 



1 


a 


V.cos(6 2 ) 




The parameter a can be expressed in terms of the sound-speed gradient and the 
medium 1 sound speed. 


a = 


g 

v 21 c l 


(F10) 


Evaluating the integral in (F8) yields an expression similar to Eq. (6.15) for the 
magnitude squared reflection coefficient: 
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1 + sin(0 2 )V 
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where 



is the sediment gain parameter. 


(F12) 


F.3 Linear increase in sound speed over discontinuous interface 


In the next case, the sediment sound speed is assumed to increase linearly 
down to a depth d 3 , where the sediment parameters are discontinuous. When 

the maximum depth of the ray in medium 3, Z max , is less than the depth of the 

surface separating medium 3 and medium 4, the loss due to the path in the 
sediment above the discontinuous interface is given by (F8). However, when 
Z max > d 3 , the ray intersects and reflects off of the medium 3 - medium 4 

interface, resulting in 
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The above equation is the same as (F8) except the upper limit of the integral is 
d 3 , rather than Z max , and the magnitude of the partial reflection coefficient |r 34 | 2 
is included. Setting a = g/(v 21 q) and evaluating the above integral, we get 



8jc ( l + sin(6 2 ) N ^ 
q ( cos(0 2 ) )j 

8jt f 1 + sin(9 2 ) > 
q l cos(9 2 ) j 
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max 



(F14) 
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